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PREFAŢĂ 


Dinamismul şi complexitatea vieții economice actuale impune. utilizarea 
unor teorii şi mijloace din ce în ce mai perfecționate privind Studiul şi luarea 
deciziilor optime. În analiza riguroasă, cantitativă şi calitativă, a proceselor 
economice, precum şi a evoluţiei lor, un rol important îl joacă metodele 
matematice care prin caracterul lor de generalitate, sunt capabile să modeleze şi 
să optimizeze o importantă clasă de fenomene economice. 

Cursul pe care îl prezentăm încearcă să surprindă acele noțiuni şi metode 
matematice absolut necesare studenţilor economişti, atât pentru modelarea şi 
rezolvarea unor probleme economice des întâlnite în practică, cât şi pentru o mai 
bună înţelegere a noțiunilor economice introduse la alte discipline de specialitate. 

Materialul este structurat pe cinci capitole şi urmăreşte îndeaproape 
programa analitică. În capitolul I sunt prezentate transformările elementare 
aplicate unei matrice şi câteva dintre cele mai importante aplicaţii ale acestora. 
Capitolul al II-lea este consacrat prezentării noţiunilor de algebră liniară şi a 
aplicaţiilor lor, în special în contextul mulțimii R" . În capitolul III sunt prezentate 
problemele de programare liniară şi metodele de rezolvare a acestora. Este tratat 
separat cazul particular al problemelor de transpori. Seriile numerice şi seriile de 
puteri fac subiectul capitolului IV; tot aici se prezintă condiţiile şi modalitatea de 
dezvoltare a funcţiilor în serii de puteri. Ultimul capitol este dedicat funcţiilor cu 
mai multe variabile şi tratează în mod special determinarea extremelor funcţiilor 
de mai multe variabile, atât cazul extremelor libere, cât şi cel al extremelor 


condiționate. 
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Capitolul 1 


TRANSFORMĂRI ELEMENTARE 


Studiul matricelor este legat de discuţia şi rezolvarea sistemelor de ecuații 
algebrice liniare şi constituie elementul de control al aproape întregii programe 
analitice a disciplinei „Matematici aplicate în economie“. Unele noţiuni de calcul 
matricial parcurse în programa de liceu sunt considerate cunoscute (operaţiile cu 
matrice, rangul unei matrice, inversa unei matrice ş.a.). Noutatea pe care o 
propunem se numeşte transformare elementară efectuată asupra liniilor sau 
coloanelor unei matrice, care o duce pe aceasta într-o matrice echivalentă din 
punctul de vedere al rangului. Departe de a fi o abstractizare inutilă, dacă se 
efectuează asupra liniilor poate duce, fără calcule dificile, la determinarea rangului, 


la obţinerea inversei sau la obținerea soluției unui sistem algebric liniar. 


Se cunosc dificultăţile pe care le întâmpinăm la discuţia sau rezolvarea 
acestui tip de sisteme. Metodele bazate pe teoremele lui Cramer, Kroneker-Capelli 
sau Rouch€ sunt greoaie şi folosesc determinanții, al căror calcul este laborios. Prin 
transformarea matricei lărgite a unui astfel de sistem se obține forma explicită din 
care se pot trage concluziile privind tipul şi soluţiile sale. O categorie specială de 


soluţii o formează cele de bază care vor juca un rol deosebit în ceea ce urmează. 


Pentru determinarea rangului unei matrice de tip mXn , dimensiunile m şi n 


pot fi oarecare. Când se caută inversa, m şi n trebuie, bineînţeles, să fie egali. La 
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rezolvarea sistemelor, numărul ecuaţiilor m va fi cel mult egal cu al necunoscutelor 


n (mn) pentru a creşte posibilitatea existenței soluţiilor. 


1.1 MATRICE 


1.1.1 RANGUL UNEI MATRICE 


Prin matrice înțelegem o aplicație A:1xJ-—K, unde I={1,2,...,m}; 


J= {1,2,0}, K o mulţime oarecare. Ea poate fi reprezentată printr-un tablou de 
elemente din K, aşezate pe linii şi coloane. 


Pentru o matrice de tipul m x n vom folosi notația: 


a l a 2 o... d, n i 
A= Mi An ee (an | Ea IN a R 
Amn an2 ir i Amn 


Suprimând din matricea A, m — r linii şi n — r coloane, obținem o matrice 
pătratică de ordinul r. Determinantul acestei matrice se numeşte minor de ordinul r 


al matricei. 


„Definiţia 1.1.1. Numim rangul unei matrice A de tipul m x n (m<n) un 
număr natural r cu proprietăţile: 
„q). în matrice există cel puţin un minor de ordinul r diferit de zero. 


b) toți minorii de ordinul r + 1 sunt nuli. 


Observaţie Se demonstrează că dacă toți minorii de ordinul r+] sunt cul, 


atunci sunt nuli ŞI toți minorii de ordin mai mare ca r + 1. 
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Vom pune în evidență operații pe care le vom numi transformări 
elementare, care se bucură de proprietatea că nu modifică rangul unei matrice» Ele 


au la bază proprietățile determinanților. 


1.1.2 TRANSFORMĂRI ELEMENTARE 


Definiţia 1.1.2. Numim transformare elementară aplicată unei matrice 
una din următoarele operaţii: 
(Tı) — înmulțirea unei linii cu un număr a=0; 


(Tz) — schimbarea a două linii între ele; 


(T3) — adunarea la o linie a alteia înmulțită cu un număr a +0. 


Definiţia 1.1.3. Numim matrice elementară o matrice ce se obține din 


matricea unitate de ordinul m (In) printr-o singură transformare elementară. 


Transformările elementare pot fi realizate cu ajutorul a trei tipuri de matrice 
elementare T; (i = 1,3) înmulțind la stânga matricea A cu aceste matrice. 


Matricea de tipul T, se obţine didie Im prin înmulțirea unei linii cu 
un număr a +0. 

Nini ca de tipul T se obţine din Im prin schimbarea a două linii între ele. 

Matricea de tipul T3 se obţine din Im prin adunarea la o linie a alteia 


înmulțită cu un număr a +0. ` 


Exemplul 1.1.1. Care din următoarele matrice sunt matrice elementare? 


3.0 10| [oa [30 
a) »b) , ©) „ d) 
01 e E0 jo 3 


Matricele a), b), c) sunt matrice elementare deoarece: matricea a) se obține 


din Iz înmulțind prima linie a ei cu 3; matricea b) se obține din Iz adunând la linia a 
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doua prima înmulțită cu 2; matricea c) se obține din Iz schimbând liniile între ele. - 
Matricea d) nu este matrice elementară deoarece se obține din: l» înmulțind linia 


întâia şi a doua cu 3, deci prin două transformări elementare. 


Exemplul 1.1.2. Se consideră matricele: 


y ->12 SSE 
b UL UMEE 340. TI 


Să se determine o matrice elementară cu ajutorul căreia pornind de. la 
matricea A să obținem matricea B. Se observă că matricea B se obține din matricea 
A adunând la linia a doua prima înmulțită cu —2. Acest lucru se poate realiza cu 
ajutorul unei matrice elementare de tipul T2, aceasta obținându-se din I prin 
adunarea la linia a doua a primei linii înmulțită cu —2. Ea va fi: 


i os 
TA = 
Ear ASEE 


1.0 
2 2 ul 


bi 1.4.2 
O Lui durea 


da 


T, = 


aa S9) 
SPRE | | 


Se demonstrează că matricele elementare sunt inversabile şi inversele lor 


sunt tot matrice elementare. 


Teorema 1.1.1. Transformările elementare nu modifică rangul unei matrice. 


Demonstraţie. Fie matricea A de tipul m x n şi matricele T; (i = 1,3) care se 
obțin din matricea unitate Im. Înmulţind la stânga matricea A cu matricele T, 


obținem matricele A; = T;A . Presupunem că rangul matricei A este r}. Vom avea; 
TA, < min (r; TĘ ) = Min (m, TA) 2 TA (1.1.2) 

Din A =TA, cum T; sunt inversabile găsim A = TA, şi deci: . . 
TA, < min (tr Ta) = min (MTA) = TA, si (1.1.3) 


Din (1.1.2) şi (1.1.3) rezultă că T4, = T}, i=1,3. 
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Definiţia 1.1.4. Două matrice A şi B care se obțin una din alta prin 


transformări elementare se numesc echivalente în privința rangului şi scriem A ~B. 


1.1.3. FORMA GAUSS - JORDAN A UNEI MATRICE 


Fie A matricea unui sistem de m ecuaţii algebrice liniare cu n necunoscute 


(m <n). 


Definiţia 1.1.5. Matricea A se spune că are forma Gauss-Jordan dacă 


conţine r (r < m) coloane ale matricei unitate de ordinul m. 


Dacă aceste coloane sunt jy jz, ...; je (<< = <ir, j=1) atunci forma 


matricei este: -: 


(1.1.4) 


Sub linia în scară toate elementele sunt nule. Deasupra liniei în scară pot fi 
elemente diferite de zero doar pe coloanele care nu sunt coloane ale matricei 
unitate, S-ar putea ca ultima treaptă să fie sub matrice, adică r = m după cum 


putem avea j, =m. 


Teorema 1.1.2. Orice matrice nenulă poate fi adusă la forma Gauss- 


Jordan, printr-un număr finit de transformări elementare. 


Demonstraţie. Fie deci matricea A de tipul m x n a unui sistem liniar de m 


ecuaţii cu n necunoscute. 
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Qi Aa. Mj Cn 

do d» aj An 
: Ci aiz ..s aij > sei Ain 

Ami Am2 si: Amn irn ann 


Fie a, = 0. Dacă ay =0, atunci printr-o transformare de tipul (T2) putem 


aduce în locul lui ay un element a, 20; i= 2,m care 'există deoarece prima 
coloană are cel puţin un element diferit de zero, în caz contrar sistemul nu ar avea n 
necunoscute. | 
Printr-o transformare de tipul (Fr) cu a ait realizăm în poziția (1,1) un 
| Gu 
element egal cu unu. Prin m—1 transformări de tipul (T3) în care linia întâi este 
succesiv înmulțită cu numere convenabile: —a,;, —431, igo Şİ adunată de 
fiecare dată la liniile: a doua, a treia, . . . a i-a linie, . .. a m-a linie, realizăm pe 
prima coloană toate elementele nule cu excepţia primului element care este unu. 
Dacă dorim să realizăm pe coloana j toate elementele nule cu excepția 


elementului din poziția (6) ce care să fie egal cu unu procedăm astfel: 


a) ne asigurăm că a; + 0; 


b) aplicăm o transformare de tipul (T.) cu a aN: obținând în poziția (ij) 
: a. 
“ij 


“Sun element egal cu unu; 


c) pa succesiv iA transformări de tipul CI, înmulţind linia i cu: 


=j j Şi adunând- -o la liniile: întâia, a doua, . a m-a 


di 2 ma ADR 
linie obținem pe coloana j toate elementele nule cu excepția elementului 


din poziţia (ij) care devine unu. Algoritmul se opreşte când pe liniile 
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r+l,m nu mai avem elemente nenule. Elementul nenul 'aij dat se va 


numi pivot, iar ansamblul transformărilor elementare necesare 
transformării matricei A într-o matrice cu coloana'i din matricea unitate 


de ordinul m în coloana j se va numi pivorajul cu element pivot dij. 


Din forma Gauss-Jordan a matricei A putem observa că rangul matricei este 
r, acesta fiind de ordinul zeita ră COT cu primele r linii şi coloanele 
Îi» Î35..],. Toţi minorii de ordinul: r+1 sunt nuli. În practică nu este neapărat 
necesar să obținem coloanele matricei unitate în ordinea lor naturală, numărul lor 
ne va da rangul. Putem alege pivot orice element diferit de zero oriunde s-ar găsi el 


de preferinţă cel care este egal cu 1. 


Exemplul 1.1.3. Să se aducă la forma Gauss-Jordan şi să se determine 


rangul matricei: 


2 2 l|  Luăm pivot elementul încercuit. Prima linie 
—l —2 1 1 —2| înmulțită cu 1 se adună la a doua, apoi cu —2 şi 
4 3 


oadunăm la a treia, obținând matricea : 


| 2 2; 2 — 1f =c Luăm pivotelementul încercuit. Linia a doua o 
| 0 0 © 3 al înmulțim cu 1/3 şi apoi înmulțită cu —2, o 
A cra io să iai | 


adunăm la prima, obținând matricea: 


5/3| Luăm pivot elementul încercuit. Linia a treia o 
„ înmulțim cu —l, apoi înmulțită cu =l, o 


0 
NGI 
CD adunăm la a doua şi obţinem: 


Rangul matricei este egal cu 3, 
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Exemplul 1.1.4. Să se determine rangul matricei: 


PĂI i 3 
—-] —2 1 —2 
lwt da val 


Făcând pivotajele cu elementele pivot încercuite obţinem succesiv: 


22 ul leul Pita? aa 
2y e Dat- dg ai INA AJ 
pape” iti, pipa ing Aaa 


5.4 Baste Om Sle D000 (ji Pa a Data 
HAO 0 6|=|1 00 NVU 0 0 
3 i: hajl poni pei) o 1 E 


Rangul matricei considerate este 3. 


1.1.4 INVERSA UNEI MATRICE 
Fie A'o matrice pătrată de ordinul m şi B=[A : 1) matricea obținută prin 


alăturarea matricei unitate de ordinul m la matricea A. Se cunoaşte că se numeşte 


matrice inversă a unei matrice pătrate A, o matrice, notate cu A!, astfel încât: - 
AA =AA=l | | i E1.1.5) 


Presupunem pentru moment că există matricea A. Înmulțind la stânga 


matricea B cu A”! obținem: 
B=A!B=|A ALA 1]=EiA | (1.1.6) 


De la matricea B putem ajunge la matricea B cu ajutorul transformărilor 


elementare. Făcând pivotaje în matricea A în vederea obținerii celor m coloane ale 
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matricei unitate de ordinul m, în dreapta barei vom găsi tocmai matricea A”. Dacă 
coloanele nu sunt în ordinea lor naturală, putem schimba liniile între ele obținând 
în stânga barei chiar matricea unitate. Dacă nu reuşim să obținem coloane ale 


matricei unitate în număr egal cu ordinul matricei A, atunci matricea nu are 


inversă. 


Exemplul 1.1.5. Să se găsească inversa matricei: 


1 "9 13 
A=|—-l —l —2 
a tă 


considerând în matricea B pivotaje cu elemente pivot încercuite avem: 


Cip 2979 EEE TE LE” r i 3 
Bl = —2 0 Vo Sa (Di: 110| 


rar NE -L 0-10 
l —1.0î —2 30 [1 00: -1 -1 1 
e-tc.ro. eg ia atat Bt Aha e 
OAP + mm Mpa > ce asi să dale ada be ek 
ko 0 Sl d. 
~io 10 pP E oniri 
001 CE, ME 


Inversa matricei A este: 
=] -l 1 

Aa 2 „| 

Gai 
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[PI iți 
Aeormanirda_A 
um gh 


Vom avea: 


Dar 3 


B=R 2 LB = 


1.0 3/4 
—5/4 
Dada auf 


[i 


Matricea nu admite inversă. 


PEPI RE e 
20 1 0|=lo (4-5 
20 01| [o 4-5 
1/2 1/4 0 
—12 14 0 

braul 


- Exemplul 1.1.6: Să se găsească inversa matricei: 


100]. 


= 2 1 20|” 


JI 0 


| 
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1.2. SISTEME LINIARE DE ECUAȚII ALGEBRICE 


Un sistem liniar de ecuaţii algebrice cu m ecuaţii şi n necunoscute (m< n) 
se prezintă sub forma: 
Ayx Faa t. F ună = be 
Aaa Halay tot dona Sb (1.2.1) 


Amii ai Am2 X2 Tg Ann Xn ai bu 9 


unde a 


pb; (i=1,m, jÆ ln) sunt numere reale. 


“Dacă există b; = 0 sistemul (1.2.1) este neomogen, în caz contrar este omogen. 


Fie B matricea lărgită a sistemului (1.2.1): 


(122) 


„Dat sistemul (1.2.1) matricea B este unic determinată şi invers dată matricea 
B sistemul (1.2.1) este determinat. Să presupunem că rangul matricei A a 
sistemului este r (r < m) . Prin transformări elementare putem aduce matricea A la 


forma Gauss-Jordan, matricea B căpătând forma B. 


(1.2.3) 
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Sistemul corespunzător matricei B este:. 


Xi NI, Ke anul +a aj i +.. -Fan = b, 
, aTi aie Hajat Înv Hea dop = 5, 


s Haa a +... + dă = b, (1.2.4) 
0 = bi 
0 T bu 


Deoarece matricele B ŞI B sunt echivalente, sistemele (1.2.1) şi (1.2.4) vor 


fi echivalente. Este suficient pentru aceasta să observăm că transformările 


elementare. aplicate liniilor matricei B, pentru a fi adusă la forma B au drept 


consecințe asupra sistemului de ecuații: 
(Tı) — înmulţirea unei ecuaţii cu un număr diferit de zero; 
(T2) — schimbarea a două ecuații tre ele; 
(T3) — adunarea la o ecuaţie a alteia înmulțită cu un număr diferit de zero. 
„Acestea reprezintă transformări elementare aplicate unui sistem de ecuaţii şi 


au evident proprietatea că îl transformă în . unul echivalent, Transformările 


elementare nu se vor aplica asupra coloanelor matricei lărgite deoarece nu s-ar 


obține un sistem de ecuaţii echivalent. 


Teorema 1.2.1. Condiţia necesară şi suficientă ca sistemul (1.2.1) să fie 


compatibil este ca b; =0, i=r+lm. 


Demonstrație. Să arătăm mai întâi necesitatea condiţiei. 


Presupunem că sistemul (1.2.1) este compatibil. Aceasta înseamnă că 
sistemul (1.2.4) este compatibil. Din compatibilitatea acestuia din urmă rezultă 


b; = 0 pentru i=r+l,m. 
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Să arătăm acum suficiența condiției. Presupunem că b, =0, i=r+l,m. 
Rezultă atunci că sistemul (1.2.4) este compatibil admițând o" ” soluţii şi deci 
sistemul (1.2.1) este compatibil. 

In baza acestei teoreme uşor pot fi demonstrate teoremele lui 
Kronnecker-Capelli şi Rouchâ. Pentru r = m = n se obţin rezultatele cunoscute de 


la sistemele de n ecuaţii cu n necunoscute. 


“Exemplul 1.2.1. Să se studieze “compatibilitatea sistemului Şi în caz 


afirmativ să se găsească soluţia lui. 


— 


X + Mda Ay = 
| A Mă 4209 2%, | E 
3 x F8 —3xj = 5 
— dy F dz Sa dag = =3 


5 


Matricea lărgită a sistemului, după o serie de pivotaje cu elemente pivot 


încercuite, devine: 


O 1 -1 1M 


Sistemul considerat este echivalent cu sistemul: 


a H2 1a = 372 
salisi LIY 2 = lD 


care este compatibil având soluția: 
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x, = 342.— MM 2ak4/24, 
x, = —1/2+3/2a—3/26, X; =Q, -X4 =b; a, b ER 


Exemplul 1.2.2. Să se determine compatibilitatea sistemului: 
=w A X + X3 ze X4 = 3 
Matricea lărgită a sistemului după o serie de pivotaje cu elemente pivot 


încercuite devine: 


[Da 3 aja [ui 2 PIL A h 0s Wapa mAN 
B=|-1 1 1 —1| 3G 4 3] S|-|oLI 4/3 1 
—3 5.6 110| -10 9 sol af-l0stor a 0-04] —1 


Cum b.==—lz 0, sistemul este incompatibil. 
3 p 


1.2.1. EXPLICITAREA SISTEMELOR LINIARE 
DE ECUAȚII ALGEBRICE s 


Fie sistemul (1.2.1). Presupunem că rangul matricei A este m < n. 


Definiția 1.2.1. Vom spune că sistemul (1.2.1) este explicitat în raport cu 
un grup de m variabile (necunoscute), dacă matricea A a sistemului conține m 


coloane ale matricei unitate de ordinul m. 


Deoarece cu cele n variabile ale sistemului putem forma CẸ grupuri diferite de 
câte m variabile, rezultă că un sistem liniar, care poate fi explicitat în raport cu cel 


puţin un grup de m variabile, va avea cel mult CẸ forme explicite. Acest număr 


maxim va fi atins dacă explicitarea poate avea loc în raport cu oricare dintre cele C 
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grupuri diferite de câte n variabile şi toate formele explicite sunt distincte. 

Într-un sistem liniar explicitat sunt două tipuri de variabile (necunoscute): 
unele,. ai căror coeficienţi formează coloanele matricei unitate, altele..ai căror 
coeficienţi formează celelalte coloane ale matricei sistemului. Ca de obicei, primele 
vor fi denumite variabile principale, iar celelalte variabile secundare. În teoria 
programării liniare variabilele principale mai sunt denumite bazice: iar variabilele 


secundare nebazice. 


Definiţia 1.2.2. Numim soluţie de bază a unui sistem liniar de ecuații 
algebrice, o soluţie a sa, obținută prin egalarea cu zero a variabilelor secundare şi 


atribuind celor principale valorile ce rezultă dintr-o formă explicită. 


Din definiția dată rezultă că dacă sistemul are cel puţin o soluţie de bază 
atunci va avea cel mult CP soluţii de bază. 

Dacă într-o soluţie de bază avem exact m variabile diferite de zero atunci 
spunem că soluţia de bază este nedegenerată, iar dacă are mai puţin de m variabile 
„diferite de zero spunem că este soluţie de bază degenerată. 

Soluţiile de bază ale unui sistem liniar se clasifică şi după un alt criteriu. Dacă 
orice variabilă are într-o soluţie de bază ò valoare nenegativă, vom spune că soluţia de 


bază este admisibilă. În caz contrar, vom spune că soluţia de bază nu este admisibilă. 


Să considerăm că r = m. În acest caz din (1.2.4) obținem o soluţie de bază: 
restul până la n egale cu zero. Dorim ca variabila principală x, să devină 
n Ji 
. . . x s y . 8 n 90, $ ingiž Ja 
secundară iar variabila secundară xJ = js ---,) să devină principală. Pentru 
aceasta fie şi a; = 0. Efectuând în matricea B un pivotaj cu element pivot a; 


găsim o altă soluţie de bază: 
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Şi Au pane la n rapa cu zero. 


“Putem scrie variabilele principale concentr at sub forma: 
y E b; > . 
= b, = i dj => pentru k = J 
ij y 
ajj = 0 


x =h , pariu E=j 


sa 


ahli2) 


Exemplul 1.2.3. Să se determine toate S explicite ş şi ap opri de bază 


corespunzătoare pentru sistemul: 


Xi o ka =0, 


Matricea lărgită a sistemului este: 


lun ile zale a Oulu 
Bu | 
zi Adria CI 
Pivotajele: i ~ | | E 
A W lsb nO Pi dloneaalenlute0 d Oct uide 0 
PES. 12] |0 01 
Oi a aţa. = A 
b) 
2 1|2 aur 1 0 1]2] 
Sale 2123 zi Clol | 1p ta 
c) laS) mpe TE 
-INET TT A an A A O Pa 
a Ni 71 0a T =E. POO [i T 1 010 
A amia ud 2] o A OANA] Joao 04 2 
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N 
n 


l 1 00; |—-1 1 —1 0[0| |l-i 1 -1 010 
la nf la a 10 1l 2 h œ 3 "Ura 
conduc la formele explicite în raport cu grupurile: 
(x „A ), (x, X4 ) (x> X3 i (x sX), (33 3) 
| o PRE =2 qi ti =2 
i | r uE? i ý ac) +y = 2 
“> + x, = x +X, — X3 =0 
oha Pat, 2 ad Xp 


e) X] => Na + A = 0 
X + 2%, + Na = J 


Soluțiile de bază vor fi: 

(S) A =2, x =2, x =0, x, =0; soluţie de bază nedegenerată admisibilă 
(52) x1=2, x =0, x =0, x = 2; soluție de bază nedegenerată admisibilă 
(53) x50, x% =2, x= 2, x= 0; soluţie de bază nedegenerată admisibilă 

RA A 3 4 
(S4) x, =0, x =0, x, =0, xy = 2; soluție de bază degenerată admisibilă 
(S) x, =0, x =0, x= 0, x, =2; soluţie de bază degenerată admisibilă 

5 l 3 4 8 


Sistemul nu poate fi explicitat în raport cu grupul de variabile (x,,x) şi 


admite doar 4 soluţii de bază distincte. 


pin, rate) OURS 1- „Ér 

A T IORN h aR 

şi testul până ln regale cu eră. :sihueţirg dă óga ni stidilgxo S si pubio 
Pule în ienie variubilefe prunci Pb EI ANALE) e q$ 


BY dr baad ab alt 
b m sI 


TERT El ig ath XE. a TE Svi x il 

BI țozaiy ubh Bono Dear ab gi = a Osa E= is Eia (2) 
| n da iin pipe N iei ala 21 imi) su 
Xe pi catei isa (+8) 


| Arti > pe 


Ar aia ral AN 


4 j w p J | T | E a 
r 4 w Lă 
h Li 
4 a E F 
m: Ga, Mie 
| J 
fà a 


y x} ie pee eA E ENS 


k N i >, i ji E it PTa oera 
AI i du | SiN har RE wT a 
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G i | | Di 
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SPAȚII LINIARE 


Inzestrând o mulțime cu operații având anumite proprietăți obținem o 
structură algebrică. Unele structuri au fost studiate în liceu (grupul, inelul, 
corpul). Spaţiul liniar se defineşte în. acelaşi mod, elementele sale (vectorii) 


generalizând noţiunea = vector liber utilizat în fizică sau matematică, 
În orice spațiu liniar se introduce noțiunea ia actului liniară, precum 
şi cea de bază. Aceste consideraţii se fac şi aici, dar se reduc la spațiul liniar R” 
Funcţiile de un ini tip care se stabilesc între diferite spații liniare se 
numesc operatori. În cele ce vom prezenta ne ocupăm cu operatorii liniari care au 


numeroase aplicaţii în diverse domenii de activitate umană. 


In studiul problemelor economice se impune problema determinării valorii 
maxime sau a valorii minime a unei funcţii. Un ajutor prețios în acest scop îl aduc 


formele pătratice care generalizează funcția de gradul al doilea. 


2.1. DEFINIŢIA SĂ Ra YU! tii 


Fie A o mulțime oarecare şi K un corp având elementele” neutre Og, 


respectiv Ik. Presupunem că pe A se pot defini două operații astfel: 
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a) (V)x,y EA, (3)un element notat cu x+y€A, astfel încât 
(xy) x+y. 
b) (V)xeA, (V) XEK, (3) un element notat AxEA, astfel încât 
(XX) Ax. 

Prima operaţie, de tip aditiv, este internă, a doua, de tip multiplicativ, este 
externă pentru A. 

Definiţie. Mulțimea A formează un spaţiu liniar peste corpul K dacă sunt 
satisfăcute următoarele proprietăţi: | | 

Li) (V)X, y, ZEA, XEO HZS X EYy)+ Z, 

L2) (V) X, YEA, U xp yEy, 

La) (JO, EA astfel încât x+0, =x, (V) XEA, 

La) z € re 3) XE A astfel încât eri (—x)=0,, 

L) (V) NEK, (V) Xy EA Ay, 

Le) (Y) AEK, (v)pueK,.(V)xeA, (XDA AX Hya 

L7) (Vu E K, (XE A, ~ NUX = pu) ALX, 

Ls) (V) E A e Spa 

Elementele „lui K se numesc scalari, iar cele ale lui A, vectori. Prima 
operaţie se numeşte adunarea vectorilor, iar cea de a doua, înmulţirea cu scalari din 
K. Atunci când nu este pericol de confuzie, vectorul nul Q4 se notează simplu 0. De 
asemenea, elementele neutre ale lui K se scriu 0 şi L. 

i Consecințe. Din axiome decurg următoarele proprietăţi: | 


a) (V)xeA, 0x=0, 
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Într-adevăr, din Le) obținem- (A +0)x =)x+0x, adică Ax =\)x+0x. 
Axioma L; ) conduce la 0x =0. | 

b) (V)NEK, X.0=0. | 

Plecăm de la Ls): X(x+ 0) = = AO. Obţinem Ax = " + AO. 
Aceeaşi axiomă L3) arată că: AO 0. 

c) (V)xEA, (-Dx=-—x 

Folosim a), adică 0x = 0 şi A dă ti O prin 1+(—1). Din Le) obţinem: 


Ix + (—L)x = 0. Lg) conduce la x+ (— Dx = 0, iar Lig) arată că (—1)x i p 


Exemple: 
1) Considerăm A= R" = RxR R= (x = (Appel X ER) şi 


K=R. Adunarea elementelor din A revine la adunarea matricelor cu o 


-singură coloană sau cu o singură linie. 
(VIX = (ep not) YO Ya) din A, 
XEY = (x Yo at Yare Xa EA 


In mod similar produsul unui scalar din R cu un element al lui A se 


defineşte prin: 
AX — Xn) = Ax, AX, ÎN, Axa) 


Verificarea celor opt proprietăți este imediată. 


2) Fie sistemul omogen MX =0, în care M EMnn(R) XEMa(R), 
0E Mni(R). Mulțimea S a soluțiilor acestui sistem formează un 


spațiu liniar peste corpul K=R. 
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Într-adevăr, SD == (Xe Mu (R) MA = 0) permite definirea sumei a două 
soluții X’, X” şi produsul unei aiii X’ cu numărul Teal A 
(V) X',X" ES >= X'+ X” es, 
| VAER, VX’ ES re s 

Într-adevăr, X’, X” € M (R) implică X'+ x” EM, (R) şi 

M(X' +4+X")= :MX +MX' =0+0= 0 | 
phi X eM,, (R), A ER implică XEM (R) $ i 

MOX') = MĂ =0=0. 


Cele opt seci dia sunt verificate imediat, aarti Os find soluția 


Bankia ali stiai, 


Majoritatea aplicațiilor care “urmează sunt realizate în R” “privit ca o 


mulțime de matrice cu o'singură coloană, R” = M (R). 
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2.2. DEPENDENȚĂ ŞI INDEPENDENȚĂ LINIARĂ 


Fie A un spaţiu liniar peste corpul K şi fie x, e A, €K, i=lh. 
„Definiţie 2.2.1. Se numeşte combinaţie liniară a vectorilor x, cu scalarii 
A; expresia: 


E=AĂ, + Xa Heret AX ne | (2.2.1) 


Consecința a) din definiția spațiului liniar ne arată că atunci când 


N =0,(V)i=lh = E=0. Reciproca acesteia nu este întotdeauna adevărată. 


Definiţie 2.2.2. Vectorii x; = A, i=], h se numesc liniar independenți 


dacă E= At AX ke FAX 30 Ni = Az == =i., | (dd) 


In caz contrar, dacă E = 0 cu măcar un scalar nenul, vectorii x; se numesc 


liniari dependenți. 


Exemplu: 


Fie spaţiul liniar A = R” = M,,(R) şi vectorii: 


neie E E T. Ae 3 
X1 = (Xp Bigat- Xin) Xa = (Xap, X22511 Xa) scăpa Xa Xp) 
vectori din A. Condiția E = 0 conduce la sistemul liniar şi omogen: 


X, PR i, e n, ENA Xa = (2.2.3) 


Xi Ñin rič Aa Xon wiria Ah Xin = =0 
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Vectorii sunt liniari independenţi dacă sistemul (2.2.3) admite numai soluţie 
banală şi liniar dependenți în caz contrar. 
Din teoria sistemelor omogene independenţa liniară este caracterizată prin: 
Popi ir (2.2.4) 


unde matricea ME M,„(R) este formată prin elementele vectorilor X; 


(componentele vectorilor x; ). 


Dacă (2.2.3) are soluții nebanale, ele se exprimă prin h—r (M) parametri 


reali. Înlocuite în (2.2.2) aceste soluţii conduc la una sau mai multe relații de 


dependenţă liniară. 

Pe baza acestui exemplu vom realiza un mare număr de aplicaţii. 
Proprietăţile familiilor de vectori liniar independenți şi al familiilor de vectori liniar 
dependenți sunt date în următoarele teoreme: 


Teorema 2.2.1. O: familie de vectori liniar independenţi nu poate 


conţine vectorul nul. 


Demonstraţie: Presupunem prin absurd că familia de vectori independenți 


F E, A ..>X, + ar conține vectorul nul si că x; = 0. Combinația liniară : 
PS h 


E=0:x HOX. kd: XA = 0a cua scalarul > nenul. 


Contradicția justifică teorema. 


Teorema. 2.2.2. Orice subfamilie a, unei . familii. de vectori liniar 
independenţi este formată din vectori liniar independenți. 
Demonstraţie: Tot prin reducere la absurd presupunem că deşi Xp Xose Xy 


sunt liniar independenţi, subfamilia (X; Xpo Xi J} k <h. Este formată din vectori 
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liniar dependenți. Combinația liniară KI BRALA LAU PX, =0 cu măcar un 
| scalar nenul. 

Rezultă că: 

AXi PAR Xi, FAX + PLR =0 


jEI 


cu măcar un scalar nenul ceea ce contrazice ipoteza teoremei. 


Teorema 2.2.3. Într-o familie de vectori liniar dependenţi cel puțin un 
vector se serie ca o combinaţie liniară a celorlalți. 

Demonstraţie: Definiţia familiilor de vectori liniar dependenţi arată că: 

AXi + AX te HAX t+ AX = 0 


cu măcar un scalar nenul. Presupunem că acesta este \;. Cum orice spațiu liniar 
este grup abelian, regulile de calcul dintr-o astfel de structură permite să scriem: 


NIX; = NIX, SAN > SNAN a Xp = NXg 
In corpul K, X, este inversabil iar dacă X; | este inversul său 
ae m X : zi > miN g 
X= AA AR DDA XI Ai AT Aaa (2.2.5) 


şi teorema este demonstrată. 


Când corpul de scalari este R” inversul >; ! al lui X, este 1/X; şi obţinem: 
A; A 

tss Xp Xa e aa aana (2.2.6) 

Teoremele prezentate au consecințe importante cu valoare teoretică şi practică. 


Astfel, modificarea componenței unei familii de vectori liniar independenţi prin 


suprimarea unora conduc totdeauna la o familie de vectori liniar independenţi. 
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Dacă adăugăm noi vectori unei astfel de familii, rezultatul nu este previzibil. Din 


(2.2.3) rezultă că (2.2.4) nu poate fi realizată dacă h>n, deci numărul maxim de 
vectori liniar independenţi din R” este n. 
Există efectiv mulțimi de vectori liniar independenţi. Una dintre acestea este: 


F= (epezs.:e,) | | (2.2.7) 


t i : t NNA iin 

unde: e, = (1,0,0,...,0), e, = (0,1,0,.:.,0),...,e, =(0,0,0,...,1Y. * 
Dacă prin transformări elementare modificăm matricea M a vectorilor 
familiei F, obținem alte familii formate tot din n vectori. liniar independenţi. 


Numărul acestor familii este infinit. 
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2.3. BAZE DE VECTORI. a pată X 
COORDONATELE UNUI VECTOR ÎNTR-O BAZĂ 


Am luat în studiu în paragraful precedent familiile maximale de vectori 


liniar independenți. Rolul acestora va fi precizat în cele ce urmează. 


Definiţie. Familia B = (ua) formează o bază în spaţiul vectorial 
V dacă: 


a) u,up,...,uy Sunt liniar independenți 


b) (V)ueV, familia F= (uu uzi), este formată din vectori 


liniar dependenţi. 


Familia B este aşadar o familie maximală de vectori liniar independenţi. 
Când spaţiul V este R”, această familie conţine n vectori, deci k = n. 


numită bază 


n 


Spațiul R” are cel puțin o bază: cea formată din e,„e,,...ue 
canonică. Aşadar: 
B. = GEZAK 


Din aceasta prin procedeul indicat anterior putem obține o infinitate de baze 
toate având n vectori. Numărul n, al vectorilor oricărei baze din R” se numeşte 
dimensiunea acestui spaţiu. 

Condiţia b) din definiția bazei şi una din proprietățile fundamentale ale unei 


familii de vectori liniar dependenţi ne permit să formulăm şi să demonstrăm: 


“ Teorema 2.3.1. Orice vector din R” se exprimă în mod unic printr-o 


combinaţie liniară a vectorilor unei baze B. ; 
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Demonstraţie. Teorema 2.2.3 ne arată că cel puţin un vector al familiei 
{u uu,- 0} în care (uz...) formează bază în R” se exprimă 


printr-o combinaţie liniară a celorlalţi. Acest lucru este valabil şi pentru u, deoarece 


în caz contrar, din: 
Au +A u AU, koku, =0 6) (2.3.1) 
în care X = 0 ar rezulta: 
| CA «n A e =0 (2.3.2) 
şi unul din scalari fiind nenul s-ar contrazice condiţia a) din definiţia bazei: 


Aşadar A =0 şi obținem 


dr N 
U= ——u, ——u, — u RIA 
aT Ty i 
A A A 
Notăm py = ==, u= E. u, sj scriem 
Hi òg H3 y Hn e Ş 
U = pt, pia FF Hate | (2.3.4) 


Prin reducere la absurd demonstrăm că scalarii p; din (2.3.4) sunt unic 


determinaţi. În caz contrar are loc relația: 
u= Vu, tau e HvU (2.35) 
“(vu PANS Vu e (a 2 du 20 (2.3.6) 
Cum u Ud, u, sunt liniar independenți, toți scalarii din (2.3.6) sunt nuli, deci: 
ra VAa = Varu esa = Va 


ceea ce este în contradicție cu ipoteza. Teorema este demonstrată. 


Observaţie. Scalarii și, (i= 1,n)} se numesc în mod uzual coordonatele 
vectorului u în baza B = (u uzi): Unicitatea lor rezultă şi din teoremă, dar 


şi din faptul că (2.3.4) conduce la un sistem Cramer (cu soluție unică). 
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Vom nota mulțimea coordonatelor vectorului u în baza B cu u g» adică: 
| i 
Ug = (His Has.. Hn) | 2-3.7) 
Un calcul elementar ne arată că în baza canonică B, = fe,,6y,....e,) 
coordonatele oricărui vector u coincid cu componentele sale, adică: 


| 
U= (Up Uau) ue, true, +-:-+uie, 


Exemplu: 
Să se arate că vectorii u, = (1,1,1)', u, = (2,0%, u, = (1,1,3)! formează o 
bază în R? şi să se determine coordonatele lui u = (6,1,8)! în zeui bază. 
Relaţia (2.3.4) devine: AA, 
U = piu + pU, + pat 
şi conduce la sistemul: 
ÎN hat u3=6 
Hi +2 + pa = 
Hit Hz +3p3 =8 
Prin transformări elementare de linii în matricea lărgită a sistemului obținem: 


1 F-1j6 Ailey kang 1): O dell 1 0 0| 10 
1 2 1| 14071 TISNO 1'01—si-l0 i 0-5 
1 7 31) PU-U a 2] |U 0 a 27 109] l 


care corespunde sistemului cu r(M)=3 ) deci (u4, uz, u; formează baza). 


Hy =10 
u> = —5 
p3% Í 


acestea fiind valorile coordonatelor căutate. 
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2.4. SCHIMBAREA COORDONA TELOR UNUI 
VECTOR LA SCHIMBAREA BAZEI 


Dacă în iz bazei B= era „u şi, a sk liniar R” Panoorin o 
altă bază B' = (vu pa va} având aa număr de vectori, n, ne aşteptăm ca în 


relaţia (2.3.4) coordonatele să se schimbe. A să aflăm legea de schimbare a 
acestora în cazul general şi într-un caz particular. 
Să considerăm că în baza B (numită baza veche) vectorul u are 


coordonatele ci S 1, Gel deci: 


E „A Ne ASP 
şi că acelaşi vector are în baza B’ (numită baza nouă) coordonatele Bs HE , adică 


Pe de altă parte fiecare din vectorii v; ai bazei noi are în. baza veche 


coordonatele a., unic determinate 


ji? 
Vi =at + azt, H.. Fa, = Sau i (2.4.3) 


Înlocuim (2.43) în (2.4.2) şi comparăm cu (2.4.1). 
Obținem: 


A > ai Haba Het aab i=ln (2.44) 


Dacă notăm cu A matricea: 


SPAȚII LINIARE 39 


numită matrice de trecere, cu up matricea (aaaea) şi cu up: matricea 
(B,»Bas++--8,.) atunci (2.2.4) se scrie matricial : 


Up = A'uy (2.4.5) 


Cum vectorii lui B’ sunt liniar independenți, A este inversabilă şi 
t zl 5 


Am dedus astfel legea de transformare a coordonatelor unui vector atunci 


când se schimbă baza. 


Exemplu: 
Se consideră vectorii u, = (1,2,1)', u, =(2,1,0)!, u, = (1,1,0)! din R*: 
a) Să se arate că Reg? o bază în R? 
b) Să se determine coordonatele lui u = ( 1,1,—2)' în această bază. 
c) Să se determine coordonatele lui u în baza formată din 
IV E Up HU, Va =U i =U, V3 = uyt 2u 
Soluție: 
a) Matricea formată cu componentele vectorilor U,Uu>,u are rangul trei, 
egal cu dimensiunea lui R?, aşadar cei trei vectori formează o bază. 


= b) Vectorul u se scrie u=a,u, +a,u, +a,u, care conduce la sistemul 


AT 2A F aS 
W FAE a =l 


cu soluția: a, = —2, a, =—2, a =7,deci: 


up =(-2,—2,7),B = {u,,u,,u;} 
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c) Matricea de trecere de la B la B' = fvi, Va, v3} este : 


—] 1 0 
Ax |] =] 
01 $: 
iar inversa A! 
—3/5 2/5 —1/5 


(A'y!=| 2/5 2/5 -1/5 
"1/5 1/5 2/5 


Obţinem în final: 


Up = (— 1,—3, 2)! 


Observaţie: Dacă vectorii noii baze nu sunt exprimaţi de la început prin 
combinaţii ale vectorilor vechii baze, coordonatele lor (deci elementele matricei de 
trecere A) trebuie calculate separat. În această situaţie este mai simplu să calculăm 


direct up, folosind definiţia. 


În unele cazuri bazele B şi: B’ diferă printr-un singur vector. Formula de 
calcul a noilor coordonate va avea evident o formă mai simplă. Să presupunem că 


în baza B, vectorul u; s-a înlocuit printr-un vector v; care are coordonatele : 


Aj a jo see dijo eo anj adică 


Yj z aju; + aju teret aju; Ti: Qnjila i & (2.4.7) 


Ceilalţi vectori ai vechii baze rămânând neschimbaţi, matricea de trecere A 


devine: 
| 0 0. 0 
0 l 0. 0 
ajj aj aij aij 
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După determinarea lui (Ah)! şi aplicarea formulei (2.4.6) deducem noile 


coordonate: 


Formulele (2.4.7), greu de memorat pot fi utilizate prin intermediul unui 
alzoritm de calcul bazat pe transformări elementare ale liniilor matricei 
coordonatelor vectorilor u,u,,u,,...,u;,...u, şi vj. Ilustrăm mai jos modul de 


organizare a algoritmului. Etapa iniţială are forma : 


Coloanele tabelului conțin coordonatele fiecărui vector în baza de pe prima 
ETTER Ele se presupun cunoscute. pia coloană reprezintă coordonatele lui u 


în baza veche B. 
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Etapa finală are forma: 


şi conţine oboi ala tuturor vectorilor în noua bază. 

Se observă că linia i din ultima fază se obține prin împărţirea liniei i din 
prima fază prin a; (care este nenul, altfel matricea de trecere A nu este 
nesingulară). 

O linie oarecare k a matricei finale se obține înmulțind linia i cu —a,j după 


care se adună cu linia k a matricei iniţiale. Ca efect al acestor operațiuni, pe ultima 


coloană se obțin coordonatele f;,3,(k zi), pe coloana v; se obține coloana „i“ a 


matricei unitate I, iar coloana lui v; se modifică în mod corespunzător. 


Algoritmul poate fi aplicat înlocuind cauzele cu efectele şi invers. Astfel, 
„putem să aplicăm transformările elementare de linii în matricea fazei iniţiale pentru 


a transforma coloana lui v; în coloana vectorului u; din faza iniţială. Unii numesc 
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această operaţie „Lema substituției“ şi aplică diverse metode mnemotehnice de 
utilizare. Noi o considerăm o simplă aplicare a transformărilor elementare ŞI o 


folosim succesiv înlocuind toți vectorii bazei iniţiale B prin vectorii bazei B}. 


Exemplu. În baza canonică B.= (e exe.) se consideră vectorul 
u = (1,1,—1)'.:Să se determine coordonatele sale în baza B'= {vis V2>V3) unde 


vı =(1,2,3)', v, = (1,—1,1)', v, =(0,0,D'=e,. 


Baza e  & e vi v2 u 
d e e Pa m Gl pa 

e, 0 l 2 -l l 

e; 0 0 l 3 l -l 

V, l 0 0 l l l 

e, —2 1l 0 0 -3 -1 

e; e ACUL l 0 -2 4 


vi 1/3 NN 0 l 0 2/3 
Va WA =1B3 0 0 l 1/3 
e; = -23-1 0 0 —10/3 


Operaţiile — extrem de simple — au urmărit crearea matricei unitate pe 


coloanele vectorilor noii baze B' = {v}, Vase]. 


În prima etapă am completat tabelul cu coordonatele - coincizând cu 
componentele = tuturor vectorilor în baza canonică B.. S-au aplicat transformări 
elementare de linii pentru a obține pe coloana lui v; prima coloană a matricei 
unitate LI. În ultima etapă, prin transformări elementare de linii se obțin pe coloana 
lui vz cea de a doua coloană a matricei unitate 1. Coordonatele lui u în noua bază 


sunt înscrise pe ultima coloană şi deci: 
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pkg igt “iiy 
e adie $: 
Observație: Chiar dacă pare prea complicat în comparaţie cu aplicarea 
definiţiei în determinarea coordonatelor lui u, algoritmul are avantajul principal că 
în fiecare etapă coloanele conţin coordonatele în baza curentă. În plus, ca avantaj 
secundar este posibilitatea de i-l folosi în probleme de programare liniară. El va fi 


aplicat consecvent în capitolele care urmează. 
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2,5. OPERATORI LINIARI 


2.5.1. DEFINIȚIE. GENERALITĂȚI 


Considerăm două spații liniare R” şi R” şi o aplicație L'R" R’, 


Definiție. Aplicația L se numeşte operator liniar dd 
La) L(x tX )= L(x tL) (V) xx e R” 
b) L(kx)=kL(x) (Y)k E R, (Vx e p™ 
Prima proprietate se numeşte aditivitate,iar a doua omogenitate. Un 
operator liniar estè aşadar o funcție aditivă şi omogenă între două spații liniare. 
Din definiție decurg imediat două proprietăți esențiale ale unui operator liniar: 
a) dacă On» 0, sunt elementele nenule ale lui R™, respectiv R”, atunci 
L(a J= 07. 
b) Dacă — x este elementul pije (opus) lui X în R", atunci 
Ax) > LX) 


Mai observăm că proprietăţile a) şi b) din definiţia operatorului liniar se pot 


condensa într-o proprietate mai generală: 


c) L(ox,+6x2)=aL(x,)+BL(x2), (V)oBeR, (V) xx, ER” 


Echivalenţa dintre c) cu a) se poate demonstra imediat. Proprietatea c) se poate 
generaliza prin inducție matematică la o combinaţie liniară de k vectori din R™: 


pentru (Y) a ER, (Y)x ERT, i= l,k. 
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Exemplu: Să se arate că aplicaţia L: R° — R?, definită prin: 
Pa Herse A Malot SO fe, pp pa n... 
L(x) = bn] ja (ap t a tta a) 
este un operator liniar. 
Soluție: Aplicăm proprietatea c) care conduce la: 


Ek BEDS L aali in 43) FB Xa) )= 


=a(X1 +X + 2x 2X — x2) + Bla + X A2 2X1 — Xa) = aL(x)+BL(x) 
Observație: Caracterul de operator liniar se poate deduce din calitățile de, 
funcţie liniară şi omogenă în X}, X25... Xa vectorului L(x) e R". 


2.52. NUCLEUL ŞI IMAGINEA UNUI OPERATOR LINIAR 


Să considerăm un operator L:R" — R". 


Definiţie: Se numeşte nucleul operatorului liniar L mulţimea Ker LCR” 


definită prin : .. 


Keri = {x ER" La) Amen HUE RANA 
Se observă că Ker L nu este vidă, cel puțin 0, fiind un element al său. Mai 


mult, fără dificultate se poate demonstra că Ker T este'un subspațiu liniar al lui R™ . 


Definiţie: Se numeşte imaginea operatorului liniar L mulțimea Im L C R” 


definită prin: ImL=fy ER" / IxER” aî. L(x)=y} (2.5.2) 
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„Nici ImL nu este mulțime vidă, cel puţin 0, fiind element al său. Și aici se 


poate demonstra că ImL reprezintă un subspaţiu liniar al lui R”, 


Exemplu: Se consideră operatorul liniar L:R?— R°- definit prin 
xi oaie a e Pe 
L(x) = ATEEN je (2x1 =X Xa Ay Haa = 23 ka). 
Să se determine nucleul şi imaginea operatorului L. 


Soluţie: Se observă mai întâi că toate componentele vectorului L(x) sunt 


funcţii liniare şi omogene de variabilele x,,x>,x deci operatorul este liniar. 
Pentru determinarea nucleului trebuie rezolvat sistemul 
L(x)=% 
adică: 
2x] == Xa +X = 0 
Xi + Xa = 0 
-2 + xy = 0 
care are o infinitate de soluţii exprimate prin: 
A | 
Wa, m=R Wa, AER 
1 2 3 
3 3 
„Un vector oarecare x = (x1, x3,X3) al nucleului Ker L are aşadar forma: 
t t 
l i l 1l 
A ae A A A =a]——,=,ll,acR 
3 3 Eur) 
In acest caz, nucleul Ker L are dimensiunea 1, deoarece orice vector din 


t 
Ker L se exprimă numai funcție de v = ia) 
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t 
Ker L=4x €R? safira] wer (2.5.3) 


Determinarea imaginii operatorului, Im L, revine la rezolvarea sistemului: 

i L(x) =y, 

De fapt căutăm valorile parametrilor y); y2, y3 pentru care sistemul admite 
soluţii. Plecând de la matricea lărgită a sistemului şi aplicând transformări 


elementare de linii obținem succesiv: 


2 a E Te EA s AE elit 
M=|1. 1-0 Ya |n eirda O] Mr 
] —2lu'bpiypzpie asa NR 
canti ati > (ră 0 ý; 
H gan Dă oa e a, a 3 e A BEI e i 
0 —3 1 —Y2 h Y3 0 —3 l —Y2 RN Y3 | 


10 1/3 1433, +1/3y, 
O o —173 TYN 27 3y 
0 0 0 Wit YaF Y3 


Sistemul corespunzător acestei matrice este compatibil dacă şi faimă] dacă: 
=I T P F Yat | 
ceea ce înseamnă că y = 04 y; =B iariy =a +8, a; BER: Um'vector oarecare 
y = (Yp Y2, y3) din Im L se scrie sub forma: 
y=464B,a4+f)!=a(1,0,1)! 80,1, a, BER. 
 Coønchidem că ImL ca spaţiu lar are Al Ape egală cu 2 şi că o bază 


în acest spațiu este formată din: . 


XI (1/0; 1)t, #5 (0,141)! oala 
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Observaţie: Se poate arăta că pentru un operator! L: R™ = R”, are loc 
proprietatea: 
dim KerL + dim ImL=m (2.5.4) 


care în exemplul precedent este evident verificată. 


2.5.3. MATRICEA UNUI OPERATOR LINIAR 


Să considerăm un operator liniar L:R™%=R® şi două baze: 


B = (Ups bb, B' = {vis Va... v} ale acestor două spatii liniare. Oricare ar 
e2 m 1372 n 


fi vectorii x€ R, y €R" au loc descompunerile: 


Xx= dau, (2.5.5) 
i=] 

y= Bv (2.5.6) 
J= 


"Cum L(u;)e R” putem serie: 
L(u;)= 2 ajv Vis Lm. (2.5.7) 


Numerele reale a; formează o matrice A cu m linii şi n coloane numită 


matricea operatorului corespunzător celor două baze. Ea este evident unic 


determinată şi se schimbă după anumite legi odată cu schimbarea bazelor. 


Rolul acestei matrice este evidenţiat de următoarele consideraţii. În baza 


proprietăţii a) din definiţia operatorului, 


L(x) = L[Seu]=; a L(u,) = a, Banys = 
i=l 
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Comparând rezultatele şi folosind unicitatea coordonatelor deducem: 
m p — | ; 
1= , 
sau, matriceal: 


Yy =L(x)y = A'xXp (2.5.9) 


deci matricea operatorului defineşte modul în care. acesta acționează asupra 
elementelor lui R”. 
Exemplu: Să se determine matricea operatorului liniar L: R? — R? definit prin: 
L(x) SLG 22, 3) (a ra — 332 Fa) în bazele 
B= fu, = (1,1,1)',u; = (0, Diu 3 (0,1,1)' } și 
E eina e 


Soluție: Evident, operatorul este liniar. Pentru determinarea matricei 
calculăm L(u,), L(u,), L(u3). Obținem: 


LD = (UV), LU) 
= L(00, AR 02), LUN =LI) = 0,1) 
Pentru determinarea matricei operatorului trebuie rezolvate sistemele: 


2 
j=l 
Putem aplica metoda transformărilor pentru cele trei sisteme simultan: 
Apiata foma TEE E 007, 9 pe a e d 
“lu 2|3 2 1 (2 1|1 2 0) lo —3|-5 —2 —2 


i Di/A. 2/2 regia 
0 ili 35.213 293 
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iia a UTIN. Ù apt? Pr 
Obținem: TTV A az dou Tia, ATT. ADT Se N 


deci matricea operatorului în cele două baze este: 


| 


> 
Il 

D= win wj= 

wjn II wolu 


şi se poate verifica formula (2.5.9) prin calcul direct. 


Observaţie: Dacă în locul bazelor B şi B’ aa bazele canonice ale celor 
două spaţii, evităm calculul coordonatelor care în acest caz coincid cu 


componentele vectorilor. 


În exemplul considerat, matricea acestuia Ac este dată :. 


Î 2 
A=| l 1 
—1 0 


iar verificarea formulei (2.5.9) este imediată: 


II = 
K . X3 
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2.5.4 VECTORI ŞI VALORI PROPRII PENTRU 
UN OPERATOR LINIAR 


Considerăm operatorul liniar L: R” — R” şi vectorul nenul x din R”. 
Definiţie. Vectorul x se numeşte vector propriu corespunzător valorii 
proprii XE R dacă: rii 


L(x) = Ax | (25.10) 


În baza formulei (2.5.9.) obținem sistemul . 

Ax = Ax mz (2.5.1) 
care se poate serie şi sub forma: i | 

(A'—M,)x=0 ui | fad 2) 

Se remarcă faptul că este vorba despre un sistem liniar şi A care 
trebuie să aibă o soluție nebanală: 


E Es Oaa T 


ce reprezintă vectorul propriu. Condiția de existență a acestei soluții nebanale este: 


rang (A'-N,)<n (2.5.13) 


care este echivalentă cu: det (A! — A) >0. 


Consideraţiile de mai sus permit, practic, determinarea valorilor proprii í şi a 
vectorilor proprii corespunzători acestora.. Astfel, din ecuaţia (2.5.13) numită 
ecuație caracteristică se obțin valorile proprii A, numerele reale sau nu, distincte 
sau nu. Pentru fiecare valoare proprie X astfel determinată din sistemul (2.5.12) se 


obțin vectorii proprii corespunzători. 


Observaţie: Deoarece operatorul liniar este complet determinat prin 
matricea sa A corespunzătoare unei perechi de baze ale celor două spații nu este 


exagerat să spunem că valorile proprii şi vectorii proprii aparţin matricei pătratice A. 
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Exemplu: Să considerăm operatorul L:R? =R? a cărui matrice într-o 
pereche de baze este: 


332 0 


Să se determine valorile proprii şi vectorii proprii cu ai acestui operator. 
Soluţie: Ecuația caracteristică (2.5.13) este: 
3—X —2 0 
-DJEN mE 
0 0 55) 
care conduce la (|—) (5— ÎN = 0 
Valorile proprii sunt >, = =5, Ms. 
Pentru >, =5= A, sistemul (1.1.12) devine: 


„—2X%, — 2% =0 


care are soluția Xx =- =a, x =ß, deci vectorii proprii x corespunzători 
valorii proprii X = 5 sunt de forma: 
x= (a) = a (—1,1,0)'+6(0,0,1) a,BeR. 


Pentru sistemul > =1 sistemul (2.5.12) devine: 


2x44 =0 
= 2 hate. =0 
3x = 


care are soluția: x, = 34, x3 =, x3 =0, deci vectorii proprii care corespund valorii 
proprii X=| au forma: 


x = (0) = (10), BER. 
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Proprietăţile vectorilor proprii şi utilitatea acestora sunt date de următoarele 


teoreme: 


Teorema 2.5.1. Dacă x este un vector propriu corespunzător valorii 
proprii > pentru operatorul liniar L, atunci kx pentru orice keR are 


aceeaşi proprietate. 
Demonstraţie. Prin calcul direct, L(kx) = KL(x) =k: ox = Ag (kx). 


Teorema . 2.5.2. Mulțimea soluțiilor sistemului L(x)=ħọx este un 


spaţiu liniar. 


Demonstraţie. Cele două teoreme anterioare justifică rezultatul din această 


teoremă. 


Observaţie: Mulțimea S), a soluțiilor sistemului din teorema 2.5.2 se 
numeşte subspațiu propriu corespunzător valorii proprii ^o.. Acest subspațiu 


conține vectorii proprii, dar şi vectorul nul care nu este vector propriu în sensul 


definiției. 


Teorema 2.5.3. Doi vectori proprii x,,x» corespunzători la două valori 


proprii distincte X, şi A, sunt liniar independenți. 


Demonstraţie. Prin reducere la absurd presupunem că x, şi x, ar fi liniar 
dependenţi prin urmare: aX; +&2X, = 0 cu măcar un scalar nenul. Aplicăm acestei 
relaţii operatorul L. Obţinem, luând în considerare şi calitatea de vectori proprii: 
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Din sistemul: 

aix ta xX =0, 

AX HAALX = 0 
se obține în cazul a, = 0 


Ay (A — X3)X, = 0 care nu poate fi adevărată. : 


„Observaţie. Proprietatea din teoremă se poate generaliza la mai mult de doi 
vectori proprii corespunzători la valori proprii distincte. În particular, dacă L are n 


valori proprii distincte X}, Nz,.. A, atunci „vectorii xX ES), sunt liniar 


independenți şi formează o bază în R". Are loc astfel, următoarea 


„ Teorema 2.5.4. Într-o bază formată din vectori proprii a lui R” 


matricea operatorului L are forma diagonală. 


Demonstraţie. Din ipoteză rezultă : 
LAE Nx FO xp HOF FO, 
L(x) = 0-x, FAN + X3 = mite pi i Aa 
L(x,) = 0-x, +0-x, + A3ă%3 +e 0- xX, 
şi în baza formulelor (11.7) obţinem matricea A, în baza formată din vectorii 


proprii By = (x, X2.: Xp } 


>, Fotin. 24 
Ap = 0 > . 0 . 
ER 
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2.6. FORME LINIARE. FORME PĂTRATICE | 


In cele ce urmează, vom considera spațiul liniar R” organizat peste corpul 


de scalari R. 


2.6.1. FORME LINIARE 
| Definiţie. O aplicaţie f: Kori R se numeşte formă liniară dacă satisface 
condițiile: 
a). f(x Hx =f x+ fa) (VX XER", 


b) FAV =M A), (xE R” NER: 


= Dacă privim R drept spațiu liniar peste el însuşi, atunci forma liniară f 
devine un operator liniar căruia îi corespunde o matrice:A din Mar, deci o matrice 
cu o singură coloană. | 

A E (PARIU pego 
şi din (2.5.9) obţinem expresiă formei liniare: 


f(x) = At, AER". Si, at = Opl) 


Exemplu: Să se arate că: f : R? — R, definită hes 

EEST h x)= 244x + 6 
este o Pi, liniară şi să se cre matricea sa corespunzătoare bazei A îi a 
lui R? 

Soluție: Ca operator liniar de la R? la R! f este evident o formă liniară 


deoarece singura componentă a lui f(x) este liniară şi omogenă. Pentru 


determinarea matricei formei, calculăm f(e,), f (e2), f (e3): 
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f(e,)= (0,0) aaa 

f(e,)=f{(0,1,0)' 

f (e3) = f { (0,0,1)' 

Matricea A. a operatorului în baza canonică este aşadar: 
Aga D046), 


şi formula (2.6.1) se poate imediat verifica. 


2.2.6. FORME PĂTRATICE 


Considerăm în continuare spaţiul liniar R" definit peste R, precum și 


funcția Q:R" x R” — R definită prin: 
Q(x) = 292 Ajj- Aj aij ii a ji (V) i,j = Ln 5 (2.6.2) 


Această funcție s se nlimeştă pătratică şi generalizează funcțiile de gradul al 


doilea omogene într-o singură variabilă. 
2 


Matricea A E M; (R) care are dierentele ajo 1,]=1,n este matricea formei 


ŞI se depre la exprimarea acesteia. Prin verificare directă se poale arăta că formula: 


Q(x)=x'Ax. (2.6.3) 


Exemplu; Să se scrie matricea corespunzătoare bazei canonice a formei 


pătratice Q: R? —>R definită prin: 
? s] 3 
QJ =2a — 0 3x3 — Xe + 2. 


Soluţie: Folosim (2.6.3) sau (2.6.2) elementul aj al matricei A fiind 


şi a. =a. M A 


coeficientul lui xx; în expresia acesteia. Cum XX, = xx şi aj=aj 


J 


coeficienții produselor de variabile diferite se înjumătățesc. Aşadar: 
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2 —1/2 9 
A =|-—1/2 SIAL 
0 | 3 


şi formula (2.6.3) se poate imediat verifica. 


Observație: Matricea A depinde de baza în care se lucrează. Ea se schimbă 
odată cu schimbarea bazei. Este simetrică conform condiţiei a; = a;;. 
Dacă în matricea A numai elementele diagonalei principale pot fi nenule 
spunem că forma pătratică are expresia canonică. Ea se poate obține modificând baza 


în care se determină matricea A. Expresia canonică a unei forme pătratice are tipul: 


QAYLA HNA e A 2 e A (2.6.4) 


O formă pătratică Q este nedegenerată dacă expresia canonică are exact n 
termeni. Când numărul lor este mai mic forma se numeşte degenerată. 

Expresia canonică are importanță deosebită în studiul semnelor acestora 
având în vedere că Q(x) e R. Există forme pătratice care păstrează acelaşi semn 
pentru orice argument x R” >Nu este greu de verificat că: Q(0) = 0. Următoarele 


definiţii pot fi relevante: 


Definiţie. Forma pătratică Q:IR" — R se numeşte pozitiv (negativ) , 


definită dacă: 


QSO (003) <0) V)x A RE ul SFS 


Definiție. Forma pătratică Q:1R" — R se numeşte pozitiv (negativ) 


semidefinită dacă: 


Q(x) Z0 (Q(x)<0) (V)x ER" (2.6.6) 
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Orice formă pătratică definită este o formă pătratică semidefinită. Reciproca 


nu este adevărată. 


Dacă o formă pătratică este adusă la expresia canonică caracterul de formă 
pătratică pozitiv (negativ) definită (semidefinită) poate fi mult mai simplu apreciat. 
În acest spirit, încercăm să aducem matricea formei pătratice la expresia superior 
triunghiulară folosind metoda transformărilor elementare de linie. Metoda aceasta 
se numeşte metoda lui Gauss. 

aji a sie ain 


] 7 
0 422 ... aag 


RO 26.7. 
O D Da PRS 
5 a Ciana 
Cu transformarea: 
R. p NT e ala 
Yı = an ta fe ani antn 
WA Talart hean (2.6.8) 
Sn 53 anna 
se obține expresia canonică a formei pătratice: 
a] l | i 
O ra 7 mr AP to ape aa (2.6,9) 


A â23 nn 


In acelaşi scop se poate utiliza şi metoda lui Jacobi, dar este necesar ca 


lanțul minorilor principali 


piata A _ mate (2.6.10) 
Arun- AR 


A Sân A, = 


să fie nenul. Expresia canonică este în acest caz: 


| pl wijange Aaly hy 
AT osia “ue DB e yi (2.6.11) 
‘5 l 2 


n 
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„ Exemplu: Să se găsească expresia canonică a formei pătratice: . 
Q(x) = 2x — 2x + 232 — 233 + 333. 
“Folosim mai întâi metoda lui Gauss. Obţinem succesiv: 
EP TE zu aj cz 
A=|=1 2 = 2 103 —2ho 
(e mt a a D i ae | Fc K gt OY a 
1 1 —l 1 1 -=i 
~IO 1 4|=|0 1 4 
0 —1 3 0 07 


Expresia canonică este aşadar: 


l 
C tata 


Prin metoda lui Jacobi obținem: 


MENE T demne 
A =2 a= i Avé A =|-1. 2 -l=3 


Expresia canonică astfel obținută este deci: 


2 d a 
lat f- qg 
Q(z) zá ga 7.3 


Observaţii: 

„1%. Expresia canonică nu este unică. Prin orice metodă s-ar obţine numărul 
de termeni cu coeficienţi pozitivi şi numărul de termeni cu coeficienți negativi este 
acelaşi. 

2%. Forma pătratică din exemplul precedent este pozitiv definită deoarece 
toate pătratele au coeficienţii pozitivi şi valoarea nulă se obține doar când 


variabilele y; (sau z;) sunt nule. 


Capitolul 3 


ELEMENTE DE PROGRAMARE LINIARĂ 


3.1. PROGRAMARE LINIARĂ 


3.1.1. MODELE ECONOMICE CE CONDUC 
LA O PROBLEMĂ DE PROGRAMARE LINIARĂ 


a) Consumul lunar de materii prime 


Consideram un proces de producție în cursul căruia se realizeză m tipuri de 
produse = Im, pentru care se utilizează n tipuri de materii prime M;, j= In. 
Dintr-o unitate de materii prime M; se pot produce a; unităţi de produs P;. Se 
preconizează ca produsele de tipul P, să fie realizate lunar în cantităţile b,. Costul 
unei unități de materie primă M ; este cj. “Se cere a determina consumul lunar de 
materii prime, încât să se realizeze producția planificată, iar cheltuielile legate de 
materiile prime să fie minime. 

Să vedem care este modelul matematic al acestui model economic. Va 


însemna cu Xx; j=l,n cantitățile de materie primă ce urmează a fi utilizate. Va 


trebui să scriem sub formă algebrică restricțiile economice. Aceste restricții sunt de 


trei feluri: restricţii legate de planificarea producției, restricții impuse de sensul 
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economic al variabilelor şi restricția de minimizare a cheltuielilor. Aşa cum a;X) 
reprezintă cantitatea de produse de tipul P, ce se obține dintr-o cantitate x; de 


materie primă M Ja deducem că: 


n 
aX td. TE x, = di a;X) 
par 


reprezintă calitatea totală de produse de tipul P,, ce se obține când materiile prime 
M, se utilizează în cantitățile X; ,j=1l,n. Dar această cantitate nu poate fi mai 


mică decât b; şi deci: 
D ayay biT = m) FOA A) 


Cum x; reprezintă mărimi economice rezultă: 


Es E TES. (3.1.2) 


cheltuielile totale obținute ca sumă a cheltuielilor unitare vor fi date de funcţia: 
n | | T | | 
fa 2C jjs AG Hante a takoki (3.1.3) 
J= 


care se cere a fi minimalizată. 

Din punct de vedere matematic, problema poate fi formulată astfel: să se 
găsească în mulțimea soluţiilor sistemului de inecuaţii (3.1.1) care satisfac 
condiţiile (3.1.2), o soluţie pentru care funcţia f ia valoarea minimă. Problema are 
sens dacă sistemul (3.1. 1) are cel puțin o soluție. Inecuațiile (3. 1.1) reprezintă 
restricțiile problemei, (3.1.2) condiţiile de nenegativitate, iar f poartă denumirea de 


funcţie obiectiv. 
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b) Folosirea optimă a resurselor 


Intr-un proces de producție în care se pot realiza n tipuri de articole 
AI, j=l,n sunt utilizate m tipuri de resurse R;, i=l,m, care sunt limitate de 
cantitățile bj. Pentru producerea unei unităţi de articol de tipul A, se consumă o 


cantitate a; de resursă de tipul R,. Profiturile unitare pentru articolele A, sunt 


Cj, J5l,n. Să se determine ge la x; de articol Aj J= In care urmează a se 


realiza în procesul considerat, cu resurse limitate, asfel ca procesul să aibă drept 


rezultat profitul maxim. 


Limitările resurselor conduc la: 


n — 
2 ajx; Sba i=l,n (3.1.4) 
j=l 


Necunoscutele x;, j=1,n, trebuie să satisfacă condițiile de nenegativitate : 


20, nth (3.1.5) 
Profitul total este dat de funcţia liniară : 
FSE ex (3.1.6) 
j=l 


care trebuie să ia valoarea maximă. 


ŞI acest model constituie un exemplu de problemă de programare liniară. 


c) Problema amestecului optim 


Să presupunem că pentru realizarea unui amestec se pot întrebuința n materii 


prime M;,j=,n, disponibile în cantități nelimitate, ce conțin m substanţe 
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S;i= lm. Dintr-o unitate de materie primă M j» trece în amestec o cantitate lij de 
substanță S;. Substanța S; trebuie să intre în amestec în, cantitatea. b;. Costurile 
unitare ale materiilor prime M, sunt. C; unităţi băneşti. Să se determine cantităţile 
A j=ln de materie primă. M j necesare realizării amestecului cu compoziția: 
prescrisă şi la un cost total minim. 
Restricțiile sunt şi acum de „trei feluri: restricţiile de compoziţie a 
amestecului, pisi ic: datorate sensului economie al aero jE şi restricția de 


minimizare a a costului total. 


Tr da algebrică a primelor două categorii este: 


i Xj i ele (3.1.7) 
j=l i 
xj 20, j=ln | (3.1.8) 


loga. total al materiilor prime necesare pentru realizarea unei i unităţi de 


amestec este dat de funcția: 
n y 
j=l 


care se cere a fi minimizată. 


3.1.2. FORMULAREA GENERALĂ A UNEI PROBLEME 
"DE PROGRAMARE LINIARĂ | 


Se constată că modelele matematice, ataşate unor probleme economice 
prezintă o serie de asemănări, fapt ce permite înglobarea lor într-un model general. 
Să considerăm un sistem care poate conţine GC iapa inecuații de un semn sau altul, 


deci de forma: 
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aa HAEA Fei, Sb, 
a HX Eaa ke 


"<H 


i (3.1.10) 
anm HX + Ama F ++ mn Xn 5 b m 
sau 
n a 
2 ajx; S bn ixl (3.1.11) 


unde semnul S înseamnă semnul =, ori <, ori >, în diferitele condiții (3.1.10). 
Considerăm şi sistemul de condiţii: 
x; 20, j=la i | (3.1.12) 


Introducem şi o funcţie obiectiv: 
f(x)= Zcjă, (3.1.13) 
j=l 


Problema generală a programării liniare se formulează astfel: în mulțimea 
soluţiilor sistemului (3.1.11) ce satisfac condiţiile (3.1.12) să se găsească o soluţie 


ce dă valoarea minimă sau maximă a funcţiei (3.1.13). 


3.1.3. FORME ALE UNEI PROBLEME 
| DE PROGRAMARE LINIARĂ 


a) Forma matriceală 


Considerăm matricele: A de tipul m x n (m < n) cu elementele lj, X 


matrice coloană de componente X; j= n, b matrice coloană de componente 


b., i= l,m, C matrice linie de elemente C; j=l n, atunci putem scrie problema 


de programare liniară sub forma: 
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sau 
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Ax Sb, 
x>0 4 | (3.1.14) 
(min/max) f(x) = C 


b) Forma vectorială 


Dacă introducem vectorii: 
T ; 7 I 
P, al] , P2 = [0,3 nasr. e tai pe Da REL Dei a 


Po = b. 
Atunci problema poate fi scrisă sub forma: 
xF +P, FF Xp Pa S h» 

x >0 (3.1.15) 


(min/ max) f(x) = Cx 


c) Forma canonică 


O problemă de programare liniară este sub formă canonică dacă ea se scrie: 


Ax >b, 
x>0 | (3.1.16) 
(min) f(x) = Cx 


x>0 (3.1.1) 


Ax <b, 
(max) f (x) = Cx 


d) Forma standard 


O problemă de programare liniară este sub formă standard dacă restricţiile 


ei sunt sub formă de epalităţi. Avem deci: 
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Ax =b, 
x>0 sau 
(min/max) f(x)=Cx 


xh +P, e Xa Pa F h, 
x>0 


> (3.1.18) 
(min/ max) f(x)=Cx 


Observaţie: Formele (3.1.16), (3.1.17), (3.1.18) ale unei probleme de 

programare liniară sunt echivalente dacă ținem seama de următoarele: 

l. Sensul unei inegalităţi se schimbă prin înmulțirea cu —1; 

2. O inecuaţie de forma ax, + aX +... + au, < b;, poate fi scrisă ca o 
ecuație aX, Fapt. dim + = b;, Prin adăugarea a unei noi 
variabile (variabilă de compensare) x; >0. 

| a. O inecuaţie de forma: ax, + aX fst a, > bj, poate fi scrisă ca o 
"ecuaţie aX tan Fe tan =—xXi =b, scăzând o variabilă de 
compensare x; >0. 


4. O ecuaţie poate fi scrisă ca două inecuaţii de semn contrar. 


5. Problemele în care se cere (max) f(x) pot fi reduse la o problemă de 
minim pentru funcţia f=-—f. 


„Având în vedere cele de mai sus, ne vom ocupa de probleme aflate sub 


forma standard cu cerința de minim pentru funcţia obiectiv. 


3.1.4. MULȚIMI CONVEXE 

Fie spaţiul vectorial R” şi M o submulțime nevidă asa MC R". 

Definiţia 3.1.1. Mulțimea M C R” se numeşte convexă dacă (V)x Xa EM 
şi (A €]0,1] avem: 


x = x, +(1—2)x, EM (3.1.19) 
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Definiţia 3.1.2. Un vector x EM se numeşte punct extrem sau vârf al 
mulțimii convexe M, dacă nu există x, = x, din M şi nu există ^ €(0,1) astfel 
încât 


x = Àx t (1 A0x2 (3.1.20) 


Definiția 3.1.3. Numim combinație liniară convexă a vectorilor 


x; EM, i=1,k o expresie de forma: 


JF 


. 


Să l | 4 à k i 
AX, cu A EIM] TA El (3.1.2) 
ie] 


Observaţie. Se demonstrează că orice punct ce nu este punct extrem se 


scrie ca o combinaţie liniară convexă de puncte extreme ale mulțimii convexe. 


3.1.5. SOLUȚII ALE UNEI PROBLEME 
DE PROGRAMARE LINIARA 


“Fie o problemă de programare liniară aflată sub forma standard. 
Ax =b, 


x >0 y (3: 1.22) 
(min) f(x) = Cx | i 


Definiţia 3.1.4. Un vector x€ER" ale cărui componente Xp J=l,n 


verifică sistemul de restricții şi condițiile de nenegativitate se numeşte soluție 


admisibilă a unei probleme de programare liniară. 


Mulțimea soluțiilor admisibile o vom nota cu S}. 


Teorema 3.1.1. Mulțimea soluţiilor admisibile S, este o mulţime 


convexă. 
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Demonstraţie: Fie x, = x, ES, două soluţii sataiaitiile- 

Atunci avem: i | 

AX, Sb, x 20; Ax, =b x2>0. 

Să considerăm un vector x de forma:, 

x =x ((—dx2, cu A €[0,1]. 

Să arătăm că x €S, . Într-adevăr vom putea serie: 

Ax = A[Ax, + d-ho] = MAX, +(1— 1)Ax, = Ab+(1— A)b = b 
şi deci verifică sistemul de teich. 


Cum efo, L iar X1, X2 20, rezultă Ax *(1—1)x, =x >0 şi deci x este 


soluţie admisibila a problemei de programare liniară. 


Definiția 3.1.5. O soluție admisibilă care are cel mult m componente strict 
pozitive, iar vectorii formaţi cu coloanele matricei A, corespunzătoare acestor 


componente sunt liniar independenţi, se numeşte soluţie admisibilă de bază. 


Dacă soluţia admisibilă de bază are exact m componente strict pozitive, 
spunem că este nedegenerată, iar dacă are mai puţine, degenerată, Mulțimea 


soluțiilor admisibile de bază nu este o mulţime convexă. 


Exemplul 3.1.1. 


Fie problema de programare: 


2x — N A Za 
Xj T khaa n 


(min) f(x) = 2x) +x + 2x4 
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În baza definiţiilor (3.1.4) şi (3.1.5) constatăm că: 


x4=(243, 1342/3, 513)" este soluție admisibilă 
x2=|0, 1, 0, f este soluție admisibilă de bază nedegenerată 
x3=|0, 0, 2, of este soluție admisibilă de bază degenerată 


T r M PE 3 
X „ll 0,1, 0] este soluție admisibilă, dar nu este de bază 


(vectorii Pı şi P3 nu sunt liniar independenți) 


Teorema 3.1.2. Orice soluție admisibilă de bază a unei probleme de 
programare liniară este un punct extrem al mulțimii soluțiilor admisibile S, şi 


reciproc, orice punct extrem al mulțimii soluțiilor admisibile este o soluție de 


bază a problemei de programare liniară. 


Demonstrație: Să presupunem mai întâi că soluția admisibilă de baxă x nu 
ar fi punct extrem al mulţimii S, . 

Aceasta înseamnă că există x, = x, ES, şi A€ (0,1), astfel ca: 

i e UA (3.1.23) 

Fără a restrânge generalitatea, putem presupune că primele r(r <m) 
componente ale vectorului x sunt strict pozitive şi deci vectorii Pı, P2, . . .„ P, sunt 
liniar independenți. Notând cu x, şi x; componentele vectorilor xı şi xz din 
(3.1.23) găsim: 

xpt A Axa =0, i=r+ln (3.1.24) 


Cum xa şi X; >0 şi à € (0,1) din (3.1.24) rezultă că 


Xil = Ai = 0 „i = r+ ln (3125) 
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Pe de altă parte x, şi x2 verifică sistemul de restricții scris sub formă 


vectorială şi deci avem: 


Fa: Xa +... Fă = B 


3.1.26 
Xa th +. + Xa, R ( 


Cum vectorii xı şi x2 sunt liniari independenți, relaţiile (3.1.26) au loc 
numai dacă i 
ora ia u i 9 (027) 
Din (3.1.25) şi (3.1.27) rezultă că x, = x2 ceea ce AREN presupunerea 
făcută că x nu ar fi punct extrem al mulțimii S. Pentru a dovedi partea a doua a 
teoremei va trebui să arătăm x € S, » punct extrem al ei este soluția admisibilă de bază. 
Să prosper că primae s componente ale vectorului x sunt strict pozitive 
ŞI că vectorii Pı, P>, PR nu sunt liniar independenti. Această presupunere din 
urmă arată că există ],, i= 1,5 nu toți nuli, astfel încât: 
-API FAP, +... +A 20 (3.1.28) 
Dar x este soluție admisibilă a problemei şi deci: 
XP Ha Pa «e see ha R ==) PY (31.29) 
„ Adunând şi scăzând din relaţia (3.1.29) relaţia (3.1.28) înmulțită cu az0, 


obținem: 


(xa + QÂ) B+ (x, +a), )P, ZA + (a, FaA JP = Pi (3.1.30) 
AQ 64 E (i — adu E F + (2 — ah EN B: 


Putem alege a suficient de mic încât: 


e leg = es (3.1.31) 


n a a a 


Do o a a 


72 MATEMATICI APLICATE ÎN ECONOMIE 


In baza relaţiilor (3.1.31) vectorii x, = (x4), X3 = (3) cu xi to), 
Xi = A +QĂ, pentru i =1,5 şi X =x =0 pentru i=s+l,n definesc două 
soluţii admisibile ale problemei. Mai observăm că: 


pia îi | 
X= x+ xX, 34.32 
Aa A ig ( ) 


relație din care rezultă că x nu este punct extrem, aşa cum noi am presupus. Aşadar 
vectorii P}, P>,..., P; sunt liniar independenţi şi cum s nu poate fi mai mare ca m, 


rezultă că x este soluția admisibilă de bază a problemei de programare liniară. 


Definiţia 3.1.6. 0 soluție admisibilă x se numeşte soluție optimă, dacă 


realizează minimul sau maximul funcției obiectiv f. 
Se poate demonstra că mulțimea soluţiilor optime este o mulțime convexă. 


Teorema 3.1.3. Dacă o problemă de programare are o soluție optimă 
finită, atunci există un punct extrem al mulțimii soluţiilor admisibile S, în 


care funcţia obiectiv ia aceeaşi valoare. 


„Demonstraţie: Fie xo o soluţie optimă finită a problemei, adică: 


f(x 9)=(min) f(x) 3 (3.1.33) 


XES, 
Dacă xo este punct extrem al mulțimii S, , atunci teorema este demonstrată. 


Să presupunem acum că xo nu este punct extrem. Atunci în S, există puncte 


extreme X;, i=l,r încât 


PALA e cu A €[0,l] - şi i | e fe d (3.134) 
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În baza relaţiei (3.1.34) şi a liniarității funcției f putem scrie: 
f(xo)=f po = Ent (341.35) 
i=] i=] 


Să notăm cu Xs punctul extrem pentru care f(x) = (min) f(x), 1=},r . Dacă 


în relaţia (3.1.35) în loc de f(x;), i=1,r punem f(x.) Şi având în vedere că 


X > 0, FN =] găsim: 
i=l 


fo > ra ata) a (3.1.36) 
i=l 


Din (3.1.33) şi (3.1.36) rezultă f(x) = f(x.) şi teorema este demonstrată. 
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3.2. METODA SIMPLEX ȘI ALGORITMUL 
SIMPLEX PRIMAL 


3.2.1. PRINCIPIILE METODEI SIMPLEX 


În baza teoremelor 3.1.2 şi 3.1.3, vom căuta soluțiile optime ale unei 


probleme de programare liniară în mulţimea soluţiilor admisibile de bază care este 
o mulțime finită având cel mult Cp elemente. Metoda simplex elaborată de G.B. 


Dantzig are avantajul că pe lângă un volum mic de calcule, permite a da răspunsuri 


la toate situațiile ce se pot ivi. 


Se porneşte de la o soluție admisibilă de bază oarecare, ce este supusă apoi 
unui criteriu de optimalitate. Dacă acesta este satisfăcut, atunci soluţia considerată 
este soluția optimă a problemei. Dacă nu este satisfăcut, se indică procedeul prin 
care se ia în considerare o altă soluţie admisibilă de bază nu oarecare, ci una pentru 
care funcţia obiectiv să ia o valoare mai mică decât pentru precedenta. Se indică şi 
situațiile în care problema admite optim infinit, din punct de vedere economic 
însemnând că nu are soluţii. După un număr finit de paşi se ajunge la soluţia 
optimă sau la situația că admite optim infinit. 

Să considerăm problema de programare sub forma standard: 


Ax € D, 


x>0 (3.2.1) 
(min) f(x) = Cx 


sau 
Ar Pb, Fe HR =Pj 
x > 0 (3.2.2) 
(min) f(x) = Cx 


Să presupunem fără a restrânge generalitatea că matricea lărgită a 


sistemului de restricții este de forma: 
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l 0 ... 0 CO ... Qin 8, 
0 J-LO Com se op bə (3 2 3) 
OO mem La ee enl is 


Vectorii Pi, P2, . . ., Pm formaţi cu primele m coloane ale matricei A, fiind 
liniar independenţi, formează o bază în R™. În acest caz, un vector oarecare 


P;, j=m+1,n se exprim în mod unic cu ajutorul vectorilor bazei sub forma: 


m 


P=0j Pi kaz Pat... +0. T E E , (3.2.4) 


Să notăm cu Cp vectorul TA format cu primele m componente ale 
vectorului C, adică coeficienții variabilelor X} x,,...,x„ din funcția obiectiv ce 
corespund vectorilor care au intrat în bază. 


Introducem mărimile: 
m PE, 
zj = CP, = co aj; j=1,n ORS) 


Dacă presupunem acum că ß;, i=1,m sunt pozitive şi dacă avem în vedere 
(3.2.3) deducem că vectorul X de componente, 
X =P 5 Base Xn Bms Im = 0e E0, 


este o soluție admisibilă de bază a problemei de programare considerată. 


Teorema 3.2.1. (Criteriul de optimalitate) Dacă diferențele 
Z;—Cj <0; j=m+l,n, atunci problema de programare liniară are optim finit 


şi X este soluția optimă căutată. 


Demonstraţie: Fie x o soluţie admisibilă oarecare de componente 


X; j=l,n şi X soluţia admisibilă de bază construită. 
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Dacă avem în vedere (3.2.3) vom putea scrie restricţiile sub forma: 


x, + 5 QX) =, i=l,m (3.2.6) 
j=m+l 


Să evaluăm diferența dintre valorile funcţiei obiectiv pentru cele două 


soluții considerate. Având în vedere (3.2.5) şi (3.2.6) găsim: 


“m n Pm | 
fX) — fa) > if d ety =2 C; Xt 22 uxi- DCN 


j=m+l Na, j=l gut 
m m m 
d c;ă + Sc; S ax; |- —Dejă, =) c 5Y QX; 4/3 Cj: (3.2.7) 
i= i=l j=m+l j=l 18] j=m+l | j>m+l 
n m n 
== di (Sea, Fa Că, = s (z; =C) 
Jz eml i=l j=m+l j=m+l 


In baza ipotezei z; —c; 20,j=m+l,n şi deoarece x; 20, j=l,n din 


(3.2.7) rezultă că: 


f(x), (V)xE S, onl (3.2.87 


adică X este soluția optimă a problemei: 


Observaţii: . 


Z; —c€; =0 pentru j= l,m 


a) se constată uşor că i 


b) In problemele de programare, în care avem cerința de maxim pentru 
funcţia obiectiv, criteriul de optimalitate va fi: 

z; Cc; 20, j=m+ln 
c) Din (3.2.7) rezultă că în cazul în care există z; —c; >0 în probleme de 


minim pentru funcţia obiectiv pot exista soluţii admisibile pentru care funcția 


obiectiv să ia o valoare mai mică decât pentru X. 


In cazul în care X nu este soluţie optimă, să vedem cum putem considera o 


altă soluție de bază. Aceasta se poate obține: dacă un vector P; din bază se 
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înlocuieşte cu un vector P) care nu face parte din bază. Să: notăm cu yk 


k=1,2,...,1—1, ji +1,...,mcomponenetele noii soluţii de bază Y . 


„Din (1.2.5) înlocuind b, cu, Bp, Api cu Œj, X CU Yy, x, cu B; obținem: 


=- Bi q 2w ay #0, ksl, 2ni h i Him 
d (3.2.9) 


Teorema 3.2.2. (Criteriul de ieşire). Condițiile necesare pentru înlocuirea 


vectorului P; cu Hlor ul P; sunt: 


a; >0, Bebe (a, >0) = (3.2.10) 
O 


ij kj 


Demonstrație. Punând condiția ca Y să fie soluție admisibilă de bază. Din 


(3.2.10) rezultă că trebuie să avem: 


A kae LEN A 
Ea t Aij 


A. 50, ECER (a, >0) 
O 
şi în felul acesta teorema este demonstrată. 


Să vedem acum ce condiţii trebuie să îndeplinească un vector P; pentru a 


„intra în bază î în locul „Aso lua P; pentru a obţine o SANItje „mai bună“. 


Teorema 3.2.3. Condiţia p pe. care trebuie să o zor men un vector P;, 


j=m+l1,n pentru a intra în bază în Jocul vectorului P; în vederea obținerii 


unei my pia „mai bune“ aste: 
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Demonstraţie: Să evaluăm diferenţa f(Y)= f(x). Vom avea: 


m m 
IOD D= z CkYk +CjYj t È Ck = Ș> REN — © CkÊk Gibi (3.2.12) 
i+ k=i+1 
Dacă avem în vedere (3.2.9), relația (3.2.12) devine: 
f(Y)-f(X)= a - Bay, +, Ba E pa -B |= Abg 
ij ! tjj k=i+l Qij O, 
sau 
m d f 
Qij k=l Qij Aj 
Cum a; >0, B; 20, pentru a avea f(Y)-f(X)<0, | va trebui ca 


Z; —Cj 20 şi astfel teorema este demonstrată. 


3.2.2. ETAPELE ALGORITMULUI SIMPLEX PRIMAL 


Algoritmul simplex primal constă dintr-o succesiune de iterații, care 
îmbunătățesc treptat o soluţie iniţială admisibilă de bază. Pe baza considerațiilor 
din aliniatul 3.2.1 putem pune în evidenţă etapele algoritmului care ne va conduce 


la soluția optimă a problemei dacă aceasta există. 


l. Se determină o soluție admisibilă de bază, fie aceasta X. 


2. Se calculează diferențele Z, —Cj. j=1,n 
a), dacă z; —c; S0, j=m+l1,n, atunci după criteriul de optimabilitate, 
soluția admisibilă de bază X este optimă. 
b) Dacă există indici j, j=m+l,n pentru care z, —C; 20, iar toate 


componentele vectorului P; sunt mai mici sau egale cu zero 


(œ; <0), atunci problema de programare are optim infinit. 
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c) Dacă există indici j j=m+l1,n pentru care 2: me >0 şi nu toți 


G,; S 0, atunci se trece la etapa următoare. 


3. Se aplică criteriul de intrare, determinând vectorul Pj care intră în bază 


ca fiind acel vector pentru care: 


Z; —C. = max (Z,—c,) 
J J Zi —C4 20 ; E 


Se aplică criteriul de ieşire. Părăseşte baza vectorul P; pentru care : 


Pi = min Be 


h ijol iLOR ag O, 


4. Se face schimbarea de bază obținând o nouă soluție admisibilă de bază 
© Y, pentru care etapele algoritmului se reiau. 


Datele problemei de programare liniară şi elementele ce au intervenit în 


teoremele 3.2.1, 3.2.2, 3.2.3 le sintetizăm în tabelul simplex: 
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„Observaţii: 


a) Valoarea funcţiei obiectiv pentru soluția admisibilă de bază este: 
m 
f(x) = cip = CpPo 
= 


b) Diferențele Z; —C;. vor fi: 


Pentru j=1,m avem: Z, —C=0. 


Dacă soluția nu este optimă (există ze, 2/0 ) alegem diferența pozitivă 
care este cea mai mare determinând în felul acesta vectorul P) care intră în bază. 
Dacă există mai multe astfel de diferențe, atunci poate intra în bază oricare dintre 
vectorii ce corespund acestora. Dacă pe coloana vectorului P; care trebuie să intre 
în bază, toate elementele sunt mai mici sau egale cu zero, atunci problema admite: 
optim infinit. Dacă există şi elemente strict mai mari ca zero, considerăm toate 


rapoartele dintre componentele lui Po şi componentele strict pozitive ale vectorului 


Pj, rapoarte de forma — cu a,,>0. Cel mai mic raport va indica vectorul care 
a; 


ij 
părăseşte baza. Dacă vectorul P; părăseşte baza, atunci noua soluție admisibilă de 
bază se obține cu ajutorul unui pilotaj cu element pivot tij- 
c) Problema de programare liniară admite soluție optimă unică în situațiile: 
Z; A c;<0 , j= m+l,n, sau (3)z, — cj = 0 dar toate elementele de pe coloana 
vectorului care poate intra în bază, sunt mai mici sau egale cu zero (&, <0); 
d) Problema de programare liniară admite mai multe soluţii optime dacă există 


Z;—C;=0, j=m+l,n, iar pe coloana vectorului care poate intra în bază 


sunt elemente strict pozitive. 
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Exemplul 3.2.1. Pentru fabricarea a 3 tipuri de produse P, j=13 se 


utilizează 2 tipuri de resurse R1 şi R2. Consumurile specifice, cantitățile de resurse 


şi profiturile unitare sunt date în tabelul: 


Limitări 


Tesurse 


Să se determine cantităţile de produse P;, j=1,3 ce urmează a fi realizate 
pentru a obține un profit maxim. 
Soluție: Notăm cu x,,x>,x3 cantitățile de produse de tipul P}, P}, P3 ce 
urmează a fi realizate. Modelul matematic al problemei este: 
Xa 2 < 304 
2m +33 FF 3 S57 
X; > 0, j=13 
(max) f(x) = 4 ka + 5 
Aducem problema la forma standard folosind variabile de compensare x, şi 
x$ şi considerăm problema de minim pentru funcția -f . 
Vom avea problema: 
Sr Xa 2x t x = 30, 
aa Zi tă tak 35, = 57 


x; 20, j = 1,3, x; >0, x5 >0. 
(min) (=f(x))= —4x x 55% 


“O soluţie admisibilă de bază este x = [0, 0, 0, 30,57]! ou; 
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Tabelul simplex se prezintă astfel: 


În prima etapă au intrat în bază vectorii x$ şi x$. Pe linia diferențelor 
Z, —c, observăm că nu toate sunt mai mici sau egale cu zero. Cea mai mare 
pozitivă este Zj-Cj = 5, şi după criteriul de intrare vectorul P3 va intra în bază. Am 


considerat rapoartele ZE = 2 =15, ii = =57. Cel mai mic raport 'este 


. Olqz 53 


—*+ , aceasta înseamnă că vectorul P; părăseşte baza fiind înlocuit de vectorul P3 . 
a | 
43 


Un pivotaj cu element pivotor 2 încercuit duce la o nouă soluție admisibilă de bază 


x= |0, 0, 15,0, 42) care nu este optimă. Un pivotaj cu elementul pivot 3/2 


E m ma, T 
încercuit duce la o nouă soluție admisibilă de bază x=[28, 0, 1, 0, 0|- care este 
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optimă. Lăsând la o parte variabilele de compensare, obținem soluția optimă a 
| Re = zeta . N | 
problemei x =|28, 0, 1] adică urmează a se realiza 28 unităţi de produs Pī, o 


unitate de produs de tipul P3 obținând un profit de 117 unităţi băneşti. 


Exemplul 3.2.2. Să se determine soluţiile optime ale problemei de 
programare liniară: 
a Xyk = ri 
Xx 4x tx Hax El, 
Xj 2 0, j pe Ik 5 f 
(min) f(x)) = —3x, — 23 4x | 


o soluție admisibilă de bază este x = [0, 0, 0,7, |: Tabelul simplex se prezintă astfel: 


Algoritmul se opreşte găsind soluţia admisibilă de bază x =|0, 0, 1, 8, 0). 


deoarece vectorul P2 care ar trebui să intre în bază are toate componentele negative. 


Problema de programare liniară considerată admite optim infinit. 
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Exemplul 3.2.3. Să se determine soluţiile optime ale problemei de 
programare liniară: 
| X T 2x3 + X4 pr 2X5 = er 
Yt G tx +2 =2, 
x;20,j=1,5 
“Umin) f(x) ==2x -3x + 4x — Sax 2x; 


O soluție admisibilă de bază este: x = [3, 2 010, 0) , 


Tabelul simplex are forma: 


Etapa I 
Etapa II 


Etapa III 


p În E ge z — : T 

In prima etapă am găsit soluția optimă a problemei x = [3, ZMO; O. 0] cu 
(min)f(x)=— 12. Observăm pe linia diferențelor că za—Cc+=0 şi Zs-c5=0. După 
criteriul de intrare vectorul P, sau P, poate intra în bază. Am continuat algoritmul 


introducând în bază vectorul P, în locul vectorului P şi cu un pivotaj cu element 


A a N ză. a, T 
pivot 1 încadrat în etapa a doua am găsit soluția optimă x =, 0, 0, 2, 0] cu 
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(min)f(x)=—12. Din nou observăm că în dreptul vectorilor P> şi P; care nu fac 
parte din bază diferenţele Z—02=0 şi zs-cs=0. Am continuat algoritmul luând pivot 


elementul 2 încadrat, introducând vectorul P; în locul vectorului P,. În etapa a 


C3, l nc PTS CX à 

treia am găsit soluția optimă x = (5, 0, 0, 0, 1] ŞI (min)f(x)=—12. Algoritmul s-a 
oprit deoarece vectorii care ar mai putea intra în bază au mai fost, Ştim că mulțimea 
soluțiilor optime este o mulţime convexă. Dacă însemnăm cu X1, X2, X3 soluţiile 
optime găsite, atunci constatăm că problema admite o infinitate de soluţii optime de 
forma: 


Soluţia generală a problemei considerate va fi: 


| dă 
x=[3A +A, +5, 2040, 242,43], A €[0,1], x =l. 


Exemplul 3.2.4. Să se determine soluțiile optime ale problemei de 
programare liniară. 
2x — 3x32 + xy >4 
A, > 0, j=1,3 
(min) f(x) = 2x xy + 2% 
Soluţie: Aducem problema la forma standard, folosind variabilele de 
compensare. Avem: 
X (dog 23 A = 3, 
DZ IX x X = 4, 
„x; 20, x, >0, xş >0, 
(min) f(x) = 2x + x + 2x, 
Sistemul nu este sub formă explicită şi deci nu recunoaştem o soluţie 
admisibilă de bază. Vom explicita deci sistemul având grijă ca termenii liberi să fie 


pozitivi. 
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Pivotajele în matricea lărgită conduc la: .. 


(DD —2 roz 
2 #83 iroda 


Problema de programare devine: 


3 


l 205 ai Hb 
4 | 


2223 it 00 alja 
Oi turbo BGÂ2ji hhi KAn? 


HA Pe 3; 
xy ai m 2x + xi =R, 
MEE TEETE 
(min) f(x) = 2x tx + 233 


O soluţie admisibilă de bază este x =f; 0, 0, 0, 2). Tabelul simplex se 


prezintă sub forma: 


Problema admite o infinitate de soluţii de forma: 

xx, A(x, 4 €[0,1], unde x, =[3, 0, 0, 0, 2] şi x, =[1, 0, 2, 0, of" 
şi deci x =[1+2), 0, 2—21, 0, 21], à € [0,1]. 

Cum x$ şi x5 sunt variabile de compensare, soluția problemei inițiale este: 

x = 1 20 x= 0, 33 = 2—2, 1€[0,1], (minf&)=6. 


Nu am mai continuat algoritmul deoarece PF care ar putea intra în bază a 


mai fost în prima etapă. 
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3.2.3. METODA CELOR DOUĂ FAZE 


Dacă sistemul restricţiilor nu este sub formă explicită, nu putem recunoaşte 
direct o soluție admisibilă de bază. Am putea utiliza transformările elementare 
trebuind în acelaşi timp să avem grijă ca termenii liberi să rămână pozitivi. Vom 
indica o metodă mai lesne de aplicat în practică pentru a obține o soluţie admisibilă 


de bază şi anume metoda celor două faze. 


Să considerăm că problema de programare liniară (3.2.1) este sub formă 
standard, termenii liberi sunt pozitivi şi că matricea A a sistemului de restricții nu 
conține nici o coloană a matricei unitate. Vom adăuga la fiecare restricție o 


variabilă nenegativă, numită variabilă artificială. Fiind m restricţii, vom adăuga m 
variabile. artificiale x, k=n+l,n+m, şi considerăm problema de programare 


liniară: 


Ax + |x= b, 
X20 Bi (3.2.13) 


i a n+m d 
(min) g(x )= 32 xk 


k=n+1 


Aici Im este matricea unitate de ordinul m și x° vectorul variabilelor 


artificiale. Se constată că orice soluţie admisibilă a problemei (3.2.14) în care 
variabilele artificiale sunt nule devine după înlăturarea acestora o soluție admisibilă 


de bază a problemei (3.2.1). 
Metoda celor două faze constă în: 


Faza întâi. Se caută soluţia optimă a problemei (3.2.14). Pentru ea 
dispunem de o bază şi cum termenii liberi ai sistemului de restricție sunt 


nenegativi, găsim pentru ea o soluție admisibilă de bază. 
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Aplicând problemei (3.2.19) algoritmul simplex primal ajungem la 


situaţiile: 


a) (min) g(x")=0 şi nici un vector corespunzător variabilelor artificiale nu este 


în bază. 


b) (min) g(x)=0 şi în bază se găsesc şi vectori corespunzători variabilelor 


artificiale. 
c) (min) g(x°)=0. 


Faza a doua. În situaţia a) soluţia optimă a problemei (3.2.14), după 
înlăturarea ‘variabilelor artificiale, reprezintă © soluție admisibilă nedegenerată 
pentru problema (3.2.1). Se iau în considerare elemente legate de problema inițială 
(vectorul Cp) şi se aplică algoritmul simplex primal. în „vederea găsirii soluţiei 
optime. Din tabel vor fi eliminate elementele ce ţin de variabilele artificiale. 

Situaţia b) conduce la o soluție  admisibilă de bază degenerată pentru 
problema iniţială însă nu vor fi eliminate datele ce sunt legate de variabilele 
artificiale, pentru care vectorii corespunzători au rămas în bază. În situația c) 


tragem concluzia că problema iniţială nu admite soluții admisibile. 


Exemplul 3.2.5, Să se determine soluţiile optime ale problemei de 


programare liniară: 


x t3 2aăy > 3, 
Să 3% Pie S l 

x; > 0, j=1,3 
(min) f(x) = 3x, — 4x, + 2% 


Aducem problema la forma standard, ea devenind: . 


x, +3x + 2x33 =3, 
30, + Bă _X3—xş =4, 

x; 20, j=1,3, x; >0, x5 >0 
(min) f(x) = 3x — 4x F2 
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Deoarece matricea sistemului de restricţii nu conține nici o: coloană a 


matricei unitate, vom întrebuința metoda celor două faze. 
In prima fază căutăm soluția optimă a problemei: 
x H3X F 2x — x +3 = 3, 
3%. t3x Xa a H =4, 
? AN 71, dă 
Xi 20, Xy X spy >0 
. a a A 
(min) g(xi ) = xp +3 


Soluţia admisibilă de bază de pornire este x = [0, 0, 0,0, 0,3, aj". 


Tabelul simplex are forma: 


2/33 -1/3 


—l/2 1/6 


1/2 —1/2 —1/2 


Am găsit soluția optimă X (1; 5/6, 0,0, 0,0, 0). Soluţia admisibilă de 


bază pentru problema iniţială va fi: x = [1/2, 5/6, 0, 0, 0|" , iar sistemul de restricții 
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cel ce rezultă din ultima etapă a algoritmului simplex pentru problema artificială. 


Vom avea; 


X =a LEA -+ Xs + 245 = 4, 
De etate AB 3, 
X, HX Xa —X3 = 2, 
x; 20, j=L6 
(min) f(x) = —2ax, — 3 — x3 + 4x4 +5% + Bag 


Problema admite optim infinit deoarece vectorul PS nu are şi coordonate 


strict pozitive. 


Exemplul 3.2,6. Să se determine soluţiile optime ale problemei. de 
programare liniară: | 
X1 — 2NF K, k2xg =4 
3 Xp x5 = 8x63 
XAS Xa -Xy Xs Xg 32 
xta =d, 6 
(min)f(x) = —2x, — 3x, — x5} + 4x {+ 5x; + 8X6 
Prima şi a treia coloană a matricei A sunt coloane ale matricei unitate de 
ordinul trei şi atunci vom adăuga la restrictia a doua o variabilă artificială $; . 


Vom căuta în prima fază soluția optimă a problemei. 
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3xa Xy kg 8X HXi =3, 
Ki Xa Xq —Xg3Xg = 2. 
Poe „A 
x, >0,j=1,6, x7>0 
. : q a 
(min)g(x')= x} 


„O soluție admisibilă de bază pentru această problemă este: x =,[4, 0, 2,0,0,0, 3], 


“Tabelul simplex are forma: 


—1/3 1/3 =- -8/3 1/3. 
4/3 -1⁄3 -1⁄3 


Soluția admisibilă de bază pentru problema inițială va fi x =[4,1, 1,0,0,0,].. 


Tabelul simplex pentru problema iniţială are forma: 


Soluţia optimă a problemei este cc, (min)f(x) = —12. 
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3.2.4. CONVERGENȚA ALGORITMULUI SIMPLEX PRIMAL 


Pentru o problemă de programare liniară prin aplicarea algoritmului simplex 
primal, constatăm că valoarea funcţiei obiectiv scade (în probleme de minim) de la o 
etapă la alta şi după un număr finit de paşi găsim soluţia optimă a ei dacă aceasta 
există. Spunem în această situație că algoritmul, simplex. primal este convergent. 
Dacă valoarea funcţiei obiectiv rămâne aceeaşi după câteva iterații succesive, este 
posibil să găsim o soluţie admisibilă de bază care a mai fost considerată. În această 
situaţie, procesul poate continua fără a ne conduce la soluția optimă cu toate că 
aceasta există. Spunem că are loc un fenomen de ciclaj, acesta apărând când intervin 
soluții admisibile de bază degenerate. Soluţiile degenerate apar atunci când în 


aplicarea criteriului de ieşire nu găsim un raport minim unic. De exemplu dacă: 


Bi Pk B, 


min—l- = = ZE 


Qjo Oyj Osj 


atunci eliminând din bază vectorul P;, în urma soluției. componentele Mog; <- žy VOL 


fi epike cu zero. Pentru evitarea fenomenului de ciclaj, se procedează astfel: 
locul vectorului Po, în criteriul de ieşire, se iau vectorii Pj; începând cu vectorii 


bazei până se obține un raport minim unic. 


Exemplul 3.2.7. Să se determine soluţiile optime ale problemei de 


programare liniară: 


Xy 2% xX + 2x56 = 4, 
X2 FSK, t x =6, 
>0,j=1,6 


(min)f (x) = —2x, +x, + 3X3 +5xX, +x; —10xg 
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O soluţie admisibilă de bază este x = [4, 6, 8, 0,0, 0f", 


Tabelul simplex are forma: 


16.0 

În prima etapă întră în bază vectorul P4. Aplicând criteriul de ieşire, găsim 
toate rapoartele egale 4/2, 6/8, 8/4. Luând în locul lui Po vectorul Pi, găsim două 
rapoarte minime 6/3, 3/4, aceasta însemnând că poate părăsi baza vectorul P> sau 
P3. Cu ajutorul vectorului P2 găsim un singur raport minim, 0/4 şi deci P3 va părăsi 
baza. În etapa a doua găsim două rapoarte egale şi utilizând vectorul P, la formarea 


rapoartelor, constatăm că vectorul P părăseşte baza în locul lui intrând vectorul Pg. 


Soluția optimă a problemei este x =|0, 0, 0, 2,0, 0) min f = 10. 
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3.3. PROBLEMA TRANSPORTURILOR j 


Să considerăm că un anumit produs este stocat în m depozite D; în cantitățile 


aj, i=l,m. Produsul este solicitat de n centre de consum C; în cantităţile 
bj, j=l,n, Costul transportului unei unități de produs de la depozitul D; la centrul 
Cj este a; unități băneşti. Se cere a se determina cantitățile x, de produs care 


urmează a fi transportate de la depozite la centrele de consum, astfel ca disponibilul 
să fie epuizat, în fiecare depozit, cererea să fie satisfăcută exact, în fiecare centru de 


consum, iar costul total al transportului produsului să fie minime. Se presupune că 
m îi | e- -y . . . . . v 

2a; =È b; , adică disponibilul din depozite este egal cu cererea totală a centrelor de 
iz] gÈ] 


consum. O astfel de problemă este denumită problema de transport echilibrată. 


Restricţiile problemei sunt din nou de trei tipuri: restricții asupra disponibilului 
ŞI cererii, restricțiile datorate sensului economic al variabilelor și restricţia minimizare 


a costului. 


Datele modelului economic le prezentăm într-un tabel de forma: 
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Restricţiile disponibilului a cererii exprimă faptul că transportul epuizează 
disponibilul fiecărui depozit D; şi satisface exact cererea fiecărui centru de consum 
Cj. Restricţiile datorate sensului economic al variabilelor sunt analoage celor din 


modelele precedente. Traducerea algebrică a celor trei tipuri de restricții conduce la: 


ec ot (3.3.1) 
xy 20, i=l, ITI n 


T 
(min)f(x) = = 2 = Ci ij» PR Xin] 
a) 9) 


„Observaţie. Dacă disponibilul este mai mic decât cererea sau mai mare, 
pentru a echilibra problema se introduce un depozit fictiv sau un centru fictiv cu un 
disponibil egal cu diferenţa dintre disponibil şi cerere și o cerere egală cu diferenţa 


dintre cerere şi disponibil, costurile fiind considerate nule. 


| Teorema 3.3.1. Problema echilibrată de transport echilibrată (3.3.1.) 
admite cel puțin o soluție admisibilă, iar o soluţie de bază admisibilă de bază 


are cel mult m+n-l componente diferite de zero. 


m pab: a 
Demonstraţie: Cum Xa, =5b; = S, considerând x, ini avem: 
i=l j=l 


și Cum, X, „20, ajungem la concluzia că X(x;) verifică aS de restricții şi 


condiţiile de nenegativitate, fiind deci solul admisibilă a tileacigi de transport. 
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Putem observa că din cele m + n restricţii doar.m+n = sunt independente. Este 
suficient pentru aceasta, de exemplu, ca din suma primelor m restricții, să: scădem 
suma următoarelor n-—l şi vom obţine atunci ultima restricție. Pe de altă parte ştim 
că numărul componentelor nenule ale unei soluţii de bază admisibile de bază este 
cel mult egal cu numărul restricțiilor independente şi în felul acesta teorema este 


demonstrată. 


3.3.1. SOLUȚII ADMISIBILE DE BAZĂ PENTRU 
O PROBLEMA DE TRANSPORT - 


Soluţiile optime pentru o problemă de transport vor fi căutate în mulțimea 
soluțiilor admisibile de bază şi di~ acest motiv vom pune în evidență două metode 


prin care poue gäsi aceste soluții. 

a) Metoda diagonalei. Se i începe cu determinarea componentei corespun- 
zătoare căsuţei din stânga sus a tabelului, GECI Xj- 

Vom lua x = min(a, b) şi vom atribui valoarea zero tuturor TEsableloi 


de pe aceeaşi linie sau coloană cu AE , în funcție de următbărele situaţii: 


1. dacă a, <b,, vom lua x, =a, x; =0, j=2,n 

2. dacă a, >B, „Vom lua x = bj x, =0, i=2,m 

3. dacă a, = b,, vom lua x, =4 =b, iar xp sau x, =0, 
Și restul componentelor de pe fii şi coloana lui xy. egale cu zero. În ultima 
situație obținem soluții degenerate. 

Modificăm apoi a şi bı înlocuindu-i cu af = a~ Xi b = bu. 
Procedeul se repetă la pasul următor pentru. calculele rămase necompletate care 


formează un tablou cu m- 1 linii şi n coloane ; m linii şi n—1 coloane ; sau m-—l 


linii şi n —1 coloane după cum primul pas s-a desfăşurat î în situaţiile 1), 2), respectiv 3). 
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< Prima componentă ce va fi determinată acum va ocupa aceeaşi poziție (colţul 
stânga — sus), în noul tabel. Procesul se termină în cel mult m+n -—1 paşi, la fiecare pas 
determinându-se complet o linie (situația 1), o coloană (situația 2), sau o linie şi o coloană 
(situaţia 3). La fiecare pas se determină o singură componentă nenulă. Componentele 
nebazice nu se completează în tabel, căsuţele corespunzătoare rămânând libere pentru a nu 


se confunda componentele nebazice cu eventualele componente bazice nule. 


b) Metoda costurilor minime. În această metodă se va determina mai întâi 


componentă xij pentru care c, 


j = miney, luând x; = min(a;,b,) “ca în metoda 


diagonalei... . 


“Exemplul 3.3.1. Să se determine o soluţie admisibilă de bază pentru 


problema de transport: i 


a) Metoda diagonalei. Problema este echilibrată. Vom căuta o soluție cu 
m+n—l=3+4-—1=6 
35,0 
40,25, 15,0 


20, 10,0 ` 
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Pasul întâi: deoarece min(a,,b,)>a=35, punem conform situației 1) x,, = 35, 
X2 > 13 =X = 0. Modificăm pe a, şi b, înlocuindu-i cu a = a, S S0: 
bh =b =x 15. 

Pasul al doilea: se pleacă de la tabelul format prin înlăturarea Aha linii: 
deoarece min(a, bb! =15, punem (conform 2)) x21 = 15, x4 =0 (33, componentă 
nebazică). Modificăm pe a; şi bý înlocuind cu.a, = a; —x3j =25, bl'= bl = xy =0. 

Pasul al treilea: a rămas de completat tabelul format din liniile a doua, a 
treia şi ultimele trei coloane. Cum min(a,,b,)=b,=10 punem X22 > 10, Xp =0 
E EPONA nebazică) înlocuindu-i cu ap =a- Xy = 155 bj = b= x 0, 

Pasul al patrulea: a rămas de completat tabelul format din ultimele două. 
linii şi coloane. Cum min(a5,b;)=a”=15, punem (conform 1)) x =15, x, =0 
(X24 componentă nebazică). Modificăm pe a şi b, înlocuindu-i cu 
az= 0) — X 20, by by — pp NO 

Pasul al cincilea: a rămas de completat tabelul format din ultima linie şi 
ultimele două coloane. Cum min(a3,b1)=10 punem x33 = 10. Modificăm pe a, şi 
b; înlocuindu-i cu dy = a, — x =10, by = by — x =0. 

Pasul al şaselea: obligatoriu vom pune x, = 10. 

S-a obținut o soluţie admisibilă de bază: 

xi 35, X = 0,33 = Gais 0, Xa E 15, 39 2103 = 153, =0, 

x31 = 0, x3 = 0, x33 = 10, x34 = 10 


costul total al transportului fiind 300 u.b. 
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b) Metoda costului minim. Vom considera din nou tabelul simplificat: 


40, 30, 0, 


20, 0 


20 0 0 0 


T Pasul întâi: mincj=c = cp = = else 3: Putem dvd a oricare 
dintre componentele x), te Mr plc Luărn Xn =min(a2,b2)=10, X2 = X32 = 0 
(componente nebazice). roca pe a, şi b, cu a, =a, — x = 30, b; =b; — xy =0. 

Pasul“ al doilea: . min Cij= C= e ipe 2. Luăm  xp=min(a5,b,)=30, 
X23 = X24 = 0 (componente nebazice). Înlocuim pe a, şi b, cua/=a'—x,, 5 0, 
b =h- xn =20, 

Pasul al treilea: min c;;= cy3= c35 2. 

Luăm Daaa min(a,b,) = 25, X33 = 0 (componentă nebazică). Înlocuim a, 
şi b; cu a, = ay =% =10, b? =b -x =0. 

Pasul al patrulea: min C;j= Cu4= Cs Ca3= 3 . 

Luăm x, = min(a3,b) = 20, x, =0 (componentă bazică), ENE! 
(componentă nebazică). Înlocuim a, şi b! cu 3 = 03 — x4 30, b'=b—x, =0. 


Pasul al cincilea: obligatoriu x 4310! 
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Am obținut pentru problemă o soluție admisibilă de bază degenerată 


(x, = 0 componentă bazică). Ea este indicată în tabelul de mai sus. 


Costul unui asemenea transport este 220 u.b. Metoda costurilor minime dă 
de cele mai multe ori o soluţie admisibilă de bază mai bună decât metoda 


diagonalei, în sensul că realizează o valoare a cheltuielilor de transport mai mică. 


3.3.2. SOLUȚII OPTIME PENTRU O 
PROBLEMĂ DE TRANSPORT 


Soluţiile optime (soluţii admisibile ce realizează minimul funcţiei obiectiv) 
ale unei probleme de transport le vom căuta în mulțimea soluţiilor admisibile de 
bază nede generate. Al goritmul pe care îl von. prezenta presupune aşa ceva. Asta nu 
înseamnă că o soluție admisibilă de bază degenerată nu poate fi soluţie optimă a 
problemei. Să considerăm că am determinat pentru problemă o soluție admisibilă 
de bază, fie ea x(x) . Pornind de la aceasta, să vedem cum putem obține o soluție 
admisibilă. = $ 

Convenim ca o pereche de indici (i, j) să o numim celulă. Când vom spune 
celulă liberă înțelegem că ea se referă la o componentă nebazică, adică nu se 
găseşte printre cele m+n-1. Celule ocupate vor cele care se referă la componentele 


bazice. Să considerăm o celulă liberă (i, J). 


Definiția 3.3.1: Un şir de celule care începe cu celula liberă (i,j) şi se 


termină cu aceasta, conținând în rest numai celule ocupate şi anume câte două din 


aceeaşi linie şi coloană se numeşte ciclul celulei libere (i, j). 


Cel mai simplu ciclu al unei celule libere (i, j) se prezintă sub forma de mai 


jos, unde în dreptul celulelor am scris şi componentele corespunzătoare ale soluției. 
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(Lp (E 
cu ARE Ai 


Să atribuim componentei X; (componenta nebazică), o valoare 6 >0.Vom 


obține atunci un vector care are m+n componente diferite de zero, acesta trebuind 
să verifice sistemul de restricții. 

Acesta din urmă este satisfăcut dacă vom scădea şi aduna succesiv valoarea 
9 la componentele care au intrat în ciclul celulei libere (i, j). Indicăm pe acelaşi 


ciclul acest lucru. 


Gi e Gi Ji) 


E E 


Pentru a fi verificate condițiile de nenegativitate va trebui ca 0 să fie luat, 


astfel încât x; —0 >0 şi Xij 020. 
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Odată îndeplinite aceste condiții, putem spune că am obținut o nouă soluţie 


admisibilă x(x,) unde: 


x =0, ue —9, Bz, 1-0 


-ni (3.2.2) 
Xij = 3,50 Xe = Xp pentru zii; j 2 j j 


Dorind acuma ca noua soluție admisibilă să fie şi admisibilă de bază, este 


suficient să considerăm: 

0 = min(x Xj) (3.3.3) 

În acest fel se determină care componentă bazică devine nebazică. Relaţia 
(3.3.3) reprezintă forma criteriului de ieşire din algoritmul pentru găsirea soluțiilor 
optime a unei probleme de transport. 

Să vedem acum care sunt condițiile ca soluția admisibilä de bază x să fie 


soluție optimă. să evaluăm diferența f(x )- f(x) } 


Vom avea: 


m n 


f(x)-f (x)= E Lre le — Ñe) = =c; 0+c,8—c,0+c,j8 


zi că 
sau 

(x)-f(x)=68, (3.3.4) 
unde 

> ey Cu Cap FC“ | (3.3.5) 


care reprezintă suma algebrică cu semn alternat a costurilor celulelor 


corespunzătoare ce au intrat în ciclu şi care începe cu —C;j. 


Din (3.3.5) se observă că dacă pentru toate celulele libere (i, j) vom avea 


ô; <0, şi aşa cum 0 >0, soluţia admisibilă de bază x(x,) considerată reprezintă 
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soluţia optimă a problemei de transport. Rezultă în felul acesta că criteriul de 


optimalitate este 


<0 | (33.6) 


pentru toate celulele libere (i, j). 
Tot din (3.3.4) se constată că dacă există 6, > 0» soluția x poate fi 
îmbunătățită; Pentru aceasta fie: 


dal max PE, (3.3.7) 


3,20 
Atunci componentei nebazice x, îi vom atribui o valoare 0>0. Relaţia 
(3.3.7) reprezintă forma criteriului de intrare în algoritm. naprije 6, joacă rolul 


diferențelor z, CE, din pr oblemele de kaa d, N de care ne-am ocupat. 


Sintetizând cele de mai sus, în practică aii în felul următor: 


a) Se determină o soluție admisibilă de bază nedegenerată prin una din cele 


două metode prezentate. 


b) Pentru toate celulele libere se calculează cantitățile 5, 


i daciei toate Calitățile d, e mai mici sau egale cu Zero, atunci 
soluţia este optimă; 
- dacă există 6, >0 se trece la etapa următoare. 
c) Se aplică criteriul de intrare. Vom atribui o valoare 9>0 componentei 
X; (nebazică) cea pentru care avem îndeplinită relaţia (3.3.7). 
- d) Se: consideră ciclul celulei Xj» se introduce 0 în ciclu şi se aplică 


criteriul de ieşire. Se determină valoarea lui 8 aup relația (3.3.3) indicându-se 


deci care componentă bazică devine nebazică. 


e) Etapele algorimui se reiau până se găseşte soluția optimă a problemei. 
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Observații: 


a) 


b) 


d) 


e) 


Prezența semnului de egal în criteriul de optimalitate şi de intrare indică 
existența mai multor soluții optime. i 

Pe parcursul aplicării algoritmului pot să apară soluții admisibile 
degenerate care pot conduce la un fenomen de ciclaj. Pentru. evitarea 
apariţiei soluţiilor admisibile degenerate se întrebuințează aşa-numita 
metodă a perturbării. Aceasta constă în înlocuirea cantităților .a, cu 
a, +€, şi a unei cantităţi b; cu b; + me (e o cantitate alea pozitivă). 


Când se găseşte soluţia optimă se face € =0. 


Dacă la aplicarea criteriului de intrare sunt mai multe mărimi õi pozitive 
maxime egale, putem considera oricare dintre acestea putând totuşi uneori 
să conducă la prelungirea timpului necesar găsirii soluţiei optime. 

"m 


Psobl&mele de ara = în care Sa, ze Sb, (probleme de transport 
i=l j=l 


neechilibrate) pot fi. aduse la forma echilibrată prin considerarea unui 


centru fictiv sau a unui depozit fictiv în care să se găsească diferenţa de 


cerere sau de produs, costurile corespunzătoare fiind nule. 
Modelul problemei de transport poate fi considerat şi cu cerința de 
maxim pentru funcția obiectiv cu modificările: - 


= soluţia este gáis dacă ô; > 0, 


— criteriul de intrare devine õi; = min ARĂ 
Sue <0 
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Exemplul 3.3.2. Să se determine soluţiile optime ale problemei de transport 


din exemplul 3.3 

Pentru a fi siguri că nu apar soluții admisibile de bază degenerate, aplicăm 
de la început metoda perturbării, adică considerăm: 

ay =35+8, a3=40+6, a3=20+e, b, = 50, b, =10 b} = 25,4b, = 10,+3e 


O soluţie admisibilă de bază determinată prin metoda costului minim din 


tabel este dată în tabelul: 


Să găsim cantităţile $, pentru celulele libere. Considerăm ciclul celulei 


libere (1,1) şi în dreptul celulelor scriem costurile. 


4 3 


(1, 4) 


AI ee (3.4) 
3 că 
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Vom avea: õu =-4+3-3+3=-!1 


Analog: 
õn =—2+2—34+3—3+2=-—1 
õn = —4+2-3+3—34+2= -3 
62, =—5+3-3+2= -3 
ô =—4+2-3+3=-2 

„aa =—24+2—34+3=0 


Deoarece toate cantităţile 6, sunt mai mici sau egale cu zero, am găsit 
soluția optimă a problemei. Făcând £=0 găsim: 

Xu =0, xp = 10, x, =25, x4 = 0, X= 40, x = 0, X3 = 0, x4 =0 

X = 10, x4 =0, x= 0, Aa = 0 Ax) = 210 u.b. 

Deoarece 5 =0 , algoritmul poate continua atribuind componentei x33 


(nebazică), o valoare 0>0. Considerăm ciclul celulei (3,3) şi introducem pe 0 în 


ciclu pentru a satisface restricțiile. 


25—0 £ pe ned 


CE TR (1,4) 
(3,3 je— (3,4) 
9 10+ 2—0 


Aplicăm criteriul de ieşire luând 0 = min(25, 10+ 2e)=10+2s. 
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Componentele noii soluţii sunt date în tabelul: 


Făcând £=0 găsim soluţia optimă. 
Xii =Q. Ap = 10, X13 = l Xi4 = 10, X2] =40, X22 = 0, X23 = 0, Xag =() 
X% =10, x =0, PN = 10, x34=0, (min)f(x)= 210 u.b. 


Exemplul 3.3.3. Să se determine soluțiile optime ale problemei de 


transport: 


Problema nu este echilibrată. Vom adăuga un depozit cu costuri zero. 
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Vom avea deci: 


O soluție admisibilă de bază determinată prin metoda costului minim din 
tabel este prezentată în: 


Componentele soluţiei au fost determinate în ordinea: 
%1 = 15, Xp = 25, x = 25, Xp = 20, 
xy > n X1 > 0, 2 =0, 32 = 0, X33 =0 


Pentru celulele libere calculăm cantităţile NAN 


AEE, WES E IE ES 
8n = —4+2—3+1=—4 
ôy =—0+2—1+0=1 
ô =—0+3—1+0=2 
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Soluţia nu este optimă. Avem că ôn = max 6;j=2. Considerăm ciclul 
33 J 


6,0 


celulei (3,3) componentei x, (nebazică) îi atribuim o valoare 0 >0 Şi introducând 


pe 0 în circuit avem: 


25 + 0(30) 5 —0(0) 


(bi) e (LA 


TE | POS ee S-A 
15—0(10) 8(5) 


După criteriul de ieşire luăm 0 = min(,15)= 5. Noua soluţie este 


prezentată în tabelul: 


Pentru celulele libere vom avea: 


Oi 2341-0402 29 
SEEEN 1= 0-0 —2 

da = —4+2—1+0—0+1= 32 
533 = —0+2-—1+0=1 
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Soluţia nu este optimă. Atribuim componentei x; o valoare 9>0.: 


30+0(40) 20 = 010)... 


(1, JĀ (1;2) 


GDE 
10 — 0(0) 0(10) 


Noua soluție este prezentată în tabelul: 


Această soluție este optimă deoarece avem: 
fa lb = 36 336 E 
Soluţia problemei iniţiale va fi: 
| Xu = 40, x2 = 10, x3 =0, x =0, xx ZY, N 25. (min)f(x) = 85u.b. 


Se observă că cerințele datelor C2 şi GC nu vor fi satisfăcute î în totalitate. 


Capitolul 4 
SERII NUMERICE. SERII DE PUTERI 
——— Arena 4.1.1 (ie carocierizire d iririley come) O OOOO 


Este bine cunoscut faptul că, în prezent, fenomenele economice devin din 
ce în ce mai complexe şi interdependente; mai nou, tendința de globalizare a 
economiei (şi nu numai) duce la o concurență acerbă pe mai toate pieţele şi în mai 
toate domeniile de activitate economică, fenomen care se simte din ce în ce mai 


pregnant şi în economia românească. 


Aceste situaţii cer din partea specialiştilor în economie, cunoştinţe cât mai 
precise în vederea observării, înțelegerii, modelării şi rezolvării pe baze riguros 
ştiinţifice a problemelor complexe cu care se confruntă. Această abordare nu mai 
este astăzi posibilă, fără recurgerea la un aparat matematic cât mai adecvat. 

Modul mai puţin riguros, chiar empiric, de luare a deciziilor economice doar pe 
baza unei simple experienţe, este total depăşit, astfel încât cei care nu se vor adapta 
noilor condiții şi cerințe economice, vor fi eliminaţi încet, încet de pe piaţă. 

În: acest capitol se prezintă, într-o situaţie accesibilă studenţilor economişti, 
câteva elemente de bază din teoria seriilor numerice şi a seriilor de funcții (cazul 
particular al seriilor de puteri). Aplicațiile practice ale acestei teorii se întâlnesc în 
multe domenii de cercetare din matematică, fizică, chimie, biologie etc. În particular, 
aplicaţiile cu caracter economic ale teoriei seriilor sunt multiple şi permit modelarea şi 
rezolvarea unor probleme economice cu “caracter. atât determinist, cât şi aleator 


(probabilistic). 
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Elemente şi aplicaţii ale teoriei seriilor sunt regăsite în Teoria 
probabilităților şi Statistica economică, Micro} şi Macroeconomie, Gestiunea 
stocurilor şi a fenomenelor de aşteptare etc. 

Pentru parcurgerea şi înțelegerea noțiunilor prezentate în acest capitol sunt 
necesare noţiuni de teoria limitelor şi a Şirurilor studiate la Da DM Matematică din 


clasa a XI-a. 


4.1. SERII NUMERICE 


4.1.1. ȘIRURI DE NUMERE REALE 


Definiţia 4.1.1: Numim şir de numere reale, o aplicaţie: 


F:N—R 
fhm ali ei 


Vom nota şirul cu simbolul: (a, en 3an SE numeşte termenul general al șirului. 


» 


Definiția 4.1.2: Spunem că şirul (annen este convergent- către un număr 
a ER, dacă orice vecinătate à lui a conţine toţi termenii şirului cu e xcepția unui 


număr finit de termeni. 


“Vom nota cu a —a sau [ic =a: 


n=oo 
Numărul „a“ se numeşte limita șirului şi dacă există este unică. Dacă şirul 
(avem Nu este Convergent, atunci acesta se numeşte şir divergent: Un şir 


divergent sau nu are limită sau limita sa este +oo sau +00. 


Exemple: 
“DĂ n? ÎI, reda 
a) a= i PS 3 Ca =|l+—|. sunt şiruri convergente 
n +l 2n"+3 Na ea: dt 


(cu limitele 0, — respectiv e=2,71...; 


5 |—= 
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Khi TaT 3 > ră ă ` j : k 
b) a, = Vn bi CPAN ; Ca =N] =n sunt şiruri divergente cu-limitele 
n 
+00, +00, respectiv —00; 


c) a, =(—1)", b, =(—1)" n(n? +n +1) „Ca =sin(n +1) sunt şiruri 


divergente care nu au limită. 


“Teorema 4.1.1 (de caracterizare a şirurilor convergente) ` 
Şirul a, —a dacă şi numai dacă: (V)e > 0,(3)n. e N astfel încât 


u(Y)n>n,.. 


Teorema 4.1.2. (Cauchy) 


Șirul (a, „em este convergent dacă și numai dacă: (Y)e > 0, (3)n, EN 


astfel încât : japik —a,|<e, pentru (V)n >n, şi (Y)k EN. 


Observatii: 

i) Cele două teoreme, ca şi definiția, nu sunt folosite în mod curent în 
aplicaţii practice; cu ajutorul lor se pot demonstra rezultate a căror 
aplicare este mult mai simplă (de exemplu, teorema de convergenţă a 
şirurilor monotone şi mărginite, teorema cleştelui etc.). 

ii) Aşa cum se observă, aplicarea teoremei 4.1.2. nu presupune 
cunoaşterea apriorică a limitei a, în enunţul acesteia intervenind numai 


termenul general a, respectiv a, ; acest lucru este deosebit de util în 


aplicaţiile care conţin Suzi recurente sau în care limita şirului este 


dificil de deteriminai. 


114 MATEMATICI APLICATE ÎN ECONOMIE 


4.1.2 NOȚIUNEA DE SERIE NUMERICĂ. CONVERGENȚĂ. 
PROPRIETĂȚI GENERALE ALE SERIILOR 


După cum se ştie, oricărei mulţimi finite (a aces] de numere reale i 
se poate asocia un număr real S=aytazt...ta, numit suma acestora. Se pune 
întrebarea firească: dacă mulțimea conţine o infinitate de numere reale 
ga Bega ed goe) (cu alte cuvinte un şir de numere reale, (a, )uem.) i-am putea 
asocia un număr S€ R care să generalizeze conceptul de sumă finită? 


Kie sirul: (2, ennio A aE 


Vom asocia acestuia, un ¿noù şir de numere reale, 


Ereni So» SS eSa definit astfel: 


So = ao 
Sya ayti 
S, =ap tata» ir (Ala) 


n...canoecsreasseseonssese 


Si Fapta trap P.ta, = dar 


pe care îl vom numi şirul sumelor parțiale ataşat șirului (a, )nen -" 


Definiţia 4.1.3. Numim serie numerică, de termen general a, expresia 

def Ri | 
a. = lim S, =agtaj +a, +.. P EIE (4.1.2) 
n30f 9 ae E 


Observație: In unele cazuri numerotarea termenilor seriei poate începe de 


la n= Žan) sau de la n=ng oarecare fixat (+. (vezi şi proprietăţile 
n=] ; 


n=no 


seriilor numerice). 
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[e9] 
Definiția 4.1.4.: Spunem că seria Y a, este o serie: 
n=0 


a) convergentă (C ), dacă şirul (S, en este şir convergent; 
„„b) divergentă (D)., dacă şirul (S, )ney este şir divergent, 


Cu alte cuvinte problema convergenţei unei serii numerice este echivalentă 


cu problema convergenţei unui şir de numere reale şi anume, şirul sumelor parțiale 
(Sa en = 
Observaţii: 


0o 
i) Dacă S, —S, vom spune că seria Y` a, (C) şi are suma S (notăm: 


| j. lai l ' 
ii) Dacă S, — +00, spunem că seria $` a, (D) şi are suma to notăm 
n=0 


a = to); 


=) 
ll 
© 


(8.9) 
iii) Dacă (S, aer nu are limită, spunem că 5 a, (D) şi nu are sumă. 
n=0 
In studiul seriilor numerice, trebuie cercetat (Ca şi în cazul şirurilor) două aspecte: 
- natura seriei (dacă este convergentă sau divergentă); 
- suma seriei (în caz de convergenţă). 


„Dacă pentru şiruri stabilirea convergenței şi determinarea limitei acestora 


este în cele mai multe cazuri una și aceeaşi chestiune, în cazul seriilor situaţia se 
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schimbă radical; de cele mai multe ori este relativ uşor să stabileşti convergența 


seriei, dar foarte greu, ba câte odată chiar imposibil de determinat suma acesteia. 


De obicei, cu metodele specifice analizei numerice, se determină o 
aproximaţie S’ a sumei S a seriei (S'=S). `~ 

Pentru a se cerceta natura unei serii, nu se poate folosi întotdeauna definiția, 
ci adesea se utilizează condiții echivalente cu aceasta sau numai condiţii suficiente 
de convergență numite „criterii de convergenţă“. Vom prezenta unul dintre aceste 


criterii în cele ce urmează. 


Teorema 4.1.3. (criteriul general de convergență a lui Cauchy) 


Seria Sa, este (C) dacă şi numai dacă: (Ve > 0,(3)n, EN astfel 


wi 


încât: lags Ea ++ anik] <e, pentru (Y)n > n., (Vk EN. 
Această teoremă se obține ca o consecință a Teoremei 4.1.2, considerând în loc 
de şirul (a, nem» şirul sumelor parțiale (S, en - Într-adevăr, conform (4.1.1) avem: 
Sue —S,| =i|(aj +4.art any Eta) A +... +a, )| == 


= |ant +-+ ans] 


4.1.3. EXEMPLE REMARCABILE DE SERII 


n=0 


1) Seria geometrică S aq”; cu aqER| 


Termenul general al acestei serii are forma a, = aq” (*) şi depinde de doi 


parametri (constante reale) a, q. Evident, este logic ca natura seriei şi suma acesteia 


în caz de convergență să depindă de aceşti parametri. Astfel, dacă: 
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i) „a = 0, atunci a, =0, (v)qgehR, deci seria este convergentă şi are 
suma S$=0. 
il) az0, avem: 


(*) = 
Sa = ao +a; +a; +...a, =a Haq Haq’ h... Haq” 


| 2 
=a(l+q+q"+...+q") 
deci S, reprezintă suma unei progresii aritmetice cu rația q (de aici şi 
numele seriei) şi conform rezultatelor de la acest capitol (vezi manualul 
de Algebră, clasa a X-a) vom avea: 


Dă n+l 


a 
S= 1—q 


n 


> ql 
amt) galne 


Trecând la limită în relația de mai sus, obținem: 
—— ; gel) 
I—q 
lim S, = 12004 qe[l,00) “(depinde de semnul lui a) 


n= 


A) ; qe(—o0,—1] 


cu alte cuvinte, pentru: 


qE (-11)> Să gg = —— p converge cu suma S= an 
n= lmg l=q (4.1.3) 
qe(-o,L]U[L,oo) = a aq"(D) (cu suma S = +œ sau nu există ) 
n=0 


Observaţie: De fapt, această serie cuprinde o infinitate de serii particulare 
(pentru a, q având valori fixate) a căror natură o cunoaştem; foarte important este 


că, în caz de convergență, ştim şi suma seriei. 
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Exemple: 
e | (1) fst e 
a) 23- =} |z| „este serie geometrică cu a =], q=—şiiconform 
îti 2? n=0 s ' 25 
(4.1.3) rezultă că: 
& sa] Pare | 
—=l+- +-+... +—+...=2 (C 
er 2 2% P ; ) 
e TE ETT | ama ba 
b). = 2i „ este serie geometrică cu a= 1 şi q= a$ ŞI 
n=0 4/0" n=0 


conform (4.1.3) i că: 


y 3 
a (C) 


n=0 
5 co (6 | 5 
6) tp =) a) este serie geometrică cu a = 5, q => şi conform 
Da 2 3 n=0 4 
+, pa n+l 
(4.1.3) seria diverge e fina acest caz J` = +00]. 


n=0 2 
lee 5 Gin (v2); = (a D(- -V2) 'este serie geometrică cu 
n=0 n=0 


—1, q= — V2 şi deci seria este divergentă; mai mult suma nu există. 


2) Seria armonică generalizată Sa ER 
usl n 
Evident, ca şi în şirul seriei geometrice, natura (şi în caz de convergență 
suma) seriei armonice generalizate, va depinde de valoarea parametrului œ 
Aplicând, Teorema 4.1.2 se poate arăta că seria: 


| (EN o ihai 
| È 5 (4.1.4) 
n=} N 


I 
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Observaţii: 

i) În caz de convergenţă (a > 1), nu există o formulă generală de stabilire 
a sumei S în funcție de parametrul a, asa: cum s-a stabilit în cazul 
seriei geometrice; 

ii) Evident, în cazul în care seria diverge (œ <1), suma seriei este 


S= +œ (ca serie cu termeni pozitivi). 


Exemple: 
cond E ce] erp n 
a) '5, =D) deoafece o E] POS == P-n H0) 
n= N n=} N PAS, AN du 
i o ] - or, ] l | l 
b) 5 — (C), deoarece a=2>1 4 me el AT S a 0: i ae Ş 
n=1 Nn” n=} N 2 
e = Bi œ ] 
c) = dag (D), deoarece a= <1 
n=l Jn? n=l 3 
n: 


Seria BE — se numeşte seria nitel e (simplă) sau seria lui Riemann, 
n 


n=l 
deoarece termenul general an este media armonică a vecinilor săi, adică satisface 


t "A 
+ 


Vi 2 
relația Era 


oo 1 
3) Seria telescopică 
i 5 i) 


Această serie este convergentă. Într-adevăr, scriind termenul general 


F per] 
sub forma (*)a_ =————, obti t ] l af atul 
IE ma, Ei i v m general al șirului 


sumelor parțiale CEIR de forma: 
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Sa Sa ta, +...+a lipse, g E Laza 
d mi ae Mei 


+ f EPE | = (după simplificări) 


A 2] 
l 
= 1 
n+l 
Atunci Ș lim S. = limaa A SR LO 
tunci 9” = lim = um|l———|, deci = A 
n n(n+l) n=% " n-o n+l m n(n+l) 
s y l i Tie ] 
Dacă notăm cu A, = —, atunci termenul general al seriei a, = , se 
i m n(n +1) 
poate scrie sub forma: 
ap E Op Oi | | 5 (44.5) 


Seriile numerice a căror termen general a, au forma de mai sus, :se numesc 


serii telescopice și se poate observa uşor că seria dd (EN (0, hai (E). 
n=l 


Mai mult, seria va avea suma S=0, —Q, unde œ, este primul termen al 


șirului (Œ due» iar œ= lim 0, i 


Vom enunța în continuare câteva din proprietăţile generale ale. seriilor 
numerice, care se pot demonstra relativ uşor, ținând cont de definiţia seriei şi a 


operațiilor cu şiruri convergente. 


4.1.4. PROPRIETĂȚI GENERALE ALE SERIILOR NUMERICE 


ate 5) | 
(P1) Dacă unei serii Ş: an „i se adaugă sau i se suprimă un număr finit de 
n=l ' 


termeni, natura acesteia nu se schimbă. 
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Observaţie: Evident, în caz de convergenţă, suma S a unei serii se modifică 
corespunzător prin adăugarea (sau scăderea) sumei termenilor finiţi care se adaugă 
(sau se suprimă). 


(P2) Sacă unei serii convergente Y a,» i se asociază termeni în grupe 
n=l 


finite, cu păstrarea ordinii lor, atunci se obţine tot o serie convergentă şi cu 


aceeaşi sumă. 


Observaţii: 


i) Proprietatea de mai sus, poate fi scrisă şi sub forma: 

| Au E | . 95 

(P,) Dacă ăn =a,taz+...+a, ps, atunci Şi: 
n= 


(a, +a, +... Han, ) H Cangi at ap) (a. a Ea ) e. = S 


ii)  Asociind termenii unei serii divergente, în grupe finite cu păstrarea 


ordinii, se pot obține serii convergente. 
Exemplu: 


Fie seria $ (—1)"=1—1+1—1+...+(—1" +... (termenul general fiind 
n=0 


a, = (— 1)" ). Atunci: 
Sg si SĂRI ak0, S, 1-1412, SETE 


ge... 


deci: Sa, =1 şi San =0; evident şirul sumelor parţiale (Sa nea diverge deoarece 


are două subşiruri care converg la limite diferite. Rezultă că seria 5(-1)" (D) ŞI 
i n=0 


nu are sumă. 
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Dacă asociem termenii în grupe de câte doi (sau de câte patru, şase etc.) 


obținem: 
(1 —1D)+(1-D+ (1-D+...(1—-D+...20+0+...+0+...=0 


deci noua serie obținută converge şi are suma S = 0. 


oo 09 
(P3) Dacă seriile X a, şi E b, sunt convergente cu sumele S}, respectiv, 
n=] n=l 


2, atunci seria > (a, +b,) converge şi are suma S; + S>. 
n=l] 
Cu alte cuvinte dacă: 


aa 3 F (acA S, +S, (4.1.6) 
: i n=l. i , 


(P4) Dacă NER , atunci seriile Ş: a, şi È àa, au aceeaşi natură; mai 
n=] 


mult în caz de convergenţă, dacă Ywa, =S atunci Y. Aa, BAS 
n=] n=] ; 


Teorema 4.1.4. Dacă seria S a, (O), atunci lim a, =0. 
n=1 n00 


Demonstratie: 


E * . 
Fie Sp Sayta2 +... tan+(V)neN . Deoarece seria converge, avem: 


e 


= lim, =S (|) 


A n= 


n 


n 


"Dar a, = EA =S,-13 trecând la limită în această relație obținem: 


lima, = lim (S, — Sp) > lim S, — lim S, Tie -S=0 


n—yeo n=» n—>oo 


q.e.d. 
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Observaţie: Condiţia din teorema de mai sus (a, —>0) este o condiţie 


necesară, dar nu şi suficientă pentru convergența unei serii; deci, dacă termenul general 


a, —0 acest lucru nu implică obligatoriu că seria a, (C). Ca exemple se pot da 
n=l 


; j N ) TA At! ; l 
cazurile particulare ale seriei armonice. Astfel 3 — (D) chiar dacă a, = — —> 0. 
n=i N n 
O formulare echivalentă a teoremei 4.1.4. des folosită în aplicații practice; este: 


Teorema 4.1.5. Criteriul general de divergență (C.G.D) 


Fie seria D ano a, ER. Dacă lima, 20 atunci seria 5 a, (D). 


n=1 1—00 n=1 

Exemple: 
a) S -is (D), deoarece lima, = lim 4 ai, 

n=1 2n neo n>» 2n+] 2 
b) ŞS (=) (D), deoarece: 

n=] n 

n n 
lira 2 liră (=) SA +4) PEETI, 
n—oo n—»oo n n—>oo n 
IPE Yhn 7 | l+n— | 

c) $ ———~— (D), deoarece lim a, ca, ju — = -xv z0 

n=l 2n l n= n= ao s 


În concluzie, dacă: 


| | j i kaSi 
Y lima, 20>% a, 


n=] 


ii) lima, =0 = poate fi (C) sau (D) (trebuie făcută o analiză 


n= 
n=l 


suplimentară în acest caz). 
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4.2 SERII CU TERMENI POZITIVI 


Definiţia 4.2.1. Numim serie cu termeni pozitivi, o serie de forma: 


8 


în cu 2,20, (V)neN. 


naj 


u 


Teorema 4.2.1. (de caracterizare a convergenței seriilor cu termeni pozitivi) 


Seria cu termeni pozitivi "Vans an 20 converge dacă şi numai dacă 
n=l 


şirul sumelor parțiale (S) en este mărginit. 


Demonstrație: 

Avem 

Sns TSn = (aptag feet an) A a aSa 20. 
deci Sas ZS, + adică şirul (S, ae este monoton crescător. 


Dacă (Sem este şi marginit (superior), rezultă că este convergent 


(conform teoremei de convergență a şirurilor monotone şi mărginite), deci Xa, (C). 


n=} 


n=i 


Reciproc, dacă Sa, (C) = (S, nem (C). Dar un Şir convergent este 


obligatoriu mărginit, deci teorema este demonstrată. 
q.e.d. 
Observaţie: În aplicațiile practice, Teorema 4.2.1. nu este folosită foarte 
des deoarece, în general, este relativ complicat a se demonstra mărginirea şirului 


(Sa nem s într-adevăr, trebuie demonstrată relația: 


|S,|=|a,+a,+...+a, (VY)neN' 


cu M > Oconstantă care nu depinde de ne N". 
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Mult mai utile sunt următoarele rezultate (numite criterii de convergență), 


care pot fi clasificate în două categorii: criterii de comparaţie şi criterii cu limită. 
1) Criteriul de comparaţie de specia I (cu inegalităţi simple) 


| 35 ü | 
Fie seriile: Oa, şi > b, care satisfac relaţiile 0 < a, Cba, (v)neN'. 


n=l n=] 


Atunci: 


i) Dacă Sb, (O > Sa (O) 
n=l 


n=l 


i) Dacă Sa, (D) > Sb, (D) 
n=] n=l] 


A „8 dl 4 DE 
Exemplu: Fie seria Y —. Termenul general a, = — satisface: 
n n n 


n=I N n 
28S rmm “INS i TT Z M pot 
0 <a, = L.L. -<= 2=2 =b, 
n-n:n..0/ nen n n nn n 


Dar: 
© % 2 æ ] : : A 
Lb = 43 (C) (ca serie armonică generalizată 
n=l n=1 H n=i N iii i 
œ |] 
Dute ag a I a D) 
n= N 


an Aaa ai a 
Atunci conform punctului i) din criteriu N Tao): 
; n=I N 
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II) Criteriul de comparație de specia a III-a (cu limită) 


CO CO 
Fie seriile Sa, şi Sb, cua,, b, >0.Dacă există: 
n=] n=l na: ' : | i 


PB 
lim — = k € (0,00) (4.2.1) 
n-o b, stea 
atunci seriile au aceeaşi natură. 


oo l l 
Exemplu: Fie seria ———— cu a =———>0 
i > 2n+3 PRF 


l 
im 202 ir 20 Ba ag Monod creeate 1 e (0,%) 
n=cob, n= d i n=o2n+3 2 2 

= | 


(9,9, (2,9) 
deci conform criteriului, seria "Y a, are aceeași natură cu seria 3) b, , deci diverge. 
n=] n=] 


Observații: 
i) Criteriul de comparație la limită, tratează şi cazurile când k=0 sau k = +00, 
"dar s-a renunțat la prezentarea lor pentru simplificarea expunerii. 


ii) Aplicarea criteriilor de comparație, nu este în general foarte simplă, 


deoarece trebuie satisfăcute simultan următoarele două condiții: 


a) alegerea seriei de comparație (de exemplu, D b, ), funcțe de seria 
n=l 


i OD 
dată (de exemplu, 3-a, a cărei convergență vrem să o stabilim); 
n=l 
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À CO 
b) trebuie cunoscută aprioric natura seriei 2- b, (convergentă. sau 
l gis 


divergentă). 

Prima condiție este influențată foarte mult de abilitățile matematice ale 
rezolvitorului, iar cea de a doua de „experiența“ sa. Pentru această din urmă 
condiţie, seriile geometrică şi armonică pot furniza o infinitate de serii particulare 
(a căror natură o cunoaştem) posibil de folosit ca serii de comparație. 

iii) În schimb, următoarele criterii (numite criterii cu limită) se aplică mult 

mai uşor, formularea acestora depinzând exclusiv de termenul general 


OO 
al seriei = du a cărei natură vrem să o determinăm. 


=l 


III) Criteriul raportului (D`A lembert) 


ST i j . Se i ` . . ! a 3 i j . 3 v l 
Fie:setia => a, » a, 20 Soim AL S Ne [0,0]. Atunci, dacă: 
n=] N—OO dn 


)X<1l3 Sa 
n=l 


A> SaD) 


nal 


iii) =1 ; caz de nedeterminare (nu se poate preciza natura seriei) 


Exemple: 
1 ETN n! i n+l)! ; 
D &—> ;cu a, =— respectiv a SED g Atunci 
n n n n+l n+l 
n= N n | (n +1) 
„di PD E a nMn+l "i 
À = lim = = lim (Ea, — = ime BE E a Tahi 


noe ap  n>=(n+1)"*! n! n>e(n+1)" wD BI 


128 MATEMATICI APLICATE ÎN ECONOMIE 


n 
CATA Seria 72 fan l n! 
(după simplificări)= lim Ra BA deci seria S — converge. 


ta n=} n’ 

09 Dn n i Dn! 
2) = ‘aCi E ŞI a&i =. 

2 n“ +l mri oi su (n+1) +1 

Atunci: 
n 2 | 
z " l w: 
ASM ba Eo BEE pEr iia ii dizol n hheag e But anuar pla | 
n= aj nson? n42n+2 2 Pa ebay? 


S n 
deci, seria >, ri aa , diverge. 
nalial, 


3) $ TT papy avei. . Avem: 
aci N “pi n n+l Pei : 


ne A ; n 
lim 4 = lim — =], 
noe a, nae pi- 


l 
deci conform acestui criteriu nu putem preciza natura seriei armonice S — (despre 


l n=] N 
care ştim însă că este divergentă, vezi formula (4.1.4)). 
IV. Criteriul rădăcinii (Cauchy) 
e. e Da 
Fie seria J` a, ap >0 şi lim Y = Àe [0,%]. Atunci, dacă: 


n=] DS 


DLE a, (0) 
n=l 
i) STS Sa (D) 


n=] 
II)Ă = 1= caz de nedeterminare (nu se poate preciza natura seriei) 
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Exemple: 
o n? n? saj 
1) XS —s cu a, =— >0. Atunci: 
PEES yai 


n=] 2n 


= = lim DE lim. af i lim i ala, =5<1 deci, © seria geni aT 


converge. 


2) = aE ; Cu ah 


n=] 


| 


2 
n 
i >0. Atunci: 
n ' 


EFA PETTINE TA T r 
> adi a = [i = | = lim |1 =e >] 
A: seria a diverge. 
Sni l 
3) £ — ;cu a, = 320. AUNO 
nzi n" n 


= lim n gan = = lim mi = a ban, deci nu putem preciza cu 


n—yco n—e0 sa JA ry 


acest criteriu, natura seriei a i (despre care Ştim că este convergentă, ca 
n=l N 


seria armonică generalizată, i acu a Re. Sd) 
n=] n“ 


Observaţii: 


i) În EA din T de de mai sus am folosit limitelă Fitgteinvezilale: 
lim Vn =1 şi lim VP, (n) =1 42.2 


unde: P,(n)= a,n" + ani! +...+ an? +an-+ag polinom de grad 


k cu coeficienţi reali şi a, >0 
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ii) 


iii) 
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Se observă că aplicarea acestor criterii constă doar în calculul unei 


limite; totuşi dacă termenul general al seriei a, are o formă mai 


complicată, determinarea lui A poate fi o problemă deloc uşoară. Se 
poate apela în aceste cazuri la criteriile de comparaţie sau la alte criterii 
de o formă mai complicată care nu sunt prezentate aici (se poate 
consulta în acest scop orice curs de analiză matematică care tratează 


convergenţa seriilor numerice). 


Datorită relaţiei: “lim în 27 lim g a, (4.2.3) 
n-o a, € n>% 

rezultă că dacă un criteriu ne conduce la cazul de excepție A =, nu se 
încearcă aplicarea celuilalt criteriu, deoarece se va obține ‘aceeași 
concluzie; obişnuit, aplicarea criteriului raportului sau rădăcinii se face 
ținând cont de forma termenului general a, (astfel încât calculul limitei 
să fie cât mai simplu). 

Se observă, că există situaţii — chiar pentru serii foarte simple — în care 
aceste criterii nu pot'preciza natura seriilor (1 = 1). De obicei, în aceste 
cazuri se aplică criteriile de comparaţie, respectiv, criteriul general de 


divergență (C.G.D), sau următorul criteriu (a lui Raabe-Duhamel). 


V) Criteriul lui Raabe — Duhamel 


A 199 Fi a pri 
Fie seria DA 13.20 și limnl——-—1 =p € R. Atunci, dacă: 


i) 


næl N—00 aat 


u>1> $a, (©) 
n=ł 


i) p< 1=. 5a, (D), 
n=l 


iii)u = ; caz de nedeterminare (nu se poate preciza natura seriei) 
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Exemple: 
o |] ua mr O A .. 
I ÈT Ty j cu a, =-7 20.: Am văzut cum criteriul rădăcinii ne 


n=l n n 
conduce la cazul de nedeterminare \=1. Vom aplica criteriul lui 


Raabe-Duhamel pentru a „scăpa“ din această situație. 


Avem: 
JI 2 lui 
u= lim n An a] iniae] lim n“ Mma n higit - | 
n= an+ n=00 n nyoo n 


= lim e: =2>1, deci seria converge. 
n=00 n 


09 ] ] . v . . v A . Lă 
a RU AY > 0. Se verifică imediat că atât criteriul 
n=i 2n+1 2n +1 


raportului cât şi cel al rădăcinii ne conduc la cazul 1=1. Dacă aplicăm 


2) 


criteriul lui Raabe — Duhamel obținem: 


= lim nfet- lim | 


a] 


n=oo (2n+1 n—c0 


p= lim | An —] 


F N Ani 
deci (în acest caz) nici criteriul lui Raabe — Duhamel nu ne „Scoate“ din 


cazul de nedeterminare. 


Dacă aplicăm în schimb criteriul de comparaţie la limită, considerând seria : 


oo % | a l 
ba =} —4(D) (seria armonică), obtinem: k = lım n = lim — =—, 
> g 2 o) jobi n-ob, n-o2n+l 2 


(D). 


t i (e9) 
Deci cele două serii au aceeaşi natură, cu alte cuvinte, seria Y` E 
n=i 2n + 
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Observaţii: 


i) 


ii) 


Criteriul lui Raabe-Duhamel nu se aplică (în mod obişnuit) decât atunci 


"când criteriul raportului sau rădăcinii ne conduc la cazul de 


nedeterminare a de În cele mai multe dintre situații, acest criteriu ne 
scoate din cazul de nedeterminare: în caz contrar se va încerca aplicarea 


criteriilor de comparație. 
Totuşi, există situaţii unde datorită formei termenului general al seriei 


aplicarea directă a criteriului lui Raabe-Duhamel aduce o simplificare a 


modului de calcul a limitei. 
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4.3 SERII ALTERNATE 


Definiția 4.3. 1.O serie Su, viy ER se numeşte serie alternată (cu 
n=l 


termeni alternaţi) dacă termenii săi satisfac relaţiile: 


Uau <0, (VneN T (4.3.1) 


Observaţii: 
i) Din (4.3.1) rezultă că termenii seriei alternează ca semn; atunci 


termenul general u, poate fi scris sub una din următoarele forme: 


u = Pta, a> 0 (dacă primul termen al'seriei este pozitiv) 
respectiv: 
| ue =(Cla, > a, >0 (dacă primul termen al seriei este negativ) 


ii) Deci orice serie alternată poate avea una din următoarele două forme: 


Ñ (— Da, = ayaz + a; — agtn n (= Dia 
z al a, (4.3.2) 
|Z CDa, = anta, —a3 +a.. H (CD an +... | 


so (n= 


Exemple: 
ce ] jaf pat 18 Log 
ea ) n? 22 32 4 =i n? 


O sA a D nt.. 
n=ł 


În studiul convergenței seriilor alternate nu pot fi folosite criteriile de la 
serii cu termeni pozitivi. Vom prezenta aici un singur criteriu specific seriilor cu 


termeni alternanți. 
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Criteriul lui Leibniz 


Fie seria alternată pe l)a, cua, >0. Dacă şirul (a, net satisface 


n=} 
condiţiile: 
i) a, —0 
ii)(a, )aen este monoton descrescător (a, >a, + (n eN” i 


0o 
. . n 4 
atunci seria 2 1) a, converge. 
nz 
Observații: 
=i) Criteriul lui Leibniz oferă numai condiții suficiente de convergență nu 


şi necesare; cu alte cuvinte dacă condițiile nu sunt- satisfăcute, nu 


rezultă că seria diverge, ci doar faptul că criteriul nu poate fi aplicat. 


~ di) “Totuşi, dacă lim a, = 0, atunci De a, (D), dar datorită aplicării 


n= 
criteriului general de divergență (im tee" a, = 0) | 
n—00 n—o0o $ 


iii) În cazul în care a, —0, dar şirul (ae nu este monoton 


descrescător, criteriul lui Leibniz nu este satisfăcut, şi nu putem afirma 


nimic dsăpre natura seriei SE (- la, 
n=1 


Exemple: 
IA I 


oo | | 
a) Aba A RA PER NR Urca ABE, > Arme, == ARIE 
= n în aii Bray 


alternată) : 


b) 


"deoarece: lim u, = lim(—1) 


Gi 


subcapitol. 
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A E 


| i 
Avem şirul a, =——0 şi este monoton descrescător, deci: conform 
n al r i, 


criteriului lui Leibniz seria converge. 


œ n n l 
— 1 — ; cu a, =—— = deci rima condiţie a criteriului 
x 3n +1 SRO FP ug R dia 3 


lui Leibniz nu este satisfăcută, acesta neputând fi aplicat. Totuşi, 


nu există, rezultă conform 


n= N—09 3n+ 


C.G.D. că seria atita. 


a sin? 22 sin 


] „Cu a, = 
vapori ~E cu a Fae 


Ea, 


Ji ! azi al te ale 
Deoarece 0<sin?n<! = 0sas <>. Trecând la limită (teorema 


Vn 


cleştelui) obti inem lim a, =0, deci prima condiție din criteriul ii Leibniz este 


n= 


satisfăcută. Se poate arăta relativ uşor. că şirul (an)nen, nu este monoton, deci 
condiția a doua nu este satisfăcută și criteriul lui Leibniz nu poate fi aplicat. 
Pentru a putea totuși preciza natura acestei serii, se pot aplica criterii de la 


serii cu termeni oarecare (neprezentate aici) sau folosi noțiunile din următorul 
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4.4 SERII ABSOLUT CONVERGENTE. 
SERII SEMI-CONVERGENTE 


Evident, în cazul seriilor cu termeni oarecare È I ARE R), nu se pot 
n=l 


aplica în studiul naturii acestora, criteriile formulate pentru seriile cu termeni 
pozitivi sau pentru seriile cu termeni alternanți, enunțate în secțiunile precedente . 


4.2, respectiv, 4.3. 


Totuşi, în secțiunea 4.1 au fost enunțate două rezultate care pot fi aplicate 
tuturor seriilor numerice, şi anume: Criteriul general de convergență a lui Cauchy 
(Teorema 4.1.3), respectiv, Criteriul: general de divergență (Teorema 4.1.5). 
Aplicarea lor, în special pentru stabilirea convergenţei seriilor, este aşa cum am mai 
spus, de cele mai multe ori extrem de dificilă. De aceea, matematicienii au căutat 
şi au stabilit noi criterii, mai simple de aplicat în practică, cum ar fi: Criteriul lui 
Abel, Criteriul lui Dirichlet etc. 

Chiar şi aşa, aplicarea lor necesită abilități matematice mai deosebite, 
precum şi 0 oarecare practică în lucrul cu serii Deoarece lucrarea de faţă 'se 
adresează în mod special studenţilor. economişti, s-a considerat utilă renunţarea 
prezentării acestor. noi criterii; vom arăta totuşi că, în anumite situaţii, studiul 
convergenței/divergenței seriilor cu termeni oarecare se reduce la studiul 


convergenţei/divergenţei seriilor cu termeni pozitivi. 


[es] 
Definiția 4.4.1. Fie seria cu termeni oarecare Y a ,, a, E R. Spunem că: 
l n=l 


(8,8) œQ 
a) Seria $a ,, este absolut convergentă (A.C.), dacă seria modulelor Y` la,| 
n=l n=l 
este convergentă (ca serie cu termeni pozitivi). 


F) Seria Y_a,, este semi-convergentă (S.C.), dacă 2 a, (C) dar, S lanl (D). 
n=l n=] 


n=1 
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had LA . œo i 
Teorema 4.4.1. Dacă seria Sa, a, ER este absolut convergentă, atunci 
n=] 


0o 
şi seria a, este convergentă. 


n=] 
Demonstrație: Deoarece, seria D aa (A.C.) = S aa} (C). Aplicând 
: n=l n=] 
Teorema 4.1.3. seriei modulelor, obținem: 
(V)e> 0, (n EN ai. langi cet au= apal + anu] <e, (v)n >n, 
şi VVkeN” ) Yu (*) 


Din proprietățile funcției modul, rezultă că: 


aaar teh aanl Sapt... anu] (**) 


i Oo | j 
Din (*) + (**) rezultă că Da, (C), conform criteriului general de 
i RM n=l AR 


convergenţă-a lui Cauchy. 


q.e.d. 


Observaţii: 
1. Reciproca teoremei de mai sus nu este adevărată; cu alte cuvinte, 


convergența seriei simple nu implică convergența seriei modulelor 


a, (C) =æ la (O 


IE Au gs7 irpan wih 2 | ` 
Exemplu: Fie seria 5: aa = D (— 1)” S Evident este o serie alternată 
n=] n=l 


şi converge deoarece satisface criteriul lui Leibniz. Totuşi, seria 


modulelor, 5Y lanl = $ —, diverge ca serie armonică simplă. 
n=! n=I N 
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2 . . g9 . . 
2. Intrucât seria valorilor absolute, Sola, |» este o serie cu termeni 


n=] 
pozitivi, pentru studiul naturii acesteia se vor folosi criteriile enunțate 


în secțiunea 4.3 


(2.9) 
În cazul în care seria modulelor Y lan] (D), am văzut că despre seria 
n=l 


IE. | 
simplă 3 a, , nu se poate afirma nimic. Totuşi, se demonstrează că în 
n=! P A ; i 


: (3.9) 
cazul în care divergența seriei modulelor p la, a fost stabilită cu 
z n=! 3 


Criteriul raportului sau cu Criteriul rădăcinii, atunci şi seria (simplă) 


Me 


n diverge. 


3 
i 


Pentru seriile cu termeni pozitivi È Bo a 0) noțiunile de 
n=l j 
convergență şi “convergență absolută coincid (deoarece 
=: ad 
Da, = 5 la,|). 
n=l 
In anumite situații, este utilă şi următoarea formulare echivalentă a 


X | 
Teoremei 4.4.1. „Dacă seria Y a, a, ER diverge, atunci şi seria 


nT" ; 


8 


Pi i s ce 
aa] diverge. 


n 


A [2,94 
In concluzie, în studiul naturii seriilor cu termeni oarecare, Sa, a, ER, 


ihl 


se poate aplica următorul: 
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Mod de lucru: 


| oo 
i) Se ataşează seriei cu termeni oarecare, seria modulelor Sa, (care 
n=l 
este serie cu termeni pozitivi); 
11) : Se studiază natura seriei Y la, folosind criteriile enunțate pentru serii 
n=l 
cu termeni pozitivi; 
iii) Dacă: 
o OD 
a) D akO S E anO) e 
n=l n=] 
N 9 .. .. . b [] Å . Lă . 
b). D lanl (D) (divergența fiind stabilită cu Criteriul raportului sau 
n=l 
. . . .. g9 
Criteriul rădăcinii > J a, (D) 
n=l 
f g | | RI . . p .. 
c) D fanl (D) (divergența fiind stabilită cu alte criterii) = despre 
‘n=l | - ! i 
. LI 09 . . 9 | s . 
natura seriei X-a, nu se poate spune nimic. În acest caz trebuie 
: n=] 
folosite alte criterii (Abel, Dirichlet etc.). 
Exemple: ; i 
EI 2" sin” a 
D De, cu aE [0,27] parametru real. 
f=l Vn 


Evident seria este o serie cu termeni oarecare (datorită lui sin” a ) şi nu 


cu termeni alternanți. 


„e œ 2" [sina 
Ataşăm seria 5 a,]= Di mit şi aplicăm criteriul raportului. 
n=l n=] n o- 


Atunci, obținem: 
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ai local À ' 
pee [Pr sin al ali | 
= — rm ————— = lim 2 SAS al 
n=% a, n= Jn da] | n-o ngA 


şi deci: 


XE 


| A= 2|sin a] 


Atunci, pentru: 


i) AZ 2|sin a] <1, seria Sla, COS Sa a, (C) (conf.T.4.1.1) 
n=] i= À 


li) A = 2|sin a|>1, seria S lanl (D) > SS a, (D) (conf.Obs.3) 
n=] l 


odii) N= 2|sin a|= lnu putem preciza natura seriei > ja, | şi deci a 
i n=l 


o 
seriei inițiale Sa, . 
n=] 


Rezolvând inecuațiile trigonometrice de mai Sus, avem: 
n) [2x7n) (Lin 2" sin" a 
o] 22) [ie za] 


i) pentru a € > pat 1) FF (A.C.) 
n 


ii) pentru Lee jfa i)a $ Cea Mb 2" -sin' a (D) 


6'6 pre “n 
y 2n 7m lln tă 3 £ 
iii) pentu a«€ Z | => nu putem preciza natura seriilor 
6? 3 66 
> Rn | şi pată 


Dar, în cazul iii) cele două serii (simplă şi a modulelor) au fXpipiie: 


| 
a) 3 a” =>, (—1I re: (0), ca serie alternantă care satisface criteriul 
n=l n 


ArT- e l Fi 
lui Leibniz (şirul —= — 0 , monoton descrescător). 
n 


2) 
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b) jaj = EEE 73 (D), (ca, serie armonică yA- cu 
n=l n=I VN n=. 7 n=un' 
| 
a=—<1) 
2 


a bachi) X PV ae EE 
Deci, în acest caz, seria Y (+1) S ala este (S.C.) 
n=| 


sin” 
n 
Y >] | l-a Ha , cu a er\f5] parametrul real. 


| — 2a 


| | k se o6 I-a f ; 
Ataşăm seria modulelor: Y i | se N E | La şi aplicăm 
, "e ns] n=] n l— 2a 


criteriul Acin; 


n+l L>a n+l |l—a l—a 
= = lim n la, = = lim a aie E a 
[E za B- i n i pica NE 


n= 


Atunci, pentru: 


i) >=] iz d < K D aal (C), deci şi seria Ea (C) (conf. 
l— 2a n=l n=1 
T.4.4.1) 
AS PE TA iio > las | (D), deci şi seria Sa (D) PEN Obs.3) 
WAS ; 7 = 1 => nu putem preciza natura seriilor Shan | SEDD ə i ll 
e na] 


Atunci, rezolvând inecuațiile de mai sus, concluzionăm că: 


i) pentru X s(coso)u[ ze] = D =] h F | este (A.C.); 
a 


n=| n 


ii) pentru vefosh [1] SzH] | e | este (D); 
3 n=l n l— 2a 
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iii) pentru XE (EI => nu putem preciza natura seriei. 


Dar, pentru: 


+1)" 
„care (D) conform C.G.D. 
N I 


Intr-adevăr, 


- i ME 0 "ii 
lim a, = lim ka = lim ES =ez0 
n= n>œ| n 


N—o n 


n 


care (D) conform C.G.D., deoarece a, 40. 
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45, SERII DE PUTERI. DEZVOLTAREA 
„FUN CȚIILOR ÎN SERII DE PUTERI 


Aşa. cum. s-a observat din secțiunile precedente, noţiunea „de „serie. 
numerică“ apare ca o generalizare (extindere) firească a noțiunii de. „Şir, 


numer ic“, i 

În mod natural, aceste două noţiuni au căpătat în Analiza matematică noi 
extinderi, din ce în ce mai complexe, lucru de altfel deseori întâlnit nu numai în 
matematică ci şi în multe alte domenii ale științei. 

Astfel, pornind de la noțiunea de şir numeric: 

(2, ser: TE E: EE du + (*) 
se HORACE noțiunea de „şir de funcţii“: 


Ci :fo(x), ie S SI PIE (ART, ga "E (e) 


unde funcţiile f (x) sunt definite pe acelaşi domeniu DC R. 


Exemple: 


f:R—R 
f (3) = ix” sin nx 


a) (fa(x)) 


nene Cu termenul general: le 


adică, şirul de funcții are forma: 
. r e G e 2n + 
xsinx, 2x sin 2x, 3X 'sin3x,...,nx"sinnx,... 


ră APE cati 


b) (09) en „ unde termenul VERAS: este: H (k) = (n2.+ DIn(nx) 


deci, şirul de funcţii are forma: 


2In x, Sin2x, 101n3x,...,(n? +1)In(nx)... 
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Ținând seama că noțiunea de „funcţie“ este mult mai complexă decât cea 
de „număr“, este evident că noțiuni precum convergență sau divergență au căpătat 
pentru şirurile de funcţii valenţe noi. Astfel, un şir de funcții ERA poate fi: 


simplu (punctual) convergent, uniform convergent, slab convergent, tare 


convergent etc. 

Este natural ca în caz de convergență, şirul de funcții să aibă limita o funcție 
f(x) şi nu un număr. De asemenea, ținând cont că termenii şirului sunt iea care 
pot avea diverse proprietăți (funcții continue, derivabile, integrabile elc.), este natural 


să cercetăm în ce condiţii aceste proprietăți se „transfer ă“ la f uncţia limită f(x). 


In mod cu totul analog noţiunea de „serie iumericăt: 


e 


An = 40 +a, +a, spusa | S E) 
n=0 : : 


capătă o primă generalizare şi anume cea de „serie de funcţii“: 
= 
2 fa(x)= fo(x)+ E (X)+ (x) ++ felhs PEN 
n=0 
unde termenii seriei sunt funcții definite pe acelaşi domeniu DC R. 
Exemple: 


a) 5 (0 %4. Darctg(x" +1) = arctgl + 2arctg(x + 1) + Sarctg(x? + 1)+...un 
n=0 d: d 


:R—R 


de termenul general al seriei este: ih (x) = (n2 + Darcte(x? +1) 


es na 1! 2! Sf. 
b) 5 == + ze tei fe A 

ni 1+X" 14x L+x2 i+ 
loa 


unde termenul general are forma: 


E = 
= spa 
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Analog, în caz de convergență (simplă, uniformă etc.) seria va converge la o 
„funcție sumă“ s(x). Din nou este interesant şi important de studiat, în ce condiții 
proprietățile termenilor seriei de funcții (integrabilitate, derivabilitate, continuitate 


etc.) se transferă asupra funcţiei sumă s(x). 


4.5.1. SERII DE PUTERI 


Un caz particular de serii de funcţii, dar foarte important şi des întâlnit în 
aplicaţiile economice, este cel al seriilor de puteri, unde termenii seriei sunt funcţii 


putere. Astfel, avem:. 


Definiţia 4.5.1: Numim serie de puteri centrată în Xo» O serie de funcţii de 


forma: 


00 1 i $ 2 

2- an (X-X0) ao +4; (X-X0) Fas (x-x9) + e Han (X-X... (4.5.1) 
n=0 

unde (a, re este un şir numeric (numerele a, ER se numesc coeficienții seriei 


de puteri). 


„Observaţii: 
i) Un caz particular, dar foarte des întâlnit în aplicaţii, este cel al seriilor 
de puteri centrate în 0 (pe scurt, serii de puteri) deci când Xp = 0. Vom 


obține atunci seria: . 

„99 5 

Dax Sag Fax Haak + ra ku (4.5.2) 
n=0 


ii) Deoarece prin substituţia: 


se i-a pa | (4.5.3) 
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o serie de puteri centrată în x (Şi variabila x) se reduce întotdeauna la 


o serie de puteri centrată în 0 (şi variabila y), adică avem: 


to Dax E Say" (4.5.4) 
13 n=30 n=0 


ne vom ocupa numai în continuare numai de seriile de puteri centrate în 
0, deci de forma (4.5.2). 


(9.9) 
Definiţia 4.5.2: Spunem că seria de puteri 5 a x" este: 
neg 


a) convergentă (punctual) în x =X €R, dacă seria numerică ataşată, 


[9,9] 

=f} = K: 
2 â, Xo este convergenta; 
j= 


b) convergentă (punctual) pe domeniul Do CR, dacă converge în orice 


punct x € Do (Do se numeşte mulțime de convergenţă a seriei). 


Teorema 4.5.1. Mulțimea de convergență a unei serii de puteri, 
(o s pe | 
2. a„,X este întotdeauna nevidă. 
n=0 

. . . . 9 n . 
Demonstrație. Fie seria de puteri 2. anX cu (a, en UN şir de numere 
î= 
reale oarecare şi Do mulțimea sa de convergență. 
~ . e n . . 
Pentru x=0 seria 2 nX Æo, deci converge şi are suma S=ay. 
D= 


Rezultă că 0 € D,, deci Doz6. 


q.e.d. 
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Teorema 4.5.2. (teorema lui Abel) 


co E 
Fie seria de puteri, 5 a,x", a €R. Atunci, (E)r e [0,00] Gr“ se 


n=0 
numeşte rază de convergență), astfel încât: 
. . >g 
i) pentru x € (~r, r) = Ù a,x” converge (C) 
n=0 K | (4.5.5) 
ii) pentru x € (=00,=r)U (r, +00) => Ð a,x” diverge (D) 
n=0 

Observații: 
i) Teorema lui Abel nu precizează natura seriei în capetele intervalului de 

convergență: (-r, r). În exemplele concrete, se studiază separat natura 


seriei în cele două capete, adică seriile numerice: 


18 


n=0 


` (e) 

Dima n 

ar Şi Sa,(-r). 
n=0 
.. . . . . Dy n 
ii) Dacă raza de convergenţă r = 0 , atunci evident seria de puteri 3. aX 
n=0 
converge numai în x= 0 şi diverge în rest (pentru x € R“). 


iii) Dacă raza de convergență r = +00, atunci mulțimea de convergență a 


seriei Dg = R, deci seria nu diverge pentru nici o valoare reală. 


iv) Aşadar, folosind teorema lui Abel, dacă aflăm raza de convergenţă r, 
. | . . es n i . . 

ştim mulțimile pe care seria O a,x converge sau diverge (evident cu 
| d, r=0 i 

excepția valorilor x= tr). Următoarele două teoreme ne oferă două 


formule pentru aflarea razei de convergenţă r. 
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Teorema 4.5.3 (Teorema raportului) 


À a aapea 

° . . e e +1 
Fie seria de puteri, 5 a,x” cu a, ER şi lim = 

n Pia fanl 


raza de convergență r , este dată de relațiile: 


l 
~; 3P E (0,00) 
P 

r=4{+00 ;p=0 
O'T Eps ROO 


Teorema 4.5.4 (Teorema rădăcinii) 


=DE [0,0] . Atunci 


(4.5.6) 


L.] E. L Cad e . 
Fie seria de puteri, ý a x^ cu a, ER, şi lim ylanl =pE€ [0,00] ş 
n=0 -a 


Atunci, raza de convergență r este dată de relațiile: 
TE 
— PELOR) 
P 


r=+oo ;p=0 
0 ahoo 


(4.5.7) 


Observație: Din rezultatele enunțate anterior, se observă că în determinarea 


mulțimii (intervalului) de convergență, respectiv, a mulțimii (reuniune de intervale) 


de divergență pentru seria de puteri 32 a,x”, trebuie aplicat următorul: 


n=0 
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“Algoritm de lucru: 


1) Se calculează valoarea „p “ aplicând una din formulele: 


o= lim fans sau p = lim niga] 


n= a 


2). SĂ determină raza de convergență „r* cu relaţiile (4.5.6) sau (4.5.7) 


(care de fapt coincid); 


x 0o 
3) Se scriu intervalele de convergență / divergență ale seriei X nE: 
n=0 


4) Se studiază natura seriei în cele două capete x = +r.. 


Observație: În cazul în care seria este centrată în a a(x -xo se 
: n=0 


reduce la seria centrată în 0 È ay") folosind substituția (4.5.3) şi apoi se aplică 


n=0 
algoritm de lucru, obținându-se mulțimea de (C)/(D) pentru variabila y. Folosind 


aceeaşi substituție se scriu mulțimile de (C)/(D) în variabila x. 


Exemple: 


ca n 
DT LIRĂ; Avem: a, aiet. deci a [aa 
ha] n n n 


Folosind teorema rădăcinii, obținem: 


P= Abe t fs aid i 


(am folosit aici limita fundamentală: lim Y/n =1) 


n= 


2) 
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deci conform (4.5.7), avem: risi deci r=1. Atunci, conform 
p 


Teoremei 4.5.2. (Teorema lui Abel) obținem: 


P 7 e ae i FRE (o 
n=i n 
ii) pentru x E€ (—œ,—1)U (l,o) > > (—1)" A (D) 
n=l n 


x’ 
iii) pentru x = £1 = S (— 1)" — nu-i putem preciza natura. 
n=] 


Dar, pentru: 


x=1l obținem: Di pică =O =- a (conform Criteriului lui 


i 
Leibniz de la serii alternante) 


iar, pentru: 


x = —1, obținem: Sax: = SL (D) (ca serie armonică S DEN cua=l) 


n=l n=1 N n=l n- 


(C), dacă x €(—1,1] 
În concluzie, S- pe i 


n=l (D), dacă x € (—020,~=1]JU(1, +00) 
n+! X n+l l : l 
E 1) X avem a, =(D"—, deci la,|=— 
! n! n! 
Folosind teorema raportului, obținem: 
an n! 


= mM——— LA iezi 
n—00 lanl n—=oo ai n—con+-l 


ele A 
şi deci conform (4.6.6) avem r = +oo. Deci, seria Y (—1)*! 2 (C) 
n=0 N. 


pentru (Y)xER. 
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OD 
3) $ E (x +1)", este o serie de puteri centrată în Xp =-—l, 
n=0 
1+1 
Notăm y= N "cua, iz | (| 
ta ; 

Atunci: 

a i 2 +2 | 
sat: del pati bus si n+2 _n2 ul RASI +a | 


1 —ma 22 
nas Sanji Amn ni i 2an 2 
Atunci raza de convergență r = = = 2, Şi avem: 
P 
i) pentru y € (—2,2) = 5 niays (C) 
n=l n 


„ii) pentru y € (—0, —2uG, 99) n Zi 


iii) pentru y = +2 = ?(nu putem preciza natura seriei) 


Revenind la variabila x(x = y—1), avem: 


i) pentru x € (— -3,D => aR: L x41) © 


ii) pentru, x € (—00,—3)U(1,00)= Èr 


L (Dn (D) 


iii) pentru x = —3, x = 1 = ? (nu putem preciza natura seriei) 
Dar, pentru: 


X = —3, obţinem D a, (x+1)" = = n= (D) (conform C.G.D, 
—— n=] 


deoarece termenul general (— 1)" pri 
| n 


+0) 
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iar pentru: 
i oo no & nl 
x. = L obținem, 2- a Wah" = (D) (conform C.G.D, 
t n=l - n=l . | 
deoarece termenul general sur — 10) 


4.5.2. DEZVOLTĂRI ALE FUNCŢIILOR ÎN SERII DE PUTERI 


În secțiunea precedentă, am arătat că unei serii de funcţii (în general) 
respectiv unei serii de puteri (în particular), i se ataşează (în caz de convergență) o 


funcţie s(x) numită suma seriei. 


Vom arăta în cele ce urmează, în ce condiţii are loc reciproca, adică: când o... 


funcție oarecare f(x) poate fi considerată suma. unei serii de puteri 


(£()= Saxo) sau f(x)= $a, x?) 
n=0 n=0 


Definiţia 4.5.3: Fie f IC R-—R,I7 interval, şi Xp € I. Spunem că: 
a) fe C" (xy) (este de clasă C” în xy) dacă: | 
(3) (x6), (V) ne N (adică admite derivate de tă ordin în xo). 
a) fe C" (1) (este de clasă C” pe intervalul I) dacă: 
(3) f(x), WneN şi (V)xpe I (adică admite derivate de orice 


ordin în toate punctele intervalului 1). 


Definiţia 4.5.4: Fie f:ICR-R, xpel,feC"(xg). Numim serie 


Taylor ataşată funcţiei f în punctul x, seria de puteri: 
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co f ™ f 
tii Xo)” = f(x9)+ Boa =) 
pi f A (4.5.8) 
f 
+ a =x) +. ma S ax 0) +a 
Observaţii: 
i) În relaţia de mai sus am utilizat convențiile: 
AA 
f(x) = f(x) 
ii) dacă Ko = 0, seria lui Taylor devine: 
co 'fin) e fa e 
TEO n=O Zi EE Dan. (4.5.9) 
n=0 n! Pe n! 


Şi se mai numeşte seria MacLaurin ataşată funcţiei f(x); 


iii) Evident seria Taylor ataşată lui f(x) este o serie de puteri centrată în x, de 


forma Xa (x -Xoj Ude a, = a „iar seria MacLaurin ataşată lui 
Fe 
f(x), este o serie de puteri centrată în 0, de forma Sank" unde a, = mo j 
n= ! 
“Exemplu: 
| TAIR, ; Evident f € CR) şi fR) =e  (neN 
| Atunci scria Taylor ataşată funcției exponențiale în x = 1, este; 
D Koe" TER pes Spt. ien $! (4.5.10) 
respectiv, seria MacLaurin ataşată este: 
$ AREE r e, (4.5.11) 


n=0 N! N A n! 
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Evident, ne putem întreba: „în ce condiţii suma seriei de mai sus este 


tocmai f(x) = e“ şi pentru ce valori ale lui x are loc acest lucru ?“ 


Teorema 4.5.5. Fie I un interval simetric față de x, şi fe C”(1). Dacă 


(3)M, > 0 astfel încât: 
[OCS Mo, (Vixe şi (Wne N” (4.5.12) 


atunci are loc relația: 


PUN S DEEN 


le FN 


TTA pA 
mae 


(x= x w) i. dl 
(4.5.13) 


“Tree (cuget. = Apă (x — xo)", (V)x el 


numită dezvoltarea în serie Taylor : a funcţiei f(x) pe intervalul I. 


Observaţii: 
i) Evident luând x9=0, în teorema de mai sus, se obţin condiții 


(suficiente nu şi necesare) pentru a avea dezvoltarea unei fi uncţii în serie 


Maclaurin, adică relaţia: 


p(n) f f” f(n) 
pi d MANE E e apr Ra (4.5.14) 
! ! n! 


ii) De obicei, pentru o funcţie oarecare f(x) se caută dezvoltarea în serie 


Maclaurin (are o expresie mai simplă), dar pentru aceasta trebuie ca 
fe C”(0); în cazul când o funcție nu este definită în x =0, se 


determină dezvoltarea în serie Taylor pentru x= Xg Oarecare. 
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Exemplu: 
Considerând din nou f(x) = e* şi Xp =0 avem: 
„food =e* = f9(0)=e%=1, T 


Atunci, evident condiţia (4.5.12) este satisfăcută (V)x e R , deci putem scrie că: 


3 


n ` yn 
l eea pei CR (4.5.15) 
n! S atar n! ki 


Š © 
e = 
n=0 


Considerând în (4.5.15), x=1, obținem relația cunoscută de la liceu: 
ON ea ERS era me | (4.5.16) 
H 2r J OONA i 
- Observație: Evident pentru a determina dezvoltarea în serie (Taylor sau 
MacLaurin) a unei funcții f(x) oarecare, trebuie mai întâi ca aceasta să fie indefinit 
derivabilă (fe C”) pe o mulțime şi apoi să calculăm derivatele sale P 
într-un punct Xg; se înlocuiesc aceste derivate în formulele (4.5.13) sau (4.5.14) şi 


se determină din condiția (4.5.12) mulțimea de convergenţă. 


Vom da mai jos, câteva dezvoltări ale unor funcţii elementare, precum şi 
mulțimile de convergenţă ale acestora, lăsând la latitudinea cititorului verificarea 


lor. Astfel, avem: 


XR Sep La A aura „a E 
a) f(x)=sinx = —— TR Fu, 4.5.17 
) 109 2 Cn FII Ha ră, Poan 
(V)xER 
b) f(x)= (+x) = 5 ED pe = (4.5.18) 
n=0 n: 
EENE a e + P Caz: CGE. ad: n 


1! 2! d, n! 
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antr (V)x € G 1, 00) 
pentru [E ERIN 
care poartă numele de seria binominală. 


Dând diferite valori lui a se obțin din (4.5.18) diferite dezvoltări de funcții 


, obţinem din (4.5.18) 


5 | 


în serii de puteri. De exemplu, pentru a = 


— E Da da 
VEX =la- pă 2 BN n, a 


2! 4.2! 8-3 2".n! 
> îi fa i (4.5.19) 
pentru (V) x e(—1,00) 
Considerând în (4.5.19), x=1, obţinem formula: 
Mircia LI OE 2 e i a CERU cr AIA (4.5.20) 


71-04 e pain, 


| via parte aa ne a 
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4.6 PROBLEME REZOLVATE. PROBLEME PROPUSE 


In această secțiune vom prezenta alte câteva exemple de serii numerice şi de- 
serii de puteri rezolvate şi vom propune spre rezolvare o serie de exerciţii. Testele 
grilă prezentate sunt similare cu cele propuse pentru examen şi au scopul de a fixa 


cunoştinţelor teoretice prezentate în acest capitol. 
A) Probleme rezolvate 


A1) Determinați natura următoarelor serii cu termeni pozitivi: 


co Ah n2 n, | 2 n+l. 12 
n Sa (n! y Aver $, E. (n!) Da a 4 (n+1)! 
n=1 (2n)! (2n)! | (2n +2)! 


Atunci, conform criteriului lui D'Alambert, avem: 


nl 2 
4 ((n+ 1!) (ny j -yi 4:0 +)? 


A= lim PE! = lim S EA 
io a, n= (2n+ 2)! - “hai -(n n? n=œ (2n Fn) 
„ 2n+2 l 
= lim = 
‚n=, 2n, 4] 


deci nu putem preciza natura seriei. 


Vom aplica în continuare criteriul lui Raabe-Duhamel, deci calculăm: 


. 2n +1 . —n l 
-l|= limn —]|= lim = 
n=co |(2n+2 = 2n + 2 2 


p lim | An 


ET an+ 


= 4 (n? 


| 
Deoarece u=-—<]1= D). 
y 2 = (2n)! (D) 


1295, 23 cu azl şi pER parametri reali. 
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e 


Deoarece a21, termenul general Tap 20, deci seria este cu termeni 


pozitivi. 


Aplicăm criteriul lui D'Alembert şi obținem: 


_pn+l  -2—3p 2 
A= „lit ui cd IE lim Îi ii i = lim an|- Sai “=al 

n—0o Aa | are IP 3p (a — Dr n= l 

Avem: 

i i i . LE s : ; co (a =1) iy i 

i) pentru X=a-l<l (a€[l,2)= 5 i 40), 

NR] 
sA 
ii) pentru \=a-i>] (a e (3, to) 5 Sa mop (D) 


iii) pentru X=a-l=1l (a=2)5 nu putem preciza natură seriei. 
x i Ra a ali! | p 
Dar, în acest ultim caz, seria Y Ars) —323p Care este seria armonică 
Yi n=] =i N i 
Ba i a a. 
generalizată 5 — cu «=2-3p. 
n=1 N 
Atunci, rezultă că: 


a) pentru «=2-3p>l (adică pe - => seria Sa, (C); 
| ‘a pl 


b) pentru «=2-3p<1 (adică pe [1/3,%)) => seria Sa, (D). 
n=l 
În concluzie, seria: 


L TE d 
(C) dacă: (aclk2);nea)ula= upe| a) 


(D) dacă: (a € (2,+00),p€ RU (a = 2,p€ [173,-+00)) 
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2 n 
a n +n+l 
3) È| a:———— |; a20 parametru real. 
n= n +n+2 


Aplicăm criteriul lui Cauchy, pentru termenul general 


$ 2 
a, = A tea >0 
n +n+2 

Avem: 


n A 2 

i N +n+l] 
= lim e rem. 
n-o n'+n+2 


X= lim ya, = lim 


n= n= 


l | n2+n A2 
Discuţie: 


1). pentru A=a<l (deci ae [0,1)) > 3 apa) (C); 
as în +n+2 


| d i, 1 | 
ii) pentru À =a > I (deci a TS T baeu (D); 


n=l n2+n+2 


iii) pentru À =a =1, nu putem preciza natura seriei. 


n 
r Ema la AET]. s 
Dar, în cazul a = 1, seria Y dE | - In acest caz, nu vom 
n=l asi h4+n4+2 


încerca să aplicăm Criteriul lui Raabe- Duhamel, deoarece limita rezultată este 


dificil de calculat. În schimb, observăm că: 


nae  noe(n?+n+2 3 


n 
-(n2+n+2) DAT 
n? +n Le, l 
lim a, = lim IA AR pa n-a) i =e0% =] 
n +n+2 


Deci a, — 1 (#0) şi conform C.G.D., seria diverge în acest caz. 
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In concluzie, seria: 


(C) ; dacă a€ [0.1) 


: n 
Pa n+n+l|. 
n=l n“+n+2 
(D) ; dacă a€ [1,+00) 
œ 1 ppa a Mr 
A Dn ea T 
n=l] n 4 


? l n 
PEPA HS n 
Despre termenul general a, = Ea Fa Te y nu se poate afirma 


n 


imediat că satisface relația: a, 20. 


De aceea, vom aduce pe a, la o formă mai simplă, utilizând identitatea: 


2 2 

PP 4 y n =A RED ew 

Atunci: 

SpA mn +a] _ n antl m) 2n +1 AN 
i 4n? „4 4n „An 4n i 


Folosind, criteriul rădăcinii obținem: 


à mna ii A a at ae AA 
X S& limi a, = lim p = lim Di i 
n=00 n—co 4n n>œ 4n 2 


Prin? n) © 
n? 4 i 


Deoarece N = 


l le) 
-<l> 
2 2 
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A2) Determinați mulțimea (domeniul) de convergență (punctuală) a 


următoarelor serii de puteri: 


i n? 3 


1) E pie 


an 


iezi deci la, 


E 
n“ 4" 


TARET 


Avem: a, =(—1) 


Notăm (*) y = x —3 şi seria > a (x-3) = D ay". Calculăm: 


uh. 


=i lan] a a a i Ta ETT £ n+l (22) 
p= lim == = Du lim Ia 
n: [an st Na m n+2 


„Atunci, 


, Lari] 
Deoarece p=4= raza de convergență r=-—= 4 
p 


Tepremei lui Abel, avem: 


i) eniru e| 2.2 Say (C); 
4 4 Mal 


/ œ% 
ii) pentru ve[- au [++ %)= TDA (D); 


... l 4 . . A | n 
111) pentru y= gi nu putem preciza natura seriei $ a y”. 
n=l 


Revenind la variabila x, avem datorită notației (*), următoarele cazuri: 


i) D xeļž, z> S Sa, (3%. (O); 
4'4 
11 13 vi fi 
ii) pentru xe[-ax jurat) > IAA 0); 
n=l 


conform 


11 13 > ap $ 
iii) pentru x € | =] » nu putem preciza natura seriei J` a, (x —3)" 


SF 


nl 


162 MATEMATICI APLICATE ÎN ECONOMIE 


Dar, pentru: 


2 


1] o GERAN "HER 
a)x=— > Da (x-3"= 5 (1° (D)(cf. C.G.D, 
) 4 > ) > ) Nn+l | 


deoarece a, 4 0) 


(D) (cf. C.G.D, 


deoarece a, 40) 


2) Sp np De t/n 
nasi 


ant! AEAN oep 
4 | 
Avem a, = (opri. Ani deci la, Sola Atunci, ana: 


(*) y =x +1; obținem: s: a"(x PD" D a"y" . Calculăm: 
n=l NE. 


W n+l g 
dhupi eal EM netDwntI n+l dă e ca a. CEREN 


n= la, | n=09 qnt2 i ntn n=o0 aani 


Atunci raza de convergenţă r = 0, deci seria: 


i) Say (C), pentru y=0 
n=l 


ii) Sany’ (D), pentru. yE (—20,0)U(0,œ)= R* 
pia i 


Revenind la variabila x, rezultă din (*) că: . 


(C), pentru x= —1 


nyi Nt! t/n 


qn+i 


ota 1) 


(x +1) 


(D) , pentru x € R\{-1} 


paria 9" t Diy pai 


) 2 
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1)-9" 
Avem a, Pi aa mad] Ș 
n 


şi notăm y =x + 2 (*). 


i IO +DR o g | 
Calculăm p= lim aja, | = limp eee lim =-Yn +1 =041=0 
n= n—co n’ n=% n 
Deoarece p=0, raza de convergență r =+ . Atunei, conform Teoremei 
0o 


lui Abel, rezultă că Say” (C) (V)yERA şi deci -folosind 


n=l 


(53 Daat" (0, V) xER. 
pæl 


B) PROBLEME PROPUSE 


B,) Studiaţi natura următoarelor serii cu termeni pozitivi: 


o n+l 

a) nes R: converge 
= n +i | sm A 
00 DA 

b L -= R: diverge 
n=] Jn(n +1) p 
00 n n - 

ey ay A it R: diverge 
n=} n À 
oo ty 

d) — m R: converge 
> 2" 45" i 
; 3 n 
% l 

el 2 a = | R: diverge 
ni (N n 


R: converge 
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B2) Studiaţi natura seriilor alternante: 


œ% n i 

a) 5 (-1) R: diverge 
n=l] | Nn+l i 
oo n 

b) S(—I R: converge 
P '(n+1)! Š 
eu n+I| 

Cc): DASD =] R: diverge 
næl și 
60 n+l 

d) Feen R: converge 
n=l] (n—1)! 


B3) Determinaţi mulțimea de convergență (punctuală) a seriilor de puteri: 


a) 25 R: x €[-1,1] 
n=l I 
b) DE > > R: xe(—1,1] 
n= n“ 
œ Zi r3 
c L. ewl 2al)’ R: rejk = 
i > ; m, = ah 202 
d) MOL ar AGA R: x €(—10,8) 
næ] il 
[2.9] 12 i i 
e) SOR R: x € [0,8] 


n=1(2n)! 


Capitolul 5 


FUNCȚII REALE DE N VARIABILE REALE 
A i ii a ici zii 


Fenomenele vieţii reale sunt descrise cel mai adesea de funcţii care depind 
nu de o singură variabilă, ci de mai multe variabile indepeñdeiite. Acest lucru este 


valabil şi în viața economică. 


În acest capitol se introduc elementele de bază ale teoriei acestui tip de 
funcţii. Pentru a depăşi dificultățile extinderii conceptelor fundamentale 
(continuitate, limită, derivabilitate etc.) de la funcții de o variabilă la o funcţie de 
mai. multe variabile, am abordat numai cazul funcţiilor de două: variabile, 


dezvoltarea ulterioară făcându-se în aceeaşi manieră. 


5.1. ȘIRURI DE PUNCTE ÎN R” 


Să considerăm spaţiul vectorial R” şi să considerăm un vector oarecare 

X = (Xp Xa-:+5Xa) al său. Pentru a simplifica scrierea vom renunța la simbolul de 

transpunere şi în cele ce urmează vom scrie pur şi simplu x = (Xp X2. Xp). 

Acestui element al spaţiului R” îi asociem un punct P care are coordonatele 
4 š Dai i 

X1sX2s--:3Xņ » tot aşa cum unui vector x = (x;,x,) din R? îi asociem un punct P 


cu coordonatele x,, x, într-un sistem de axe ortogonale. 
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Asociem unei perechi de astfel de puncte P, Po  P(XpX25- Xp) 


P,(x0,x0,...,x9) un număr nenegativ numit distanță şi definit prin: 


d(P,P,) = (xix + (xx (ax. CATY, 


Evident, distanţa este nulă dacă şi numai dacă P şi Po coincid. Convenim că 
însuşi vectorul x să se numească punct din R” . În acest mod, distanța d se asociază 


unei perechi de vectori din R”. 

Definiţia 5.1.1. Se numeşte şir de puncte din R" funcția care asociază 
oricărui număr natural k un punct unic x din R”. 

k > Xk = (Xk Xk20000Xkn) 


„Şirurile (Xpy kens (Xk2 ke» «ie (Xkn)ken Se numesc şiruri de coordonate. 


„Definiţia 5.1.2. Sirul (x, en converge la punctul g dacă: 


lim d(x,.X9)=0 (5.1.2) 


În spiritul cunoscut al limitei unui şir de numere reale aceasta înseamnă: 
(Y) e > 0,3 .k(e) aî.(V)k > k(e) d(x, xp) <E AS: 1.32 
Folosind (5.1.1) deducem că: | 
Vika = X01)? + (Xa — X02)? +... + (Xa z Xom) <E 
pentru k > k(e) adică | 
sua = xal e kakek CD i=ln (5.1.4) 
suntem conduşi astfel la: 


Teorema 5.1.1. Şirul x* este convergent la x’ dacă şirurile de 


coordonate x, sunt convergente la x; pentru orice i =1,n. 
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bs iot R 2k f~ PEs 
Exemplu: Considerăm în RR” şirul Xk | ză Şirurile “de 
2k +1 Erp a. ; 
coordonate ‘sunt -x,, -5 X3 =i ŞI pentru k—oo au “limitele 0 şi 


nlan i T ati rul 
respectiv rý Aceasta înseamnă că x, are limita xp = $a i 
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5.2. FUNCŢII DE DOUĂ VARIABILE 


w . . 2 . 
În cele ce urmează vom considera spaţiul R? ale cărui puncte x au doar- 


două coordonate x = (X,,X2). Pentru ca scrierea să fie mai simplă vom înlocui 


perechea de coordonate (X4, X, ) prin perechea (x,y). 
Fie D o mulţime de puncte (x,y) din R? şi lo mulţime din R. 


Definiţia 5.2.1. Se numește funcție de două variabile tripleta formată din 
D, I şi o lege de corespondenţă care asociază oricărui punct (x,y) € D, un singur 


element ZEI. 
Exprimăm această corespondenţă prin z= f(x,y). 
Observaţie: Când D este din R” iar ICR în mod similar se defineşte o 


funcţie de n variabile. 


Dacă luăm un sistem cartezian de coordonate Oxyz şi reprezentăm 
mulțimea punctelor P(x,y,z) unde z=f(x,y) obţinem imaginea geometrică a 


funcţiei sau graficul său. De data aceasta graficul e o porţiune dintr-o suprafaţă. 


Z a 
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. . _. 2 3 ME. 
Exemplu: Să se reprezinte graficul funcției z =4/x? +y? . 


Soluție: Ridicăm la pătrat egalitatea şi obținem z? = x? + y?. Dacă facem 
aici x = 0 (intersectăm graficul cu planul yOz) obținem z = + y adică două drepte în 
acest plan. În mod similar intersectăm ` graficul cu planul xOz şi obținem  dreptele 
XS. 


Graficul este un con infinit cu vârful în originea sistemului de coordonate. 


Fig. 5.2 
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5.3. LIMITĂ ȘI CONTINUITATE 


Să considerăm o funcţie f :D C R? >R şi (xo, Yọ) un punet de acumulare 
al lui D. Aceasta înseamnă că există cel puţin un şir de puncte (x Yy en din D 


care converge la (Xo, Yo). 


Definiția 5.3.1. Funcția f(x,y) are în (Xp, Yọ) limita L, dacă oricare ar fi 
sirul (Xk Ya dhen ED, (x: Yen ~ (Xo Yo) AXL YK) = (Xo Yo) rezultă că 
f(x, Yp) >l. 

Vom scrie: lim f(x,y) =£ sau lif (gey) = g 

X—X9 


(x,y)—(X0:Y0) 
Y Yo 


Această limită se numeşte limită globală deoarece x şi y tind simultan şi 
independent către Xo respectiv yọ. Dacă xo şi y tind succesiv la xo respectiv yo adică 


dacă £, = lim (lim f(x, y)), 45 = lim (lim f(x,y)),, se obțin limitele iterate sau 
X——Xo Y—Yo Yy—Yo X—Xo 


parțiale. Are loc: 


Teorema 5.3.1. Dacă limita globală există atunci şi limitele parţiale 


există şi are loc egalitatea: 

l =A NE (5.3.1) 

Reciproca nu este adevărată. 

Demonstrație: Pentru a justifica rezultatul este suficient să considerăm mai 
întâi şiruri generale (x,y, en — (Xos Yo) şi apoi să le particularizăm în şiruri de 
formă particulară (xp, Yo)ken ~ (X0 Yo) $i (XorYk)ken (Xo Yo). Ceea ce e 


adevărat pentru general se păstrează şi pentru particular. 
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Faptul că reciproca e falsă se poate proba prin următorul exemplu. 


Să se calculeze limitele iterate ale funcției f(x, yz a. y în punctul (0,0). 
XX 


ł& = im tim | mă = liml=1 
„AO y—0 x AX x—0 xXx x—0 


N = mfi 0) im pm 
ISORO hei | Ye Ye yO 


Cele două limite nefiind egale nu există limita globală a funcției în origine. 


Observaţii: 
1°. Dacă £, = £, nu rezultă că limita globală / există, ci doar că dacă 
aceasta există este egală cu valoarea lor comună. 


20, S-ar putea ca în punctul (xo,yo) funcția f să nu fie definită adică să nu se 


poată calcula f (x0,y0) cum este în exemplul precedent. 


Definiţia 5.3.2. Funcţia f (x,y) este continuă în (xo,yo) dacă: 
lim f(x,y) = f (x0 Yo) 
X— Xo 


Yo 


Ca şi în cazul limitei globale, acest tip de continuitate se numeşte 
continuitatea globală. Se poate şi aici pune problema continuității în raport cu 


fiecare argument când celălalt este fixat, adică continuitatea parțială. 


Exemplu: Să se studieze continuitatea globală şi cea parțială pentru funcția: 


xy 
f(x,y)= pietele di O 
0 „(x,y)= (0,0) 
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Calculăm limitele parţiale: 
4 = limf(a,0) = lim0=0 
x—0 x—0 
l, = lim f(0,y)= lim0=0 
y—0 y—0 


Se observă că £, = f(0,0)= £, =0, deci funcția este continuă în raport cu 
fiecare argument. AFI 

Limita globală nu există pentru că luând şiruri (X,, VE) dle — (0,0) în care 
Yk = mXx, m=0 obținem: | | 


Xr “MĂ m 

f (xg, Yy) = =E k => A 
kk 2 pii z9 
i T O e l+m 


deci f(x, Yk) = care depinde de m. 


1+ m? 


Aşadar, continuitatea globală nu există. 
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5.4, DERIVATE PARȚIALE ȘI DIFERENȚIALE 
5.4.1. DERIVATE PARȚIALE DE ORDINUL ÎNTÂI 


Să considerăm funcţia f:D—R ŞI (Xo, Yo) un punct din interiorul mulțimii D. 


Considerăm de asemenea un punct (x, y) din D diferit de (xo, yo). 


Definiția 5.4.1. Se numeşte derivată parţială a funcţiei f în raport cu 


variabila x în punctul (Xo, Yo), numărul real 


lim, £(% Yo) -f (xo Y0) (5.4.1) 
X—X9 X —X9 


dacă limita există. Dacă aceasta este şi finită, spunem că funcţia este derivabilă 


parţial în raport cu x în punctul considerat. 
3 m Li Of / 1 
Notăm valoarea limitei cu geyo) sau f/f (Xo Yo). În mod similar se 
x 
introduce derivata parțială și derivabilitatea parțială în raport cu y în (xo, yo): 


ay oyo) = f! po Sa SA fCxory)— f (xo: Yo) 
>o Y—Yo 
Observaţie, Derivata parțială nu aduce nici o noutate față de cazul 
funcţiilor de o singură variabilă. Practic, când se derivează parțial o funcţie de n 
variabile în raport cu una, toate celelalte variabile se consideră constante 


(parametri). 


Dacă f (x, y) este derivabilă parțial în raport cu x sau/şi y în toate punctele 


unei submulțimi D; a lui D, spunem că f este derivabilă parțial pe D}. 


În acest caz pentru orice punct ăia D, stabilim corespondenţele 


(Xos Yo) = Z koyo), CA ah a e Sa) care definesc două funcții noi 
X 


Oy 
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numite derivatele parțiale ale lui f în raport cu x respectiv y în Di. Ele se vor nota: 
EE sau f’ (x, y) respectiv 00 sau f(x,y). 
Ox i Oy 


Observaţie: Calculul derivatelor parțiale (a y), SE 9) se face 


Ox Oy 
folosind regulile obişnuite de calcul pentru derivata sumei, produsului, câtului, 


funcțiilor compuse etc. 
Exemplu: Să se calculeze E a 7 F N) dacă f(x,y)=xye”. 
Ox a FAT b} 
A of Kp | 
Soluţii: Calculăm pei y) considerându-l pe y parametru constant. 
of Si 1 xy Zyn S uaay mige uyi pasi 
dp) =(xy) e" + (yo? ), = ye” +xy(ye”) = (y+ xy" )e 
BE a a 

Calculăm ZX) considerându-l pe x parametrul constant: 

Oy nar | | + 
ðf I XY xy V xy . xy 2 xy 
mai =(xy)e + (xy)(e™), = xe +(xy)( xe") = (x+ x*y)e 
Observație: 


10. În exemplul de mai sus, simetria în raport cu variabilele x şi y permitea 


calculul celei de-a doua derivate parțiale, schimbând între ele variabilele x cu y în 
əf 
je s X, . 

z Y) 


2°, Toate considerațiile anterioare se pot extinde de la funcții de două 


variabile la funcții de n variabile fără nici un fel de dificultate. 
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5.4.2. DERIVATE PARȚIALE DE ORDIN SUPERIOR 


Să considerăm din nou f :D C R? —> R, derivabilă parțial în raport cu x şi y 
pe D. Cele două derivate parţiale f; (x,y) = g,(x,y) şi f(x,y) = g2(x,y) pot fi la 
rândul lor derivabile pe o submulțime D, a lui D. Apare natural să considerăm 


derivatele lor parțiale drept derivate parţiale de ordinul al doilea ale funcției f 


0g, O Of n 

ZOL anl a i IM 

r(x y)= oo (x, j= ga ON > fa Ys 
9 = bi 9 f 

2 (x,y) zi -P (x y|- N (x,y) & f3 (x,y), 

g- ð | of 92f P 

Phe h A = — | — $ = ——— î i 

ay (x,y) ži Oy (x ») Dx0y (x i TRY), 


e RRE N Ec pat i 

Pentru funcția de două variabile f(x,y) obținem astfel patru derivate de 
ordinul al doilea. Dintre acestea se remarcă f7 (xy) şi f Cay) numite derivate 
parțiale mixte de ordinul al doilea. | | | 


„Natural se pune întrebarea dacă între cele două derivate există legături. 


Exemplul următor sugerează că răspunsul este afirmativ. 
Exemplu: Să se determine derivatele parţiale de ordinul al doilea ale 
funcției fD CR? >R, f(x,y) =x2y+2xy+1. 


Soluţie: Obţinem succesiv: 


TLAN ; Aue fai 
Ox LEWE 
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df  Ə(əf „0 9 (0 
| = ]=2y =2x +2 
3x? "0 (ž]- ' ðxðy ax HE) j 


i 23] aer A SE | 08 A 
Oyox Əy(Əx dya = Gy Edy 


Se observă că cele două derivate parțiale mixte pa E, ZA sunt egale. Nu orice 


funcţie f (x,y) are această pr oprietate. Fără demonstrație enunțăm următoarea teoremă: 


Teorema 5.4.2. (Criteriul lui Schwarz) 


Dacă în domeniul D, C D funcţia f (x,y) satisface condițiile: 


Of of 
a) are derivatele —, — 
Ox '0y” 
-2e 2 
b) derivatele mixte de ordinul al doilea — ax ay dA En sunt continue, 


atunci derivatele parțiale mixte de ordinul al doilea sunt egale. 


Observaţie: Funoţia f y) E di: spre anterior ihe condițiile 
criteriului lui Schwarz. 

Fiecare derivată parțială de ordinul al doilea poate fi la rândul său derivată 
parțial în raport cu cele două variabile. Se obțin astfel derivate parțiale de ordinul al 


treilea şi procesul poate continua. 


5.4.3. DERIVATELE PARȚIALE ALE 
FUNCȚIILOR COMPUSE 


Să luăm funcția f:DCR?—R şi să presupunem că la rândul lor, 


argumentele x şi y ale lui f sunt funcţii de variabilele u,v,w: 


f(x,y) = fa (u,v,w), y(u,v,w))= 0(u,v,w) (34.2) 
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Ne interesează derivatele parțiale ale funcției Ọ în raport cu cele trei 
variabile u,v,w. Un calcul nu foarte complicat care foloseşte. definiţia derivatei 


parțiale conduce la formula de calcul a celor trei derivate parțiale: 


op _9f Ox əƏf dy 


„du Ox 0u dy "du! 

' ap Of Ox Ox , (Li sară 
"Ov. 9x dv Oy Ov 

Op _09f Ox Orad: Oy 

ôw ôx ôw Oy ðw 


Exemplu: Să se afle dei Vatale parțiale ale funcției f(x,y) = xy2 + x2 y, 
cunoscând căx=u+v,y=u-v. 


Soluţie: Notăm e(u,v) = f(x (u, v), y (u,v)) 


dp _ of „9x f, Yog + 2x 1+ (2xy +x? Ja =? +4xy +y?, 
ðu 9x ðu ðy du 


dp ðf əx Əf ðy 
St RR 2 2xy)-1 Tsy i 
El i d Siria = (y2+ xy). +(2xy+x?). y? — x2 


Observaţie: Formula (5:4.2) se poate generaliza când f are n variabile, iar 


fiecare este funcţie de m variabile. Evident aplicarea formulei este posibilă când 
toate funcțiile implicate sunt derivabile parţial. 


5.4.4. DIFERENŢIALELE FUN CȚIILOR 
DE MAI MULTE VARIABILE 


Fie f:DCR:—R » (X0Yo0) şi (x,y) două puncte din D. Să notăm h = x - xo, 


k = y - yo. Numim creşterea totală a funcţiei în punctul (x0,y0) cantitatea Af( X0 Yo) 
definită prin: 


Af (Xo: Yo) =f (x, y)-f (Xo Yo) (5.4.3) 
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„Definiţia 5.4.2. Funcţia f este diferentiabilă în (x0,y0) dacă există două 
constante A şi B independente de h şi k şi două funcţii AXo Yo bike Aa voit Ry” 
astfel încât: 


Q(X9» Y9;0,0) =0, B(x9,Y0:0,0)=0, lim (xo, Yosh,k) > lim B(xo.Yo,h,k)=0 
» 1> 
k—0 k—0 


Af (Xo: Yo) = Ah + Bh + ho(xo, Yoh; k) + k: Po yosh,k) (5.4.4) 


Se poate demonstra următoarea teoremă: ` 


Teorema 5.4.3. Dacă f(x,y) este diferențiabilă ih (x0,y0) atunci ea este 


A 


derivabilă parțial în raport cu x şi y în acest punct şi 
əf əf | mandinkai 
— (Xo: Yo) = A, — (Xo, Yo) = B 
ax o Yo) ay. o Yo) 


Definiția 5.4.3. Expresia Ah +Bk se numeşte diferenţiala funcției f în 


punctul (x0,y0) şi se notează: 


df (x0,yo) = Ah + BK. (5.4.5) 
Dacă funcţia f(x,y) este diferenţiabilă în (x0,y0) atunci diferenţiala sa 
df(x0,y0) e dată de 
Of of -SN 
df (Xg, Yo) = (Xo Yo)h Wi Oyok (5.4.6) 
x dy 


Cum dx = h, dy = k obținem diferenţiala de ordinul întâi a funcţiei f (x;y) 


în forma finală: 


of of 
df (x0, Yo) = zy A0 Yo)dx + Dr Xonoddy 


SI rii Să se determine diferenţiala de ordinul întâi “a funcţiei 
dia jo ran punctul (1,1). 
n et. CU X ra leg Pr mud, X 
a om pai a A Vl. 
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Soluţie: f’ (x,y) = y? + 6x2y f.(ib=7 
~ x 3 Li LE 
f(x,y) = 2xy + 2x fd, =4 


Obţinem: df (1,1) = 7dx + 4dy. 


Observaţie: Diferenţiala de ordinul întâi a funcţiei f (x,y) serveşte la 


aproximarea creşterii totale a funcției TREST (Xo, Yo) conform formulei (5.4.4). 
Dacă f (x,y) este diferențiabilă + în toate PEC nyo) unui die lu 
D, € D, puten defni O > corespondență: 
ia a df (xo Yo): (Y) (X0 Yo) ED, 
care se numeşte diferenţiala lui f şi se scrie df(x, y). 


Ea poate fi la rândul său diferențiabilă. Diferențiala sa se numește 
diferenţiala de ordinul al doilea a lui f şi se notează d°f (x,y). Un calcul elementar 


conduce la formula d?’f (x, y). 


d? fayd RE (dt +22 


y orik eya. K? 


a 


Exemplu: Să se determine d?f(1,1) pentru funcţia: 


f(x,y) = xy + 2x2y7 


Soluție: 

Of s "OT 22 

— = y+ 4x — = x + 6X 

3x Pg y Oy ă 

92f 3 92f > 2f e] 
FER =1+12xy , S r 
1% i d ðxðy zi a 
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T an=4 : a apei ZEA =12 şi în fina 


d?f (1,1) = 4h? + 26hk +12k?. 


Observații: i | 

1°, Procesul poate continua obținând diferențiale de ordin superior lui 2. 

g5. Diferențiala de ordinul al doilea este după cum se poate vedea din 
exanp de mai sus o formă pătratică de variabilele h,k. Când variabilele sunt 
XisX2s:. Xp această diferenţială devine o formă pătratică de n variabile. Ea j joacă 
rolul derivatei a doua la funcții de o singură variabilă, adică serveşte studiului 


convexităţii unei funcții de n variabile. 
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5.5. EXTREMELE FUNCTHLOR DE 
MAg MULTE VARIABILE 


5.5.1. FORMULA LUI TAYLOR PENTRU FUNCŢII 
DE MAI MULTE VARIABILE 


Se cunoaşte formula lui Taylor pentru funcții de o variabilă care are 


derivate până la ordinul k+1 într-o vecinătate a punctului:xo: 


(x= sao 


F=f + z faini de ro Le 
(x — alea (k) (Xe a) (k+1) 
Enee TX DT REN 


0 fiind un punct intermediar lui xo şi x. 


Această formulă permite aproximarea valorii f în x cu ajutorul valorilor 
acesteia şi a primelor sale k derivate în punctul xo, dacă ultimul termen (restul de 


ordin k) este neglijat. 


Căutăm o extindere a formulei la cazul funcţiilor de două variabile după 


care, natural, extinderea se poate face la n variabile. 
Fie aşadar, f:D C R? > R şi (X0,y0) un punct interior din D. 


Presupunem că (x,y) este un punct oarecare din D şi notăm h = x — Xo, K=Y—Yo. 


Teorema 5.5.1. Dacă f (x,y) are toate derivatele parțiale până la ordinul 
m în (X0,Y0), iar derivatele parțiale de ordin m+1 există într-o vecinătate V a 


lui (x0,Y0) atunci oricare ar fi punctul (x, y)€ V „are loc formula: 
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f(x,y) =f(xo Yo) FR E (Xo Yo)h ră go yakla 
CIRIE F. | | V 
Zany ta AEN £ a (3.31) 
(m+1) 


O 
Z koyon y Soyo f+Ra 


m! mi 


unde : 


E (ii Vad dis Z tao yo f reprezintă puterea formală , „i“ a expresiei 


diferențiale din paranteză aplicată funcţiei f, iar R. numit restul de ordin m în 
formula lui Taylor este exprimat prin: 


(m+!) 


R, = PoE ODI (5.5.2) 


9 9 i 
(o +0h, yo + k)h hgg ta FOR Yo Gi 


5.5.2. TEOREMA LUI FERMAT PENTRU FUNCŢII - 
DE MAI MULTE VARIABILE 


În teoria diferențială a funcțiilor de o variabilă, teorema lui Fermat afirmă 
că într-un punct de extrem local din interiorul domeniului de derivabilitate al 


funcției, Du, derivata: f” se anulează. Încercăm aici o extindere la funcţii de două 


variabile care se poate apoi transpune la funcţii de n variabile. 

Să considerăm f :D C R? >R şi (Xo YJ ED. 

Definiția 5.5.1. Punctul (x, yọ) se numeşte maxim (minim) local pentru 
f(x,y) dacă există o vecinătate V a acestuia, astfel încât (V) (x,y) €V are loc: 


f(x,y) < Aay) (£(x,y)2 f(x0:Yo)) (5.5.3) 
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Teorema lui Fermat poate constitui o metodă prin care se determină 


punctele de extrem sau posibilele puncte de extrem. 


Teorema 5.5.2. Dacă f(x,y) are derivate parţiale finite în raport cu x şi 


y în punctul de extrem local (xg, yọ) din interiorul lui D, atunci: 


Of. Of 
— (Xo Yo) = — (X0:Y0)=0 (5.5.4) 
ox Oy 


Demonstraţie. Considerăm (x)= f(x,y) şi v(y)= f(x9,y). Cum 
(Xo: Yo) este în interiorul lui D, (vezi fig. 5.3) xo se află în interiorul lui (a,b) iar Yo 
în interiorul lui (c,d). Se poate aplica teorema lui Fermat pentru p(x) cu xo punct 
de extrem local şi (y) cu yo punct de extrem local. 


Obţinem: 


d Of dp Of 
(0) = 0 (tor0) =0 g 030030) =0 


ceea ce demonstrează teorema. 
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Soluţiile sistemului (5.5.4) se numesc puncte critice sau puncte staționare. 
Reciproca nefiind adevărată, punctele critice rămân posibile puncte de extrem. 


Condiţiile (5.5.4) pentru punctele (Xg, Yo) sunt doar condiţii necesare de extrem. 


Exemplu: Să se determine posibilele puncte de extrem ale funcţiei: 
f(x,y) x Fy —3xy +15 


Soluţie: Calculăm Mei 
Ox Oy 


s- 3X2 ae me = 2 ay 
ox Oy 
Rezolvăm sistemul jai =) dă, = 
ox Oy 
3x2 —3y=0 
2y — 3x =0 


Obţinem soluţiile (0,0) şi 27) care sunt posibile puncte de extrem.. 


5.5.3. CONDIȚII SUFICIENTE DE EXTREM 
(DETERMINAREA EXTREMELOR) 


Am văzut în cele de mai sus că un punct de extrem este punct staționar, dar 
din teorema lui Fermat nu putem obține şi reciproca. Căutăm aşadar condiţii 
suficiente de extrem prin care să putem aprecia dacă un punct staționar este sau nu 


este punct de extrem. 


Să remarcăm mai întâi că formula lui Taylor (5.5.1) se poate scrie şi sub forma: 


| ÎI i» 2 
Pi (aprig) 74 op EAT T f(y Ra (3. dai 
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Dacă m = 2 atunci 


(x, y= foy dt TER, (5.5.6) 


Să presupunem acum că (Xo Yọ) este un extrem local. În virtutea : 


condițiilor (5.5.4) rezultă că 
doo) = (tony )h-+ Se ayok = 0 În plus dacă X—X9=h şi 


y— Yo = k sunt suficient de mici, Rz poate fi neglijat (acolo h şi k apar la puterea a 


treia). Rezultă că pentru (x,y) într-o vecinătate suficient de mică a lui (Xo Yo): 


sen[f, est (0; Y0)]= send? Hre 


ori astfel: 


Teorema 5.5.3. Dacă forma pătratică d*f (X9+Y9) este pozitiv definită în 


punctul critic (x9,yg) atunci acest punct este minim local. Dacă în acelaşi 


punct d2f (xo, Yo) este pozitiv definită, punctul este maxim local. 


Condiţiile în care o formă pătratică este pozitiv (negativ) definită au fost 


prezentate anterior. Pentru prezenta formă pătratică, matricea este: 


f af 
gaz So Yo) Da Y0) 
Ii = y (5.5.7) 


OR of 


care se numeşte matrice hessiană şi se poate generaliza la funcţii de n variabile. 


Suntem conduşi astfel la următoarea teoremă: 


186 MATEMATICI APLICATE ÎN ECONOMIE 


Teorema 5.5.4. Condiţia necesară şi suficientă ca punctul critic (xp, yọ) 
să fie punct de maxim este ca lanţul minorilor diagonali A,,A, ai matricei 
hessiene să fie pozitiv. Condiţia necesară şi Sie lentă ca (Xo, Yo) Să tie punct 
maxim este ca lanțul minorilor principali A,, A, ai matricei hessiene să fie 
alternant cu (A, <0). În orice altă situație, (x9,yg) nu este punct de extrem. 


sau nu ne putem pronunța. 


Exemple: 
10. Să se determine extremele funcţiei: | 
f,y)=x5+y2 —3xy +15 


Soluţie: Determinăm punctele de extrem local printre punctele critice, adică 


printre soluţiile sistemului: 


i e nisa OT. yi 4 i b EA Me 
—=0, —=0. Am văzut anterior că sunt două puncte critice 
Ox Oy 
9 27 i $ 
(0,0), a . Fiecare poate să fie sau nu punct de extrem. 


of. Of ETE 
P a aa EEE AA 
ox OxOy Oy 


Matricea hessiană: 


6x —3 
H(x,y)= tn 3| 


O calculăm în fiecare punct critic: 
0 —3 l 
H(0,0) = 1 A =0, A,=—9 
(0,0) e y ȘI | 1 ir 


Punctul (0,0) nu este punct de extrem. 
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27 
sAr 
n22- ZE Pest Al LE a za 
TĂI ie ia DE vă 


| 9 27 po 
„ Punctul. TEI este un punct de minim. 
20. Să se determine extremele funcţiei: 
ERE PER 


Punctele staționare (critice) sunt soluţii ale sistemului 


a 342%, Dap o = 2x; =0, O — asha „50, 
OX, Ox, Ox3 


Obţinem un singur punct critic (0,0,0). Acesta este şi singurul posibil punct 


de extrem. 
Ar EL af fatal 6 ao unul ina 
ORE i AORN OR pe uda aa 
Harpuin tgif n = pr VO 
"0X,0X3 0x30 i 0x0x; 0x30x2 


Matricea hessiană devine 


SĂ d. 
H(x x33) =1 0 2 0 
mi md. 2 


şi are lanţul minorilor principali 
A Sg Rpa 16 


_ Rezultă că (0,0,0) este un punct de minim local. 
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5,6. EXTREME CONDIȚIONATE 


Există situaţii când se caută extremele unei funcții: f:DCR'—R iar 
variabilele x,,x3,...,X, sunt supuse unor condiții date, exprimate analitic printr-un 
număr de relaţii de legătură: 

tz îi i0 


ATE TER Pa) 


Emp Xas NART 3 


| 
o 


(5.6.1) 


li 
© 


unde m<n-—l. 


O astfel de problemă se numeşte problemă de extrem condiționat sau de 
extrem cu legături. 


De data aceasta ne ocupăm de cazul general pentru că în cazul n=2 


problema devine prea simplă ca să fie reprezentativă. 


Definiția 5.6.1. Punctul X} = (Xo X025:::»X0n) Care satisface condițiile de 
legătură (5.6.1) se numeşte de maxim ( Aim ) condiționat al functiei f dacă există 
o vecinătate V(xo) a sa, astfel încât oricare ar fi x = (x, X33., Xp) din V(x )ND 
este satisfăcută inegalitatea: 


FOX) fo) (Ea) fo). (5.6.2) 


În legătură cu restricţiile (5.6.1), trebuie observat că numărul lor nu poate 
depăşi n — 1, deci dacă f are două variabile, atunci este posibilă o singură relație, 
dacă f are trei variabile, atunci sunt posibile maxim două restricții ş.a.m.d. În plus 


e 4 i og; EA b a. . x nb 
matricea cu elementele m; = d. având m linii şi n coloane, trebuie să aibă rangul 
X. 
j 


m, ceea ce asigură independența restricțiilor. 
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Ca şi în cazul extremelor simple (libere) trebuie determinate condițiile 


necesare şi apoi cele suficiente pentru extremele condiționate. 


Metoda utilizată frecvent este metoda multiplicatorilor lui Lagrange care 
transformă. o problemă de extrem condiționat într-o “problemă de 'extrem liber 


pentru o funcţie nouă numită funcţia lui Lagrange. 
LXi Xa Naso m) = f(X Xs X Azi Caraș: Xa) (5.6.3) 
=l 


Constantele A; (i =1,m) se numesc multiplicatorii lui Lagrange. 


Teorema 5.6.1. Dacă funcția f(x{,X5,...,X,) admite în punctul 


Xo = (Xo X02» Xġn) un extrem condiţionat relativ în care 


gi (X0)=0, (i= l,m) , atunci în acest punct au loc relațiile: 


papei AN Aia e e d (5.6.4) 
OX, ôx, Ox 
A POPP - i k GES 


Gnu 1 rare KAAN 


Am obținut două sisteme care împreună au n + m ecuaţii şi n+ m necunoscute. 


De fapt (5.6.5) este chiar (5.6.1) deoarece = =g; (Xp. X,)=0, i=l, 


Fie RE = (Xop X02»: » Xon» A01 oz» Aom). O soluție a sistemului (5.6.4), 
65. 6. 5), Acesta este de fapt un punct staționar (Critic) pentru funcția lui Lagrange L. 
Deoarece 8; (Xx9)=0, i= l,m rezultă că soluția sistemului (5.6.4) este punct 
siah Mar pentru f care satisface restricțiile impuse. 


„ Căutăm acum o condiție suficientă de extrem condiționat pentru xo. 
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Să considerăm: 
EX) = L(x; 3 X27113 Xni Ao A022: A0m) 


în care Ag; (j= 1,m) sunt fixate. Diferențiala a doua a acestei funcții va caracteriza 


starea de punct de extrem condiționat a lui xo. 


PLx) = 5 Ao? 9 O LX) 


dx, dx, A | (5.6.6) 
jel Ox,0X; 


Variabilele dx; (j=1,n) nu sunt independente, ci satisfac consecințele 


diferențiale ale condițiilor (5.6.1) 


$ Tau oda t ap (oda sk. ca: ui (%0)dxa =0 

öz d | 

dea ot, eo e Xa DE (xx, =0 (5.6.7) 
Bm, = (xo)dx; rm (xo)dx; +.. E (xo)dx, =0 

Ox A Ox 


n ' 
Sistemul omogen (5.6.7) are soluții nebanale, iar rangul matricei sale este m. 
Aşadar m dintre necunoscutele dx, (fie ea dx,,dx,,..., dx, ) se exprimă prin 


combinaţii liniare de celelalte n-m. 
dx, = h dx im (5.6.8) 


Înlocuim (5.6.8) în (5.6.6) şi Site în iii o formă PATA SĂ 9 în 


“variabilele dx „ps dX m42»: dX 


m+l n 


Q= Š 8;dx; dx, (5.6.9) 


ij=m+1 


În baza proprietăților: generale ale formelor pătratice putem da acum 


condiția suficientă de extrem condiționat. 
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Teorema 5.6.2. Dacă lanţul minorilor principali (diagonali) ai matricei 

(B; ) este pozitiv, atunci xọ este minim local condiționat. Dacă lanțul minorilor 
1) a 

principali ai matricei este alternant cu (A, <0), atunci xp este un maxim 


condiționat. 


In orice altă situaţie xo nu este extrem condiţionat sau nu ne putem pronunţa 


asupra acestui fapt. 


Exemplu: Să se Jidin extremele iinei f(x,y) = xy care Tia 
condiţia x+y=6. dz 

Soluţie: Funcţia lui Lagrange L are expresia: 

L(x, y;A) = xy +A(x+y-—6) 


Punctele staționare ale acestei funcții sunt soluțiile sistemului: 


| OL OL 
—— +; =0, — = x FAF 0,... 5x + —6=0 
Ox j Oy OA ş 


Obţinem soluţia (3,3;+3). Posibilul extrem condiţionat al lui f este Xg =(3,3). 


Notăm Xp = (3,3;—3). 
Calculăm dL(%,) 


3LF) 


A 9 L) 
L(3.) = L(3,3;—3 dxd 
(Xo) = L( )= PRE Dx0y yo 


see îi die eh d ai ta 2. DL) ay 
ðy? 


= 0fdxA+ 2- 1-dxdy +0- dy? = +2dxdy 


Diferențiem condiția x+y =6. 
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Obţinem. : 
~ dx +dy=0 
adică: 
dy = —dx 


Forma pătratică Q este dată de : 
Q = +2dx(—dx) = —2(dx) 
care este evident negativ definită. Punctul xg = (3,3) este un maxim local care 


satisface condiția x + y = 6. 
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5.7 APLICAȚII PRACTICE 


I. STUDIUL FUNCȚIILOR. LIMITE. 
CONTINUITATE. DERIVABILITATE 


Fie f:DCR'—R şi P,(xp,x3...x0)e D. Pentru a studia limita (globală) 


şi continuitatea funcţiei în P, se procedează conform următorului 


I.1. Algoritm de lucru: 
l. Calculăm limitele parţiale în P,. Dacă acestea: 
a) există m sunt diferite; atunci limita globală nu există; 
b) există şi sunt egale (sau nu există), limita globală poate exista sau nu; 
c) existența limitei globale implică unicitatea ei; 


2. Pentru studiul continuității funcției f în punctul P e D (există deci f(Po)), se 


2.1 /=F(P,), funcţia f este continuă în P,. Când /zf(P,), P, este punct 
de discontinuitate pentru f. 
2.2 4+ f(P,) dar limitele parţiale există şi sunt egale cu f(P,), funcţia este 


continuă în raport cu fiecare variabilă, în parte, când variabilele 


complementare sunt fixate. 
Observaţie: Pentru calculul limitei globale se folosesc şiruri (P, „e din 


R” şi se calculează l= lim f(P,). Dacă f :D—R, şi DCR?, atunci se folosesc 


şiruri (x,y) şi se calculează lim f(x,y). 
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Pentru a arăta că / nu există, este suficient să trecem la limită utilizând şiruri 
. ? | ae. ke b 
particulare de forma: y, =mx,, yY, =mx} etc. (me R, oarecare). Dacă pentru 
diferite valori ale lui m, se obțin limite diferite, atunci limita globală 44 nu există. 


3. Pentru calculul derivatelor parţiale se derivează funcţia în raport cu o variabilă, 
considerând celelalte variabile drept constante, după regulile cunoscute. În 


calculul derivatelor parţiale de ordin doi, se ține seama de cfiteriul lui Schwarz. 


1.2. Exerciţii şi probleme propuse: 
l. Determinaţi domeniul de definiție şi mulţimea valorilor pentru următoarele 


funcții f:D—I, ISR: 


a) f(x, y)=41=x7- y7: b) FRYE aE py 
c) TE REEE | d) f(x y)=vl-x? +yl-y*. 


și y% 
R? cay D= (x, yje R?/x24+ y2< 1} adică discul de centru O şi rază egală cu 1. 
IeR, 
b) D=R,I=K; 
Cha D=Rå vi, y)e R"/x2+y2= 1} - mulțimea punctelor din plan cu 
excepția celor situate pe cercul de centru O şi raza 1, I=R; 


d) D= [>71 1]}x[-1, Ilte Re 


2. Folosind definiția arătați că: 


a 2 E, 
DR: S 


| XI = 3 
(x,y}(3.%) di 2 
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„4 4 
X+y 20 


c) 3 > 
(x.y}>(0.0) x * + y“ 


Indicaţie: 


a) Se consideră (x. yi> (1,3), şi se calculează lim f(x,y); 
b) Se consideră (x,,y,) cu limx, =3, limy, =%; 


c) Se foloseşte majorarea: 


3. Arătaţi că următoarele funcții nu au limită globală în O(0, 0):. 


2xy 2 + 4x 
f(x,y)=——=; | , 
a) f(x,y) Zap 


Indicaţie: 


a) y, =MX, X, 0, limx, =0; 


b) y, =mx,, X, #0,limx; Fag 


4. Să se cerceteze limitele parțiale şi limita globală în O(0, 0) pentru următoarele 


funcţii: 
Sia Fi 
a) f(xpye XI ATI AR: a) 01021 4muexistă; 
X.Y 
b) (x,y) = xsin b) £, =0; £, nuexistă; /=0; 
REA 


c) f(x,y) (x + y sint c) €,,£, nu există; (=0. 
DX 
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5. Calculaţi limitele globale: 


292 


. M 
l ~ R: 132. 
pu (nt) Îx2y2 +1 md a) 
D N E D LaF: 
(x.y)(0,3) x | 
(x.y)=(0.0) Xy 4 
d) li DR. dute A : d) /=0 


ım 
(xy He) x? + y? 


6. Studiați continuitatea funcțiilor: 


| x+y., 
a) f(x,y)= TITANA 


+ 
n 
© 
© 
N” 


x +y? 


b) f(xy)=i mky) (x,y) (0,0) 


pă ; (,y)=(00) 
R: a) continuăpe R”1((0,0)) 


b) continuă pe R? 


7. Calculați derivatele parțiale de ordinul întâi ale funcțiilor: 


a) f(x,y)=x°y+2xy +1; b) f(x,y)=e*”; 


c) f(x, y)=(x? +y hretg%; d) F(x, y,2)= xyz +; 
y = 


In(xy) 
a 


e) f(x, y,z)= 
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— a —— 


R: a) A aysay; Enen 
y 


Ox 

of CT, af r R: | 
b) —=2xe* 7; = = 2ye ™ 
) Ox dy i 
c) za rr S EE 

Ox y dy y 

of 2e Of | 3 l of 292 y 
d) —=yz'; —=2xyz +-—; — = 3Xy z -= 
EP, Oy i 92 aN z? 
y "ae h A + u s, 

Ox xz Oy yz Oz ză 


8. Arătaţi că funcțiile: f(x,y)=Ina]x2+y2şi g(x,y,z)=4x? +y? +z? satisfac 


relațiile: 


i  92f (22) d) SA 
2) Ox? ðy? Ox dy Oz 
9. Să se calculeze diferenţialele de ordinul întâi şi de ordinul al doilea iaiu 
funcțiile de mai jos, în punctele indicate şi creşterile specificate: 


a) fy) P, (2,1), h = dx =0.1, k = dy =0.02; 


b) f(x,y)=kx%y (funcția de producție Cobb-Douglas) ; P (1,1) cu 
h 20.01 şi k = 0.05: 


10. Să se calculeze df şi d°f pentru funcţiile: 
a) f(x, y)= e” cosy; 


b) f(x XX) RR pla i 
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2. 


f(x, y)= e* (cos ydx — sin ydy),sau 
R: a) <df(x,y)=e*(hcoy-—ksiny) ; 
d’f(x,y)=e* (h? cos y — 2hkssin y — k? cos y) 


df (x z ENN X2X30X, + XX dtn A p, sau 


b) 4df(x, k a ADR PEA E Aa | 
df(x, xx, ]=2(x,h,h, +x,h,h, + x,h,h;) 


11. Să se calculeze df în punctul P, (3,4,5) pentru funcția: 


ix yae 
(x,y,z) Y 


R: af (P,) = df(3,4,5)= 22 (3dx + 4dy =5àz) sau 
E 
df (3,4,5) =— (3h+ 4k - 51 
G4) = gh -3) 
12. Să se arate că următoarele funcții: 
a) f(x, y)= Ink? +y?) 
b) f(x, y) arctgY 
aS TR 


02f911036 și „i 
verifică relaţia PAE Jy? =0 (ecuația lui Laplace) 
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II. PUNCTE DE EXTREM PENTRU FUNCŢII 
DE N VARIABILE 


II.1. Cazul fără restricţii (extreme libere, necondiționate) 
Pentru a determina punctele de extrem local sau relativ, liber, ale unei funcții 
F(x Xa) se foloseşte următorul: 


I].1.1 Algoritm de lucru 


1. Se calculează derivatele parţiale de ordinul întâi: 
„d 9f dr 
Ox, Ox, "Ox 


n 


2. Se rezolvă sistemul: 


(%) = =0,— =0,...—=0 


Soluțiņe sistemului (*) sunt punctele critice ale funcției f. Fie acestea Pi i EE E M 


3. Se calculează cele n? derivate parțiale de ordinul doi: 


al a 
Anu OV d 
Observaţie: ‘Dacă derivatele mixte satisfac criteriul lui Schwarz „ele sunt 


n(n+1) 


egale două câte două și se calculează efectiv doar „ derivate parțiale de 


ordinul al doilea. . : 


4 Se construieşte matricea hessiană ataşată lui f: 


d'f ðf 
| iu e ăi Em 
Help, Xa) = 
9"f Of 
axa, lili: Ox: 


care, în baza observaţiei de la pct. 3), va fi simetrică (H, = HT). 
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5. Se calculează H, în primul punct critic P „adică: - 
H,(P)=[A,} i j= 5n , unde Aj= EEP) 


Ox, .0X, ji 


6. Se calcula minorii diagonalei Di, Das ues D, ai matricei H (P). 


7. Pentru a vedea dacă P, este punct de extrem, se sheldăză la condiţiile suficiente 


de extrem. Astfel, dacă: 
D A 0, Aa dz Aa > 0, P, este punct de minim local, 
b A, <O4, >04; <0,A, >0...(-1YA4,>0, P, este e) de maxim 


local; 


c) A +0, ln ŞI nu se încadrează în situaţiile a), b) atunci P) nu este. 


punct de extrem local, şi 
d) dacă cel puţin un A; =0, 1Si&n,nu se poate preciza natura lui P,. 
8. Se repetă etapele 5-7 pentru celelalte puncte critice P,,..., P, . 


Observaţie: În situația d) de la punctul 7) se aplică alte metode pentru 


determinarea naturii punctului critic, de exemplu: 
1. metoda lui Gauss de aducere la forma canonică a lui d°f (P, ); 
„2. se studiază direct semnul creşterii funcției Af (P, ) „ în vecinătatea punctului P, ; 


(AFP) =f + n xt + ho. xD ha) Ext, xx) 
[1.1.2 Exerciţii şi probleme propuse . 


|. Determinaţi punctele critice şi de extrem local pentru funcţiile: 


f:R:—>R 
a ; 
f(x, y)=x? +y? —3xy +15 
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f :(0,21)x(0,200)-—R 

f(x, y)= sin xsin ysin(x +y) 

f :(0,)x(0,%)—R | 
f(x, y)=x y a-x=yha>0 


pf: R? \{(00)}>R 
f(x, y)=xyln(x? +y“ 2 
f: R? R 
- f(x,y) > xik dz? TN 3z) 
` [f : (0, °)x (0, œ)x (0,20) => R 
l ZA ; 
y EE da Za 21 
Z 
i f: R= 
i 
AEN y)= ax” = bxy+xz+yz+l1l; ab,ceR. 
R ) -[P,(0,0) nu este punct de extrem local 
l a 
P, (1,1) este punct de minim 


b). -P;(2,3), P,(- 2,-3) nu sunt puncte de extrem local 


-pP (+ v2.2) P, -v2,2+ 43} P, (1-42,2- 3) 
nu sunt puncte de extrem local 


z -P,(l+v2,2+ V3} P „(+ /22- V3 )sunt puncte de minim 


=p, ( E Va 2)este punct de maxim relativ 
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î 
Pentru punctul critic P, (70, 70) nu se poate preciza natura 


d ; 
) (E z) este ponet de maxim local 


a a a 
e) E 2) singurul punct de maxim local.. Celelalte puncte critice nu 


sunt puncte de extrem; 


(a e (- zf punere e minim ca 
e Sa | Mai a 


(pa e = le )punčte de maxim local 
p | 


g) p(- 1, —2, 3) punct de minim ÁN | 
h) Ais l, ) punct de minim local; 
i) P(0,0, 0) este punct critic nu şi punct de extrem local; 


2. Un bazin de tablă de forma unui paralelipiped dreptunghic are volumul 
a? 
e ii Să se determine dimensiunile sale (x,y,z) astfel încât să se 


întrebuințeze cât mai puțin material (aria totală să fie minimă). 
TX ia? f z 
Indicaţie: z-înălțimea; V = Xyz = 3 ; S= 2(x2+ yz)+ xy — minimă cu restricția 
d 


a 
xyz——=0 
Lb 
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3. Să se determine un punct. M(x 4. Yọ» Zo )» astfel încât suma pătratelor distanțelor 


de la M la n puncte fixate P(a,,b;,c,), i=1,n să fie minimă. 


l la la 
R: a a, a E aa A 
Dinamani da Din i= 


4. O fermă de păsări foloseşte două tipuri de nutreţuri în hrana lor. Dacă notăm cu 
x, y cantitatea de nutreţ (în kg) dată păsărilor atunci creşterea în greutate (în kg) a 
găinilor, respectiv, cocoşilor este dată de relaţiile: 

(1) g(x y)=2x=y;e(xiy)=y-3 
Știind că profitul din creşterea în greutate (în milioane lei) este dat de relația: 

(2) f(g,c)=10-2(g -57 -(c-4} 
iar numărul de găini şi cocoşi este constant în ciclul de; producţie, să se determine 
amestecul optim de nutreţuri, astfel încât profitul rezultat să fie maxim. . 

Indicaţie: Se foloseşte regula de derivare a funcțiilor compuse (vezi [3]): 


df of. og. 9 de de ,0f. _09f 0g a QE dc 


+ —- 
Ox 0g 0x oc Ox 0y og dy ðc Əy. 


5. O societate care produce detergenți a observat că dacă cheltuieşte x, respectiv, y 
milioane lei pe “publicitate (scrisă sau la TV), profitul curent creşte după legea 


(funcția profit): 
f(x,y) = 4xy + 4x + 4y —3x2 — 2y2 -2 (milioane lei) 
Să se determine E, y, astfel încât creşterea profitului să fie maximă 
Indicaţie: Sa adiaz extremele funcției f şi se găseşte punctul P, (4, 5) care este 


de maxim şi fu = f(P,)=16 (milioane lei). 
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H.2 Cazul cu restricții (extreme condiţionate, constrânse sau legate) 


Fie f:DCR"—R, u=f(x,,X2o--X,), în care variabilele x, x,,.. X Sunt supuse 
la legăturile (restricțiile): 


R (xi )=0 
0 


pie E iatan 


„i=l,m, LENSA 
R (xpi x 
Pentru a determina punctele de extrem local ale funcției f care satisfac legăturile (*) 


se aplică metoda multiplicărilor lui Lagrange şi se foloseşte următorul: 


II.2.1 Algoritm de lucru 


l. Se ataşează problemei, funcția lui Lagrange: 
L(x in X A e) 2 E pia x) FAR x ea XE (x) 


m m 


unde A, e R sunt constante reale arbitrare (multiplicatorii lui Lagrange). 


2. Se determină punctele critice ale funcției L, precum $ ŞI air multiplicatorilor 


KR zh a rezolvând sistemul: 
Dle 24, Sl setei OL 
(2%) Ox, xə Ox 
OL Să "OERS i meri. 
— ERK 4..,xX1)=0,....— F(X: Xa) =0 
2 Kasy Da Ap oa O 
Observaţii: 


a) sistemul (**) are n+m ecuații şi n+m necunoscute (X1,X2,...Xn; A1, A25... An) 


b) sistemul (**) poate fi considerat ca fiind format din cele n derivate 
parţiale de ordinul întâi ale funcţiei L în i raport Cu X1, X2, ..., Xn la care se adaugă 


cele m restricții date de W 
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3, Se înlocuiese în expresia lui L valorile ;multiplicatorilor: A determinate la 
punctual 2, obținând: 
ră Li — p r 0 » , , > 
L(X jones daia dă ăi i» Aa) == ACI na s<-xXa) + An FiXa Xie Xn) + 
0 t . . “ai jé . j . P 
ME A Pol VE C e) 


funcție care depinde numai de variabilele th AP 


4. Se determină derivatele parțiale de ordinul doi ale funcției L de la etapa 3, adică: 
abiy j='ln 


FEP) 
0X,0X, 


5. Dacă P, este un punct critic se calculează aij = "=-aji şi nu se scrie 


diferenţiala de ordinul al doilea a lui L în punctul critic P}: 


"n 
z aij hi hj 


> 


AL(Pp = | 


6. Se diferenţiază legăturile din (*) şi se obţin expresiile a m dintre diferenţialele 
argumentelor în funcție de celelalte n-m, adică. 


dx, = fi, (AX puse dă (hu, = fi (heh) 
(x = f; (dX ascet AX), (ha 2 Fa (hash) 


po.croseeoocaseooonooooeesacoseseeenoecocneonasococeossoooeeneoeaeeoeose: 


dx, Dă PERI (rouă „dx ad (ha = = f haa m+?’ „h n)) 


A i <P 4 A f "DI, 2 
7. Se înlocuiesc dxi; ..., dXm în expresia lui d'L(P,) şi se aduce aceasta la forma 


canonică utilizând metoda lui Gauss sau Jacobi. 


8. Se determină natura punctului critic P, astfel, dacă: 


a) dL(P.) este pozitiv definită, punctul P, este punct de minim local 


condiționat, 


b) d'L(P,) este negativ definită, punctul P, este punct de extrem 


condiţionat, 
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ce) d'L(P,) nu păstrează semn constant, P; nu este punct de extrem 
condiționat. 


9. Se repetă etapele 5 - 8 pentru alte puncte critice: Pa,...,Pķ. 


Observaţie: Dacă în etapa 2 se determină mai multe sisteme de valori 
pentru multiplicatorii lui Lagrange. (şi, în consecință, pentru fiecare grup’ de 
-= multiplicatori se obțin alte puncte critice), etapele 3-9 se repetă pentru fiecare caz 


în parte. 


1.2.2: Exerciţii şi probleme propuse 
1. Să se determine punctele de extrem local condiționat pentru, următoarele funcții 
cu legăturile indicate: (m.l.c. (M.l.c.) — minim (maxim) local condiționat); 


: f(x,y) =x" +y? -x=y*] 
x+y=l ) ARE 


T'i 
R: à, =0,P, (—,—) =m.l.¢.; 
? 03 z) 


b) f(x,y) =X +y% -3x+2y +9: - 
x +y? =] 


4 DRE LE R | 
| de 2 siy. Aai 503 =m. .C. 
R: pa 
13 3 e | 

à> = -1+ —; P (=) — M.l.c. 
ii 2 j a13 VIES EAT 

e) f(x,y)=x+y +z" +11 
"az a at 23 


AR: M ==>, Po(1,1,1)—m.l;c. 
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f (x,y,z) = xyz 
„d)ix+y+z=5 
xy + XZ+ yz=8 


d2) A = 4, N, = -23 P2; 2,1), P302, 1, 2), PaT, 242) mul.c. 


„O fabrică produce 3 tipuri de piese. Profitul săptămânal fabricând x,y,z piese pe 


oră din fiecare produs este descris de funcția: 

fax, y,Z)= X+ 4y + 10z;+ 15 
Datorită faptului că săptămâna de lucru are 40 de ore, numărul de parse produse 
trebuie să satisfacă condiția (dictată de procesul de producție sau de clauze 
contractuale): i 

X+ y? +z =40. 
Determinaţi nivelul de producţie al fiecărui t en de piese în Gbaiië astie încât t profitu 


săptămânal să fie maxim. 


Ri = ahina bi pei za 5; fina 2.68. 
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